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PREFACE 


The Cryogenics Division of the National Bureau of Standards, with support from the 
American Gas Association, Inc., is engaged in a program to provide experimental data and 
computational methods for physical and thermodynamic properties of the constituents of 
liquefied natural gas mixtures (LNG). These thermophysical properties are important for 
all LNG technology. All operations such as liquefaction, separation, storage, pumping, 
transport, custodial transfer, and regasification will benefit from accurate data. As the 
commercial value of LNG depends on its heat of combustion, the densities of LNG mixtures 
and other requisite properties are especially important data for calculating the energy 


content of the liquid. 


The compositions of LNG mixtures vary widely, depending on the source and selective 
vaporization during handling. To predict properties for a wide range of mixture compositions, 
it is essential to know accurately the properties of the pure components and of selected 


binary mixtures of these components. 


The equation of state for the PVT properties of pure components is an essential tool 
for work on mixtures, because from this single formulation there can be obtained not only 
the density but also thermodynamic properties such as the enthalpy and specific heats at 


any desired temperature and pressure for which the component exists as a fluid. 


A major contribution from this laboratory for computations on the thermophysical 
properties of LNG and its components is the development of a simple, rational equation of 
state which originates on the liquid-vapor coexistence boundary, and gives a qualitatively 
correct behavior for derived specific heats, especially about the critical point. An 
earlier form of this equation was used in our report, "The Thermophysical Properties of 


Methane from 90 to 500 K at Pressures to 700 Bar,'' NBS Technical Note 653, April, 1974. 


In addition to the above publication, this laboratory has provided accurate experi- 
mental data on the densities, compressibilities, vapor pressures, melting pressures, 
specific heats, sound velocities, dielectric constants, refractive indices, and viscosities 


of compressed and liquefied methane at temperatures down to 90 K. 


The report, 'Provisional Values for the Thermodynamic Functions of Ethane" by R. D. 
Goodwin, NBSIR 74-398 (June 1, 1974) was prepared to discover the availability and 
deficiencies of experimental data needed for thermal computations, and to provide useful 
thermodynamic properties data at the earliest possible date. In the present report those 
results are modified by use of new experimental data, and by selection of a better 


computational path for entering compressed liquid states at subcritical temperatures. 


The present report makes use of many new physical properties measurements to obtain 
more accurate tables of thermodynamic functions, and provides the first of such tables 
available for liquid ethane below its normal boiling temperature (184.5 K). In this work 
we have used a simplified and more highly constrained version of the equation of state than 
formerly used for methane because the PVT data, taken from various sources, are not as 
self-consistent and precise as those for methane. Auxiliary interpolation functions also 
are improved, e.g., to accommodate the enormous range of ethane vapor pressures and 


saturated vapor densities. 
vi 


THERMOPHYSICAL PROPERTIES OF ETHANE, FROM 90 TO 600 K 
* 
AT PRESSURES TO 700 BAR 


R. D. Goodwin, H. M. Roder, and G. C. Straty 


Cryogenics Division, Institute for Basic Standards, 
National Bureau of Standards, Boulder, Colorado 80302 


The thermophysical properties of ethane are tabulated at integral 
temperatures over the entire range of fluid states at temperatures 
from 90 to 600 K along isobars to 700 bar. A new form of the 
equation of state is employed for the first time. Thermodynamic 
functions in the compressed liquid at T < T. are obtained by use of 
specific heats Cc (T) along a high-density isochore, and Cc, (7) for 
the saturated liquid. Use of these new specific heat data and also of 
new PVT data provide the first reliable set of thermodynamic functions 
for compressed liquid ethane at temperatures below its normal boiling 


point. 


Key words: Densities; enthalpies; entropies; equation of state; ethane; 
internal energies; isobars; isochores; isotherms; Joule-Thomson inversion; 
latent heats of vaporization; melting line; orthobaric densities; specific 


heats; speeds of sound; vapor pressures. 


1. INTRODUCTION 


The economic importance of ethane, as one of the major components of liquefied natural 
gas (LNG) mixtures, is well known. Our objective is to produce basic thermodynamic informa- 
tion, needed for the prediction of the thermophysical properties of liquefied natural gas 
mixtures. For the wide range of compositions encountered, it will be necessary to utilize 


accurate thermodynamic properties of the pure components. 


In [30] we computed tables of provisional values for the thermodynamic functions based 
on available specific heats and PVT data up to about 1972. To enter compressed liquid 
states at low temperatures we used a path of specific heats o along a high-pressure 


isobar [24]. 


This work was carried out at the National Bureau of Standards under sponsorship of the 


American Gas Association, Inc., 1515 Wilson Boulevard, Arlington, VA 22209. 


Nearly all of the analytical formulations of physical properties data in this report 
are new as compared with [30]. A major contribution of present work is the further develop- 
ment of a rational equation of state. As in [31], our description of the P(9,T) surface 
originates on a given liquid-vapor coexistence boundary (vapor-pressure and orthobaric 
densities equations). It yields a maximum in the specific heats C(p,T) at the critical 
point, and has only five arbitrary coefficients to be found by least squares from 


experimental P-9-T data. 


New experimental data from this laboratory, used in present work, include the melting 
line [78];\ vapor-pressures and PVT data [79]; specific heats, both for the saturated liquid 
and for single phases [66]; and speeds of sound, both for the saturated liquid and for 
single phases [86]. New densities of the saturated liquid and vapor from [37, 75] and from 


other laboratories also have been included. 


In present work we obtain thermofunctions in compressed liquid states by computation 
around the critical point (see figure 1) using new specific heat data from this laboratory 
[66] for C(t) along the isochore at 14.5 mol/%, and C,(T) for the saturated liquid at 
temperatures below ae of figure 1. Symbols and units of this report, listed in Appendix A, 


are similar to those for methane [29]. 
2. PHYSICAL PROPERTIES AND THEIR REPRESENTATION 
2.1 Fixed-Point Values 


The values selected are presented in Appendix B. 


(a) The Triple Point 


The triple-point temperature for ethane is from recent measurements by Straty and 
Tsumura [78]. The pressure is obtained from our vapor-pressure equation (2). The liquid 
density was assigned for consistency with our fitting function (3-a) for saturated liquid 
densities. The yapor density is obtained from our fitting function (3-b) for saturated 


vapor densities. 
(b) The Critical Point 


The critical-point temperature is from Douslin and Harrison [18]. The pressure 

is given by our vapor-pressure equation (2). The density was selected tediously by trial, 
based on detailed examinations of the critical isotherm and the rectilinear diameter very 
near the critical point. Our selected critical density of 6.80 mol/2 is lower than the 
value 6.86 mol/%, obtained via the rectilinear diameter analysis of orthobaric liquid and 
vapor densities of Sliwinski. Other reported critical densities in mol/2 are: Douslin and 
Harrison, 6.87; Strumpf, et al., 6.83 + 0.07; Miniovich and Sorina, 6.79 + 0.02. Older 
critical densities reviewed in [18, 21, 83], however, range from 6.75 + .07 to 7.32 mol/%, 


a spread of almost ten percent. 


(c) The Boiling Point 
At the normal boiling point (1 atm = 1.01325 bar), the temperature in Appendix B is 
from our vapor-pressure equation (2). The liquid and vapor densities are from the orthobaric 
densities equations (3-a), and (3-b). 
2.2 Melting Line and Vapor Pressures 
The melting pressures of ethane are from Straty and Tsumura [78]. In the Simon 


equation, 

Pop ho (GE) = i, (1) 
their parameters are qT, = 90.348 K, A = 2559.65 bar, and € = 2.179. The exceedingly small 
value which we use for Pe is from Appendix B. 


The vapor-pressure equation (2) for ethane is a modification of our earlier form [30]. 


Define the arguments x(T) = (1-7, /T)/(1-T,/T.)> and u(T) = (IST (Tats when the equation 


is 


Rn (P) a+ bex + ctu t+ alow + om + ome Cia), (2) 


Coefficients of (2) used in present work (for P in bar) are € = 1.30 and 


5.42844 31006 
-1.36232 70362 
0.76924 92586 


-11.38996 24306 d 
18.84522 82876 
- 7.63541 51345 


c 
Table 1 compares experimental data with calculated pressures from (2). 
2.3 The Orthobarie Densities 


(a) The saturated liquid density function is constrained to both the triple- and 


eritical-points. Define the variables-- 


iH 


x(t) = Cera ACERS Ls 
y(p) = (0-0 )/(P,-P.)s 
Y(p,T) = (y-x)/(x°-x), 


when the equation is-- 


Y=at bom © + Cex, (3-a) 


0.2818 66182 
-0.2899 37306 


€ = 1/3 b 
0.7196 84501 c 


a 
Table 2 compares experimental data with calculated results from (3-a). 


(b) The saturated vapor density function is constrained only at the critical point. 


6 


We use an estimated value for the triple point density, d, = 1.52°10 mol/%, merely to 


k 
obtain simple values for the coefficients. Define two arguments-—- 


3 


x(n) = (T-T)/(2-T,)> 


w(t) = (1,/T-1)/ (2/7, -1) 


and the dependent variable-- 


¥(0) = £n(d_/d)/2n(a/4,), 


when the equation is-- 


6 : 
Y=AsutAsx + DO pose = (3-b) 
1 2 i 
i=3 
with € = 0.35, and-- 
A, = 1.2023 48669 A, = -0.9803 17935 
A, = 0.1100 05895 A, = 1.3387 56298 
A, = 0.1371 69205 A, = -0.8077 62749 


Table 3 compares experimental data with calculated results from (3-b). 


2.4 The Virial Equation of State 


The following are taken directly from [30]. As seen below, the critical parameters 


Tos dq. differ from values in Appendix B. For the virial equation, these are merely 
arbitrary reducing constants which do not affect accuracy of the equation. They are retained 
from the original work [30] because in present work we repeatedly varied these parameters in 
developing the equation of state, and it is quite impractical repeatedly to readjust the 


virial equation. For the virial equation of state, 


Py/RE = 1+ B(T)*o + C1) “ov +. . (4) 
the second and third coefficients B(T), C(T) are dimensionless. We reduce temperature and 
density at the critical point, x = T/T; G = d/d.. Following initial research on the 
representation of these coefficients, we adopted McGlashan and Potter's formula for B(T), 
[53]; 


2 4.5 
B(T) + By + B,/x + Ba/x + B,/x af (4-a) 
where, 
~  Sbost d. = 6.74 mol/2 
Ose) Gr, Ba = -0.592 947, 
B, = -1.106 244, By = =()0VAIL GVA. 


Data from McGlashan and Potter were increased by 0.5 percent in absolute value to improve 
consistency with the data of Michels et al. [54], and of Douslin and Harrison [18] near 


300 K. Table 4 compares data and calculated values. 


For C(T) our new representation from [30] is 


3 5 
= 6 (ie 4-b 
C(T) [C,/x + C,/x + C,/x ea Ds (4-b) 
where, 
B = 305.37 K, d. = 6.74 mol/2% 
ie = Pilots Ko Cc, = 0.83253, 
Cc) = 0.24423, C3 = 0.53488. 


Table 5 compares data and calculated values. 


The virial equation of state usually is accurate at densities up to about one third of 
the critical density. In present work it has been used to obtain saturated vapor densities 
at low temperatures by intersection with our vapor-pressure equation. It also has been used 
to synthesize PVT data along an isochore at 0.4 mol/%, for use as input data in adjusting 


the equation of state. 


2.5 The Equation of State for Single Phases 


The PYT data used to adjust coefficients of the equation of state are summarized in 
table 6. From the data of A. K. Pal, reported by Pope [61], we selected only 13 isochores 
at the highest densities, runs Nos. 12 through 24 of table 12, because they were found to be 
self-consistent. We at first increased all of his densities by 0.5 percent to make these 
PVT data more compatible with data of other authors (table 6) via our equation of state. 

We finally used adjusted data of A. K. Pal kindly communicated privately by Professor R. 
Kobayashi, J. F. Ely et al. of the Chemical Engineering Department, Rice University. We 
omitted the data of Reamer et al. and of Tomlinson because those of Douslin and Harrison 


were believed to be more accurate. 


Background on this equation of state occurs in our work on methane and ethane [29, 30], 
and in summary reports [31, 32]. As described previously, this equation originates on a 
given, liquid-vapor coexistence boundary, and ensures a maximum in the specific heats, 


C(9,T), at the critical point. 


For a given density (isochore) we obtain the coexistence temperature T, (0) by 
iteration from the orthobaric densities, equations (3-a) and (3-b). Placing this in 
equation (2) yields the vapor pressures Plt (P) J as a function of density. The equation 


of state used here for ethane is described in [32], 
2 
P= P(e) + OR*[T-T,(0)] + O RT °F(0,T), (5) 
where 
F(p,T) = B(p)* &(p,T) + C(p)* ¥(p,T). 


Densities 0 = d/d. are reduced at the critical point, and the gas constant is defined for 


dimensional consistency as R = (9.0831434)*(d_), bar/K. 


The first of the temperature-dependent functions in (5) is 


1/2 


Oj) ZB x * &n(u), (5-a) 


where 


x = T/T.» W(t) = T/T, (0). 


The value of 6(9,T), and of its second derivative, 3°6/aT-, both are zero everywhere on the 


coexistence boundary at T = T (9). 


The second temperature-dependent function in (5) is ¥(p,T). It is nonanalytic, designed 


to give a maximum in specific heats at the critical point. 
We first define the arguments, 
w(p,T) = 1-O0(9)/T, 
w,(e) = 1-0(9)/T,(e), 
where 0(9) is our locus of temperatures inside the coexistence envelope of figure l, 
O(p) = T,(e)-exp[-a°f£(p)], 


£(p) 


3 3 
|o-1| @ ib) 
oO, = di/d. is a constant, and d. is density of liquid at the triple point. 


The function ¥(p,T) now is defined as the difference, 
¥(p,T) = p(p,T) - b Co), (5-b) 
such that ¥ = 0 at coexistence. Any component function b, 6) at coexistence is obtained 
from any function W(0,T) merely by replacing T with T,(P)- 
The component function selected for ethane after tedious investigations is-- 
WO) = Ofes se Glo) o|=w. <> oka) I 5 


where 0 <6 <1. Sufficiently far away from the critical point, Y behaves roughly like 


1/T, similar to the term 1/7? found in the well-known, Beattie-Bridgeman equation of state. 
The coefficients in (5) were developed for ethane-- 


= ee 3 
B(P) = B, + Bo + Bao, (5-c) 


nig at EF Ic (5-4) 


G@) a= (o-1)*(0-C.)*[C, + C5 


The parameters of equation of state (5) for ethane are-- 


fo) 

BL = 0.4875 2227 313 C = -0.4311 3918 548 
5 = 0.3319 8750 982 Cc, = 0.0 

B3 = 0.0685 4249 828 C3 = 0.0377 9460 468 


Table 7 shows behavior of the coefficients, B(p), C(p). 
6 


In adjusting the equation of state to PVT data we have obtained a well-behaved critical 
isotherm as follows. We constrain the slope (dP/0T) of the critical isochore to equal the 
slope (dP ,/dT) of the vapor pressure equation at the critical point, via the least squares 
program of R. D. McCarty [51]. We then vary the slope (dP ,/dT) of the vapor pressure 
equation at the critical point, by adjusting the exponent, €, in equation (2), to obtain a 
best "fit" of the experimental PVT data. Adjusting exponent, €, in equation (2) has no 
significant effect upon the representation of experimental vapor pressure data in table l. 
We finally explore the effect of small changes in the assigned critical density upon 
detailed behavior of the critical isotherm very near the critical point, in order to 


eliminate negative slopes of this isotherm, as demanded by equations (14) and (15). 


2.6 The Ideal Gas Thermodynamic Functions 
For use in our computations we have represented the internal energies of Chao et al. 


[12], by the following empirical power series 


fo} ) (443) /3 


9 
° as e 
Gs BPR So SO & Xu A, + (1/100 (6) 


The specific heats are obtained by differentiation, 
co(T) = dE°/arT, 


and the entropies by integration, 


As° = |/cPoay. 
P 


The coefficients for (6) are, 


A, = 21.705 718 A. = 906.218 4427 
A, = 65.498 641 Ag = -459.230 2545 
Ay = -362.011 5914 A, = 143.030 0226 
A, = 853.340 8616 Ag = -25.074 95605 
A, = -1123.601 6220 Ay = 1.897 540044 


The coefficient A, is the dimensionless integration constant for S°/R. The concise compu- 
tations are given by computer subroutine IDEAL in Appendix D. Table 8 gives comparisons 


with the data, and table 9 gives values interpolated by means of (6). 


We have compared the results of Chao et al. with those obtained earlier by Ziegler 


et al., (values in Joules, moles, kelvins), 


(0) oO 
T = 200 K We ge ce 
Ziegler et al. [93] FIBM A DON TZ. 42045 
Ghee ee ail, [2). 7257.6 210.50 42.26 


We see a difference of 24 J/mol in (HO-HP), or 0.3 percent. We have not further compared 
these independent calculations. Bier et al. [4] have found a change in the slope of c(t) 
beginning near 325 K, with values of c(t) as much as 1.6 percent smaller than from Chao 


et al. [12] at higher temperatures. 


2.7 The Heats of Vaporization 


Experimental data of Wiebe [89], Douslin [18], and Furtado [24] have been represented 


as follows, 


2 i/3 
Onna YA Yx >, kJ/mol, (7) 
i=1 


1 


where the argument is x(T) (T.-T)/ (1 -T,)> and 


T, = 90.348 K T. = 305.33 K 

A, = 12.857 784 A, = -97.520 486 
A, = 5.409 5302 A, = 100.849 76 
A, = 391095) ODI Ag = —37/ 5933) S73 


Comparisons with the data are given in table 10. 


2.8 Specific Heats c(t) at d. = 14.5 mol/2% 
As described in section 3.3, below, these data are used to obtain thermofunctions in 


compressed liquid states by computation around the critical point, figure l. 


We adjusted data of three experimental C(0,T) isochores from [66] (runs Nos. 1, 8, 9 


of table 14) to the assigned density of d, = 14.5 mol/% by means of the present equation of 


b 
state, and then represented these data with a standard deviation of 0.093 J/mol/K by use of 


the argument x = T/T. where qT. = 305.33 K, 


Cc (T,d,) =a/x+bt+cecex+t+ sox”, J/mo1/K, (8) 
= 64.02530 = 205.55275 
= -163.36173 d = -54.77598 


2.9 Specific Heats for Saturated Liquid 


As described in section 3.3, below, these experimental data are used to obtain thermo- 
functions in compressed liquid states at T < TS (see figure 1). As reported in [66], they 


were described as follows, with a standard deviation of 0.315 J/mol/K at qT, << Tos 


c,(@) =atbeT + cef(T) + d/T + e/T*, J/mo1/K, (9) 


where 
£(T) = mG =) °°. 
Le = 3K0)5)5 3}8} IK = 1.0906 40627 
= 22.6482 2683 d = 5371.2010 54 
= 0.0879 6160 711 e = -212338.994 


3. COMPUTATIONAL METHODS 
The numerical values for E and H in this report are on the same absolute basis as those 
of Tester [83], obtained by use of E. = (4.1868)+(4827.2) J/mol. 
3.1 The Homogeneous Domain 


The "homogeneous domain" of figure 1 includes all regions which can be attained along 
isotherms starting at zero density without crossing the vapor-liquid ''dome" and without 


passing close to the critical point at T > TO: 


We start our computations with ideal gas thermodynamic functions at zero density, and 


then integrate along isotherms by use of the equation of state in the following relations, 


AE =f (ra (ar/a0)}-40/0”, (10) 
ig = i - f (22/917) -a0/0" (11) 
ds = R-tn[P°/(prt)] + f [R-(aP/21)/0]-40/0. (12) 
Equation (12) is for use with initial entropies in hypothetical ideal gas states at 
fo) 
P 


= 1 atm. For the compressed liquid at T < Ty and 0 > Po» we use 


Nove —_f car/21)-40/p?. (12a) 


In each (9,T) state, reached by above integrations, we compute 


H = E + Pev, (13) 
c= C+ T+ (9P/9T)~/(8P/39)/p- (14) 
we C+ (aP/9p)/¢, (15) 


3.2 The Vapor-Liquid Transition 


As discussed below, we have used this computation only as a check against experimental 


heats of vaporization, and for closed-loop checks terminating on the saturated liquid. 


We traverse the vapor-to-liquid 'dome" of figure 1 by use of the Clapeyron equation, 


and Av = ys 


AH = T+ (dP_/dT) -Av, (16) 
AE = AH - PoAv, (17) 
AS = AH/T, (18) 


where (dP ,/dT) is slope of the vapor-—pressure curve. 


3.3 Compressed Liquid States 


As indicated by figure 1, we compute around the critical point into compressed liquid 


Sialtesiaten< Ty and 0 > Po as follows. 


(1) On isotherm Ty = 330 K use equations (10) and (12) up to qd = 14.5 mol/& to 
obtain E, and Sy of figure l, 


E(T, »d,) = 21,321.64 J/mol, S(T, 4d, ) = 165.4420 J/mol/K. 


(2) Along isochore qd. = 14.5 mol/%, obtain E and S at any temperature qT, ‘SS Ty by 
use of equation (8) for Ci (T,4,)5 


it Te 
AE = J @oche, AS = je °dT/T, 
Vv V 
uy i 


and finally, for E. and S. of figure 1 at T= ie = 256.5624 K, 


E(T,>d,) = 17,628.86 J/mol, S(T. »d,) = 152.8227 J/mol/K. 


(3) Along the saturated liquid path obtain E and S at any temperature qT, Soles qT, by 
use of equation (9) for C,(T), 


ab ae 2 
ING = fe Cae dT + = ff Po dp, /Po r) 

T T 

a a 

at 
AS = f C,°dT/T. 

at 

a 


For AE, above, Po is the vapor pressure, and Po is density of the saturated liquid. The 
specific heat C(0,T) at the saturated liquid boundary is computed by use of the equation 


of state in the thermodynamic relation, 


2 
C(0,T) = C,(T) + T+(8P/8T)+(doy/AT)/0,~- 


(4) The final path of computation is along isotherms, starting with values for E, S, 
and C_, either on isochore qd. at qT. < T,> or on the saturated liquid boundary at 


Vv 
i. Sn < Ts by use of equations (10) through (15). 


4. TESTS AND COMPARISONS OF CALCULATED AND MEASURED PROPERTIES 


4.1 Calculated Critical Isotherm, P(p) 


Table 11 gives a high-resolution examination of the P(~) critical isotherm from 
equation of state (5). It shows that the calculated slopes, dP/d9, are non-negative at all 


densities, as required by equations (14) and (15). 
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4.2 The P-p-T Compressibility Data 


The last two columns of table 6 give the rms of relative density deviations, and the 


mean of relative pressure deviations, respectively, for each source of the P-p-T data. 


Deviations of experimental densities and pressures from the smooth P(0,T) surface of 
the equation of state (5) are given in table 12. The data of Michels et al. [54] and of 
Douslin and Harrison [18] are highly precise, and the deviations generally are systematic, 
as might be expected from an equation of state with as little freedom as (5). Any 
inaccuracies in the given liquid-vapor boundary will be propagated along calculated isochores 
because the equation of state originates on this boundary. At high densities in the compressed 
liquid, the derivative 3P/dp becomes extremely large, hence pressure deviations should be 
ignored in this region. Typical deviations of densities along isochores are plotted in 


figure 2. 


4.3 Heats of Vaporization and Closure Computation 


The last column of page 1 of table 23 gives experimental heats of vaporization 
interpolated from equation (7) for comparison with values in column Q,VAP computed by the 


Clapeyron equation. 


Table 13 gives loop-closure computations for the saturated liquid. Values fcr 
enthalpy in column HLR are obtained by computing around the critical point, whereas values 


in column HLC are via the Clapeyron equation (see secs. 3.2, and 3.3). 


Heats of vaporization obtained via the Clapeyron equation are uncertain by several 
percent (up to about 500 J/mol) due mainly to uncertainty in the vapor-pressures and of the 
vapor-pressure equation, and to the vapor densities at low temperatures. In table 10, 
moreover, we see that experimental heats of vaporization may differ by over 2 percent at 
100 K, (about 350 J/mol). The deviations of about 0.5 percent (60 J/mol) seen in table 13 
therefore do not exceed known uncertainties. The heats of vaporization have not been used 


as input data for present computations. 
4.4 Specific Heats C(e,T) and C (9,T) 
ge NA Se ee eee 


Tables 14, 15 and 16 compare calculated specific heats with data of [66, 5, 24] 
respectively. As seen in tables 15 and 16, experimental specific heats oe) extrapo- 
lated to zero pressure differ by as much as + 1 percent from the ideal gas values c(t) 
used for our computations [equation (6) and table 8]. Excepting the critical region, at 
temperatures 305 + 20 K, therefore, the maximum deviations of about two percent indicate 
satisfactory agreement between calculated and experimental results. Typical deviations of 


C. along isochores are plotted in figure 3. 
4.5 Speeds of Sound, WT)» and W(0,T) 


Table 17 compares calculated speeds of sound for the saturated liquid with data from 


IGOR S65) 8i/)]\. 


iLL 


The deviations of about one percent in mid-range are quite satisfactory, considering 
the computational complexity, (section 3.3 and figure 1). At T < 256.562 K our computations 
start with smoothed experimental C,(T) data, and the deviations increased toward two 
percent as T > T- At T > 256.562 K we are entering the critical region, and the deviations 


increase toward ten percent as T > To: 


Tables 18 and 19 compare speeds of sound W(p,T) with data of [82,86] respectively. 
Except for the compressed liquid at very low temperatures, and for the critical region, we 
consider the deviations of less than two percent to be quite satisfactory. Typical deviations 
of the speed of sound along isotherms are plotted in figure 4. Approaching the critical 


point, deviations for the saturated liquid approach ten percent. 
5. TABLES OF PHYSICAL AND THERMODYNAMIC PROPERTIES 
5.1 Calculated P-p-T Isochores and Isotherms 


A selection of calculated isochores and isotherms is given by tables 20 and 21. They 
are useful to examine behavior of the surface generated by the equation of state, and to 


supplement the isobars of table 24 in obtaining P-p-T values and their derivatives. 


5.2 ‘The Joule-Thomson Inversion Locus 


Table 22 gives our calculated P-9-T locus for the Joule-Thomson Inversion, (0T/9P),, = 0. 


These results are obtained from the equation of state under the condition, T*(dP/0T) = p*°(d0P/dp). 


5.3 Thermophysical Properties of the Saturated Liquid 


Table 23 gives physical and thermodynamic properties for the saturated liquid. Column 


headings are interpreted on the first page of this table. 


5.4 Thermophysical Properties Along Selected Isobars 


Table 24 gives physical and thermodynamic properties on isobars, as computed by methods 
of section 3. Explanations for the table are given on the first page. This table is 
extrapolated beyond the pressure range (v 350 bar) of P-0-T data used for adjusting the 


equation of state. 
5.5 Tables in British Engineering Units 


Thermophysical properties in British engineering units are given in tables 25, 26, and 
27. Temperatures are given at integral values on the Fahrenheit scale, and presures at 
integral values in pounds per square inch (psi) absolute. Conversions of units in 
Appendix C are similar to those in [29] except for the greater molecular weight of ethane 


as compared with methane. 
6. ESTIMATED UNCERTAINTIES IN THE THERMODYNAMIC PROPERTIES TABLES 


The problems of estimating uncertainties in a thermodynamic network (table 24), 
affecting the relative differences between tabulated values, have been discussed by 


Kuloor et al. [46]. 
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A comment must be made concerning inherent difficulties near the critical point, which 
are independent of any equation of state. The derivatives d0/dT and d0/dP become infinite 
at the critical point. For this reason, experimental densities from different laboratories 
may disagree radically at any (P,T) point in the critical region, due merely to slight 
inconsistencies in the pressure and temperature scales used. At the other extreme, the 
derivative 0P/d9 becomes very large for liquid at low temperatures. The slightest experi- 
mental error in density will yield (via any equation of state) deviations of the calculated 
pressure greatly exceeding the accuracy of measurement. In this domain any equation of 


state should be used only to find p(P,T). 


(1) Vapor pressures of table 1 from four different sources at P > 1 bar agree to less 
than 0.05 percent on average, which we take to be the absolute accuracy. At lower pressures, 
approaching the triple point, the data estimated by.Ziegler et al. [93] may be uncertain by 


up to two percent. 


(2) For the saturated liquid densities in table 2 we believe the most reliable data 
to be those of Sliwinski [75] and of Haynes and Hiza [37]. The former estimates experimental 
uncertainty of 0.08 percent near the critical point, diminishing to 0.01 percent at 283.2 K. 
The latter estimate overall experimental uncertainty of 0.1 percent near the triple point, 
diminishing to 0.06 percent at 300 K. The results from table 2, therefore, should have an 


absolute uncertainty no greater than 0.1 percent. 


(3) The saturated vapor densities of Sliwinski [75] in table 3 have an estimated 
experimental uncertainty of 0.1 percent near the critical point, increasing to 0.36 percent 
at 283.2 K. Some of the data from Douslin and Harrison [18] in this range are several 
tenths of one percent higher than those of Sliwinski. Below 280 K we calculated these 
densities from our virial and vapor-pressure equations, estimating the uncertainty to 
approach two percent at the triple point. In summary, uncertainty of the saturated vapor 
densities increases from about 0.2 percent near the critical point to about 2.0 percent 


near the triple point. 


(4) The selected reduced second virial coefficients in table 4 are seen to have been 
represented with a mean deviation of 0.14 percent. The estimated (not experimental) data 
of McGlashan and Potter [53], however, were increased by 0.5 percent in absolute value to 
improve consistency with other data. The overall uncertainty of table 4 may be under 0.5 


percent. 


The selected reduced third virial coefficients in table 5 have been represented 
with a mean deviation of about one percent. Original publications should be consulted in 


each case for estimates of experimental uncertainties in the data for virial coefficients. 


(5) The highly constrained equation of state (5) used in present work yields 
derivatives which are qualitatively consistent with the known behavior of isochoric 
specific heats, C(p,T), including a maximum at the critical point. Because this equation 
originates on the orthobaric, coexistence boundary, it gives deviations of the density up to 


two percent in the critical region, figure 2. Uncertainty of the densities from this 


iS) 


equation in table 24 therefore approach two percent in the critical region, a + 102, 
ae 50%, but generally are less than 0.2% elsewhere for the data used in adjusting the 


equation of state, as seen by inspecting table 12. 


(6) The ideal gas specific heats c(t) recently measured by Bier et al. [4] differ by 
up to 1.6% from the spectroscopic calculations of Chao et al. [12] at temperatures above 
about 325 K. The overall effect upon E°(T) and S°(T) will be less due to integration. We 
nevertheless should for this reason set an upper bound of about two percent upon uncertainty 


of the important enthalpy differences in table 24. 


(7) Experimental specific heats C,(T) and C (0,7), used in present computations, were 
obtained by Roder [66] who observed pressure and temperature, and then computed density of 
the sample via an equation of state. Uncertainty due to this estimated density may approach 
two percent near the critical point. As seen in figure 3, deviations from our calculated 


results are greater than this near the critical point. 


The use of specific heats from [66] to enter compressed liquid states at p > pe 


and at T < T, = 330 K we believe increases uncertainty of enthalpies in this region to 


b 
about 4%. 


(8) Table 13 shows an inconsistency (amounting to about 0.5 percent of the heat of 
vaporization) for enthalpies of the saturated liquid computed via the Clapeyron equation, 
as compared with computation around the critical point. A further difficulty is that the 
specific heat C (0,7) at point (E,, S,) of figure 1 is 56.11 J/mol/K as calculated from the 
ideal gas specific heat and equation of state along isotherm Ty > whereas the experimental 
value at this point is 54.05 J/mol/K via equation (8) and the data of [66]. We conclude 


that the accuracy of specific heats in table 24 is four percent. 


(9) In summary, the uncertainty of enthalpies is four percent in compressed liquid 
domains at low temperature, increasing up to ten percent very close to the critical point, 
and diminishing toward two percent at low densities and for all densities at temperatures 
well above critical. The specific heats are uncertain by four percent. The speeds of 
sound are accurate to one percent, figure 4, except for the saturated liquid, where 


deviations approach ten percent near the critical point, table 17. 
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APPENDIX A. Symbols and Units 


Subscripts c and t refer to critical and liquid triple points. 
Subscripts g and & refer to saturated vapor and liquid. 
Subscript 0 refers to liquid-vapor coexistence (usually the liquid). 


Superscript o refers to ideal gas states. 


Oa Wy non-linear constants in the equation of state 

B(p), C(p) density-dependent coefficients in the equation of state 
BCE) 5 (Ge) second and third virial coefficients 

C(P,T), molal heat capacity at constant volume, (J/mol)/K 
ee molal heat capacity at constant pressure, (J/mol)/K 
C,(T) molal heat capacity for saturated liquid, (J/mol)/K 
d, density, mol/2% 

dis 1.52°10° mol/% for eq. (3-b) 

E(9,T), the internal energy, J/mol 

ED (4.1868)° (4827.2) J/mol, [83] 

£(0) used in the definition of O(p) 

1 (oma 5 definition in the equation of state 

(Os) 5 the enthalpy, J/mol 

J; the joule, 1 N-n, 

ue the liter, me a 

mol, 30.07 grams of ethane (cl = 12 scale) 

De pressure in bars, 1 bar = 10°N/m2 (1 atm = 1.01325 bar) 
Po?) the vapor pressure, bar 

P60) > Plt (e)l, vapor pressure as function of density 
O(0,T), function in the equation of state 

VAG) suh)) 4 function in the equation of state 

Orap? et ee the heat of vaporization, J/mol 

R, the gas constant, 8.31434 (J/mol)/K, 0.0831434 (bar-%/mol) /K 
Og d/d., density reduced at the liquid triple point 

0, d/d_, density reduced at the critical point 

S(p,T) the entropy, (J/mol)/K 

Te temperature, K, (IPTS-68), [68] 

Ts 256.5624 K, (figure 1) 

Th? 330.0000 K, (figure 1) 

T, (0), liquid-vapor coexistence temperature, K 

O(p), defined locus of temperatures 

u(o,T), T/T (0) for eq. (5-a) 

Ws 1/d, molal volume, 2/mol 

w(p,T), [1-0(p)/T] for the equation of state 

W(p,T), the speed of sound, meters/second 

se) 5 T/T, for the equation of state 
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APPENDIX A (continued) 


Gh) variously defined for other equations 
X, (0) > T,(e)/T > reduced temperature at coexistence, for the equation of 
state. 
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APPENDIX B. 


Fixed-Point Values for Ethane 


scale) 


Triple Point Boiling Point Critical Point 

Temperature, K. 90.348 184.574 305.33 
Pressure, bar 1.130810 > 1.01325 48.714 
Density, mol/% -- -- 6.80 

Vapor iL SIS4ong? 0.068381 --- 

Liquid 21.680 18.094 -- 

APPENDIX C. Conversion of Units 

Temperature, 1.8 Rankine (R) = 1 Kelvin, (K) 
Pressure, 14.69595 psi = 1 atm = 1.01325 bar 
Volume, 1 cu ft = 28.31685 liter, (2%) 
Mol. Wt., CoH.» 1 mol = 30.070 gram, (g), (Goa 
Mass, 1 1b (avoir.) = 453.5924 g 
Density, CoH, 1 1b/cu ft = 0.5327 0577 mol/2% 
Isochore slope, 1 PSI/R = 0.1241 0562 bar/K 
Isotherm slope, C, 6? 1 PSI cu ft/lb = 0.1294 290 bare%/mol 
Energy 1 BIU = 1054.350 joule, (J) 
Specific energy for ethane, 1 BIU/1b = 69.896 0223 J/mol 
Specific entropy and heat 

capacity of ethane, 1 (BTU/1b)/R = 125.812 840 (J/mol)/K 
Speed of sound, 1 ft/s = 0.30480 m/s 
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(PY (52) (>) (5) 


(>) 


Appendix D. Computer Programs for Thermofunctions. 


PROGRAM ETHRM7A (INPUT,OUTPUT) 
REVISION OF ETHANE THERMOFUNCTIONS, FEBe 1 THRU NOVe 20, 1976.6 
ENTER COMPR. LIQ. STATES VIA NBS(RODER) CV ON ISOCHORE RB» AND CSAT. 
NOTE OMISSION OF FINDSATF (TM), USE DENGASF(T), DENLIQF(T). 
NOTE THAT TSATFC(DEN) ITERATES ON DENGASF{T), DENLIQF(T). 
COMMON B25 825B855 Cis C2.CS5 ElsE2sess ER 
COMMON/1/GK,AL58ByBEsGA,DE,EP, DCRT,TCRIyPCRT ,DGAT,0TRP,»TTRP»yPTRP 
COMMON/ 3/0FOT »D2P0T2,DPSOT,DOPMDI ,DPDD,DPDR,0TSDR;DTHDR 
COMMON/4/ X81 5XB25 XCisXC25 XEL,XE23 OXBDR,OXCDOR,OXEDR 
COMMON/S5S/ ODSOT 
COMMON/Zo/ TSATs THETA, PSAT 
COMMON/7/ TAsTBsPBoRBo EAsHAsySAsy EBsHdySB 
COMMON/43/7 INoIKe PsTsDENeo EgHosSs CVsCPyCSATs WyWK 
COMMON/S/ AUsBU,CU,0Uy AS,BS,CS,D0S9E Ss ETBsSTBoSTA 
COMMON/1G/ NFE»s ETRP, EXK(9) 
COMMON/99/ TI,EZZy EZsSZ5CVZs5 HZ,CPZ 
DIMENSTON PP499) 
1 FORMAT(1I5, cFiG6.9) 
2 FORMAT(I5, 3Fi0.0) 
3 FORMAT(1Hi 114X *LOOP CLOSURE CHECK FOR SATURATED LIQUID,* / 
1 12X FENTHALFY HLR VIA NBS CV» HLC VIA CLAPEYRON EQNe>s J/MOL* 47 
2 12K3HT,K 7X3HHLR *XSHHLC 7X3HDIF &X5HQ, VAP 2X8HDH/Q,PCT ) 
FORMAT(5X 5I1G,5 F8e2) 
FORMAT (1X) 
FORMAT(I4, F7e05 F005 2F720) 
FORMAT(5X [55 SF1IG.Sy CFI0.2s FS. 2) 
FORMAT@LH1 9X FETHANE FUNCTIONS ON ISOCHORE RB = 14.5 MOL/L*¥ // 
140xX fHAT T3 =FBeSy OH» P3 =FBe35 TH» DEN =F7.23// 
140X GHEB =FiGe2s 6H, HB =Fide2s OHy SB =FI004// 
40X GHCOV =F 6e3e 16H,» CV»sXPT =F8.3/// 
4GX 7HAT TA =F9e4, OHe PS =F8es5 7Hs DEN =File37/ 
1OX4HEA =F10.25 6H HA =Fide2y 6Hs SA =FI0.4//10X 4HCV =F8-3/// 
13X SHETB =E17.9» Hs STS =E17.95 THs STA =E17.9) 
i2 FORMAT (6Fi0.%) 
44 FORMAT(1H1 18X *ETHANE ISOBAR AT P =* F8e59 4H BAR JS ) 
15 FORMAT(1H1 18X FETHANE ISOBAR AT P =* I4, 4H BAR 7 ) 
to FORMAT({i9X 1HT 6X3HDEN OX3HVOL SXSHDP/DT SXSHUP/OD 
2 &XLHE oX1LHH &XLHS BX2HCV GX2HCP SXLIHW / 
3 i15X SHDEG K &XSHMOL/ZL &XSHL/MOL SXSHBAR/K 1X9HBAR=L/MOL &X5HI/MOL 
& &2XSHISMOL 2X7HISMOL/K LX7HISMOL/K LX7HI/SMOL/K 1AXSHM/SEC ) 
417 FORMAT(10X Fibe3Sy FI0Ss FI05s FIDe4s FIidedSy 2F9019F9-392F8e2e156) 
43 FORMAT@10X FiGe3y F9e55 F903, FIDC6y FID0Ss 2FI9e1,F90392F8e2 1K) 
21 FORMAT(1H1 7X *FRODER ETHANE SPECIFIC HEATS CV(DsT)5 J/MOL/K * // 
1 8X2HID SXSHMOLYL 7X3HT,K SXSHP,BAR &XOHCV,XPTSXSHCALCDO4X4HPCNT ? 
Ze) FORMAT (Lhd 7x Seabee) EnHANE SPECHEEC RHE AIS = GRIGP aii isms /1M Ol 7 Kee amaan 
1 6X2HID SXSHP,BAR 7X3HT,K SXSHMOL/ZL BX2HCP SXSHCALCD 4X4HPCNT ) 


Coo aon - 


OWN FwWNHe 


CONSTANTS FUR EQNSTATE MARK VII, 4/12/76. 
23 WM = 30.07 $5 WK = 100000/WM $ Q = 1201325 $ GKK = Q.0831434 


24 TTRP = 96-2348 $ ODTRP = 212.68 $ PTRP = 16130842618E=-5 

2> TORT = 305-33 $$ OCRT = 6-60 $ PCRT = PSATF(TCRT) 

26 DGAT = 1.5153661E-6 $ GK = DCRT*GKK $58 EZZ = 401868*4827.2 
ey AL = 1 $3 DE = 2.0/3.0 $ ER=2 $ BB=BE=GA=EP=9 


28 B81 = 0248752227313 $ B2 = 0633198750982 $$ BS = 0406854249828 
29 EL =-0,2431159185453 $ E2 = 0.0 $ E€3 = 0.03779460468 
38 CALL PEEK 5 CALL ISOTHRM 
INTEGRATION CONSTANTS FOR CVSUMIT, SSATF(T). 
39 1382330 $ RB=14.5 $ TA=TSATF(RB) $ X=TB/TCRT $ G=ALOC(X) 
31 AU=64.0253 $ BU=-163.36173 $B CU=205.55275 $ DU=-54.77598 
32 ETS = AUFG + BUYX + CUFX*X/2 + DUEX*X*X7S 
B83 Sue) = CMW BWR sp CW & UMM ene 
34 AS=22.64822083 $ BS=0.08796160711 $ CS=1.090640627 
35 DS=5371.201054 $ ES=-212338.994 § T=TA $ TO=(TCRTI-7T)**0.4 
36°STA = AS*ALOG(T) + BS*¥T = CS*TD70.4 = DS/T = ES727771 
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Appendix D. Continued. 


C 


(Sp)i(S2) (Sp) 


(7 Ol Pk 2) 


ESTABLISH TBsTA»y EBsSB,y EAsSA, ON ISOCHORE RB = 14.5 MOL/L. 

33 T= TB $ P = PB = PVTF(T,RB,D) $$ CALL THRMOGS 

40 RC=DEN $ EB=E $ HB=H $ SB=S $ CVB=CV $8 CVX=CVXF(T) 

41 T = TA $ CALL CVSUMIT $ EA=-E $§ SA=S 

42 PSTA = PS = PSATF(TA) $ HA = EA + 100*PS/RB 

43 PRINT 105 T3sPBeRC,y E8sHBoSBy CVBsCVX>s TAsPSsDENs EAsHAsSA5CVy 
1 ETB, STB, STA $ CALL ESATFIT 


LOOP CLOSURE COMPARISON FOR SATURATED LIQUID. 
EOMP RUN One Se UO Se =I h2 0 > le = ht = 90 + 5) fy Fr = PSATE (Th) 
46 CALL SATGSTR $ OG = DEN $ HG = H $ OL = DENLIGF(T) 
47 Q = 100*T*OPSOT*(1/0G = 41/0L) $ mILG S tlG > el 
G8 IF(TelLTeTA) 49550 
49 CALL SATLQTR $ HLR = H §$ GO TO 51 
56 CALL THRMXL $ HLR =H 
51 DH=HLC-HLR $ QPCT=100*DH/Q $ IHLR=HLR $ IHLC=HLC $ IVJH=DH $ IQ=Q 
52 PRINT &, IT, IHLRyIHLC, IOH,IQ, QPCT 
55 CALL JTLOCUS $ CALL TABLIQ 
60 CALL WSATLIQ $ CALL TSUMURA $ CALL TERRES $ CALL FURTADO 


COMPARE NBS(CHMR) SPECIFIC HEATS, CV(R,T). 
70 PRINT 21 $ L=d $$ ODO 83 J=1,300 
71 READ 859 IDs PXoDXoTsCX $ IF(ID) 72,84 
72 P = PVTF(T 50X50) § fh S jlssl 8 wR Cboselo sl) 7sig 4 
vo (he OS ARSC et 
74 IFC(TeLT.TB) 75581 
75 IF(TLT.TA) 70,77 
76 CALL LIQTHRN $ GO TO 82 
77 IF(T.LT-TCRI) 78,86 
73 PX = PSATF(T) $ IFC(PeGTePX) 79,581 
73 CALL THRMOX $ GO TO 82 
80 PX = PVTF(T,DCRIT,0) $ IF(PeGTePX) 79,81 
81 CALL THRMOGS 
82 PCT = 100*(1-CV/CX) 
83 PRINT 99 10, DXefosPo CXoCV, PCT 
64 CONTINUE 


COMPARE SPECIFIC HEATSs CP(RsT).» 

BIERsKUNZEyMAURERs Jeo CHEMe THERMODYNe 8,5 857 (1976). 
SoePRING ee SE = 6 9s DO) 10/0 J=1,5/50/0 
86 READ 25 IDs TCoPoSP $ IJIFIIO) 87,101 
87 T = 273-15 + TC $ CX = WM¥CP $ L=Lti $$ IF(LeEQ.44) 88,89 
opeie= 0S PRINT 22 
89 IF(JoL7.10) 96,91 
Siete Ss) 6 6CALL TDEAL 6©$06CP=CPZ $3 DEN=0 $ GO TO 99 
91 IF(T.LT.T8) 92598 
92 IFiTeLTeTA) 93594 
93 CALL LIQTHRM $ 30 1T0 99 


94 IF(TLTeTORI) 95597 

95 PX = PSATF(T) $ IF (PeGTePX) 90,98 

96 CALL THRMOX §$ GO TO 99 

97 PX = PVTFI(T,DCRT,0) & IF {PeGTePX) 96598 
98 CALL THRMOGS 

99 PCT = 100.0*41-CP/CX) 


100 PRINT 95 IDs PsTsDENs CXsCP,y PCT 
161 CONTINUE 


COMPUTE THERMOFUNCTIONS ON ISOBARSe 
EACH ISOBAR STARTS ON THE MELTING LINE. 
ISOBARS AT FP UNDER PCRT TRAVERSE THE DOME. 


Poin = 1 S$ NI = 03) $ READ 125, (PRCT oT=1,NI) 
fae OO SOG L=EN,NT 3 WiKi. 1) S$) IP = P= PP(L) 
112 IFC(IoLT.5) 1135114 

143 PRINT 24, P $3 GO TO 125 
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Appendix D. 


(52). (SP) 


Continued. 


114 1(P) 
115 
120 
tal 
ILeZe 


PRINT 15,5 
PRINT 16 
T FINDOTMF(P) 
PRINT 17, 
IY 1/10 


3 
TsDEN 5 V9 
$ 


CALL LIQTHRM 
CPOT,OPOD, EoH 
IF(PeLTePCRT) 153052 


TWO MAIN CASES FOR P UNDER PCRT. 
130 TS IPsCNIOU SF UIP) 25 ed TS + 10 
AS We) 99 Jelo9) 5 ys 10*(1IT 
132 IF(P=PSTA) 13345133,150 


JT 


$ 
9Sy 
v0 


V 1/U0cN 


CV»CP,IW 


IW 


3 
+ J) 


MAIN CASE A FOR P NO GREATER THAN PSTAe 


133 
13+ 
£35 
ASI). 
137 
133 
139 
140 
141 
142 
143 


IF(TLTeTS) 
CALL LIQTHRM $B V=1i/0EN $8 IW=W 
PRINT 179 Ty»DENoVe ODPDT,OPDD, E>H 
IF(TeLToTX) 1375186 

T TS $ CALL SATLQTR 
PRINT 175 TsDENeVs DPOT,0PD0,% E,H 
CALL SATGSTR $ V=1i/DEN 3 IW=W 
IF(PelLT.20.0) 141,142 

PRINT 18,5 TsDENsVy DOPDT,DPOD, 
PRINT 175 T,DENyVy DPDT,DOPOD, 
T JT $ 30 TO i180 


134,136 


$ V=1/0 


EsH 
EH 


9559 CV »CP,IW GO TO 199 
EN $ IW=W 


959 CV,CP,IW PRINT 


999 
959 


CV,CP,IW GO TO 143 


CV,CP,IW 


MAIN CASE 38 FOR P GREATER THAN PSTAe 


MIP (Powis) abbals eal, Ss) 
CALL LIQTHRM $8 V=1/DEN $ IW=W 
PRINT i175 1T,DENsV» DPDT,O0POD, EsH 
IF(TeLT.TS) 154,156 
CALL THRMOX B V=1i/0EN $ IW=W 
PRINT 17, T,DENoV> DPDT»DPDD, 
IF(T»LToTX) 1575183 
T TS $ CALL TARMXL 
PRINT 17,5 Ts,DENsV, OPOT,DPDD, EysH 
CALL SATGSTR $ VW=1/0EN $8 IW=W 
160 IF(PeLT.20.5) 1615162 
PRINT 165 TsDENoVs DPOT,DPDD, 
PRINT 175 TyDENoV_y ODPOT,OPODD, 
ili JT 
CAScS FOR HOMOGENECUS DOMAIN, 
IF(JT.GI.5tu) 181,135 
K Kt DS JT 
IF(UToGTe70U) 300,183 
CALL HOMOTRH $ V=1/0EN $ 
IF(PeLTecued) 185,186 
PRINT 185 T»DENoV,y DPOT,DPOD, 
PRINT 179 T,O0ENoVsx OPDT,0POD, 
CONTINUc 


EsH 


$ 


E»sH 
E»H 


DEN 


JT + 10*K 
IW=W 


EqH 
EyH 


CASES FOR P > PCRT. 
K=L=0 $8 O69 230 J=1,99 
IF(TeLTeTA) 2025205 
CALL LIQTHRM $ V LT DEN SS 
PRINT 175 T»,DEN»yV,y OPOT,OPDD, 
IF(TeLTeTCRT) 211,206 
PX PVTF(T,OCRTgu) $ IF(P.GT.P 

CASE FOR REGION X AT TA < T < TBs 

c1ii CALL THRMOX 3 VY 1/D0EN $3 IW 

212 PRINT 175 TyDENeVs DPOT,OPDD, EsH 

CASES FOR HOMOGENEOUS DOMAIN, DEN 
IF(3T-500) 22592253221 
K Keitel . Skil JT 
IF(JT-700) 2255225,300 
CALL HOMOTRM 3 V 1/DEN $ 
PRINT 175 TsDENsV, DPDT,DPOD, 
CONTINUE 


200 S yee 
ratpe 
202 
203 
205 


Zuo 


IW 
EyH 


= 


224 
fail 
Gee 
ZED 
226 
230 


SIP sp abihslx 


= 


IW 
E 9H 


26 


V=1/DEN $ 


9S» CV,CP,IW GO TO 199 


959 CV »CP,IW GO TO 
IW=W 


9S» CV»,CP,IW 


939 163 


999 


CV,CP,IW 
CV,CP,IW 


£ DCRT, OR T 


959 $ GO TO 199 


9359 


CV,CP,IW 
CV¥»CP,IW 


10*(IT+y) $ IF(TeLT.TB) 


9S9 


W 


CV»CP,IW $ GO TO 230 
X) 211,225 

P ABOVE CRITICAL ISOCHORE. 
W 

9S9 CV,CP,IW 


<SOCR sO 


GO TO 230 
TBe 


$ 


> 


959 


W 
CV,CP,IW 


201,220 


Appendix D. Continued. 


300 CONTINUE 
999 STOP $ END 


FUNCTION CSATXF(T) 
C FORMULATION OF NBS{RODER) ETHANE DATA, 3/765 J/MOL/Ke 
C THE STANDARD DEVIATION IS 0.315. 
COMMON/9/ AUsBUsCU,DUy ASyBSeCS,DSyESy ETBesSTB9STA 
DATA (TC = 305.233) 
IF (7TeLT.TC) 299 
F = T/(TC-T) **0.6 
CSATXF = AS + BS*T + CS*F + DS/T + ES/T/T $$ RETURN 
CSATXF = 0 $ RETURN $ END 


OWN Fe 


FUNCTION CSUMF(NsT+,0A,0B) 


C DELTA CY = -T¥INTEGRAL ((02P/0T2)/D**2)*00. 


COMMON/3/DPDT ,D2P0T2» DPSDT,DPMDT » OPDD, DPOR»DTSDR,DTHOR 


2 DX = (O0B-0A)/N $ CSUMF = 0 $ DO 5 J=1y5N 
3 DN = DA + (J-005)*DX $ P = PVTF(T,0N,0) 
5 CSUMF = CSUMF = D2PD0T2*0X/0N/0DN 

9 CSUMF = 100*T*CSUMF $ RETURN $ END 


FUNCTION CUBERT(X) 


C CUBE ROOT MISSINS FROM 6600 COMPUTER LIBRARY. 


7a Pk @ a @) 


&€ = 1.0/3.0 $ IF(X) 259354 
CUBERT = -A3S(X)*F*E $ RETURN 
CUBERT 0.0 $ RETURN 

CUBERT [8 $ RETURN $ END 


Fane 


SUBROUTINE CVSUMIT 
GIVEN PeT» GET JEN AND SOME FUNCTIONS AT POINT (T5RB). 
START AT TB = 338 K WITH VALUES EB, SB, THEN, ON ISOCHORE RB = 
USE INTEGRALS OF CVXF(T) FOR DEL(E), DEL{S) FROM TB DOWN TO Te 
YIELDS ONLY DEN» EoSs CV. 
COMMON/7/ TAsTB»9PBsRBy EAsHAsSAsy EBsHB,SB 
COMMON/8/ INogIKsy Pol oDENsy EsHsSy CV5,CP,CSATs WyWK 
COMMON/9/ AUsBU,CUsDUs5 AS»BSeCS5DS5ES> ETBySTBsSTA 
DATA (TC = 335.33) 
1 xX = T/TC $ G = ALOG(X) 
ETT = AU*G + BUFK + CUFX*X/2 + DUFX*¥X*X/3 
—£ = —€B t+ TCFK{ETT - ETB) 
STT = -AU/X + BU*G + CUFX + DUFX*X/2 
S 2 Sh - Sup oc Sys 
DEN = RB §$ CV = CVXF{(T) $ RETURN $ END 


ODmavw 
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FUNCTION CVXF(T) 

ORMULATION OF NBS(RODER) ETHANE DATAy 3/765 =J/MOL/Ke 
VXF IS CV(T) ALONG THE ISOCHORE, % = 1465 MOL/L. 
ALID FROM TB = 330 K DOWN TO TA = 256.562 Ke 

TANDARD DEVIATION FOR 46 POINTS I 02093. 

COMMON/9/ AUsBU,CU,DU, ASyBS9CS5D .ES»y ETBsSTBsSTA 

X = 7/305.-33 $8 CVXF = AU/X + BU + CU*X + DU*X*X 
RETURN $ END 


FUNCTION DELTAF(T,D) 


GET (T*DOP/0T - D*DP/DD) FOR THE J-T INVERSION CURVE. 


OON OW FEF WHR 


COMMONS 3/0PUT ,D2PDT2 ,DPSOT,DPMDT,0P9D, DPDR,DTSOR,OTHOR 
OATA (DCRT=5.80) 5 (TCRT=305-233) 

IF (T=TCRT) 254 54 

DL = DENLIGF(T) $ IFCD=-DLI 353,54 

DELTAF = 1e6E+106 $$ RETURN 

P = PVTF(T,D51) 

DELTAF = ABS (T*¥OPDT-D*DPDD) $ RETURN $ END 


FUNCTION OENGASF (T) 
DEN = DORT*EXP(-YN*Y) 5 YN = LN(OCRT/DGAK) y 
X = (TCHT)/(TCHTT), Q = X¥¥1/3, Z = (TC/TH1)/S(TCSTT=-1) 5 
Y = A1¥Z + A2*XE + AS*FQ3 + AYFQK + « « o 

COMMON/5/ DOSDT 

DIMENSION AVI6) 

DATA {(NFG=6) 5 (EG=0.35) y (DGAK=1-52E=-6) 

DATA (TTRP=90. 348) »5(TCRT=305-33) » (0CRT =6~- 80) 

DATA(AV = 1.202348669, 0.110005895, 0.137169205, 

1 -0.980317935, 1.£338756298, -06807762749) 


FORMAT(LHG 9X *DENGASF = QO, T EXCEEDS TCRT. * / ) 

IF CTCRT=T) 39455 

ehh ob 85 SOP 

DENGASF = OCRT $ ODSOT = 1-0&+10 $ RETURN 

XN = TCRT-TTRP $$ ZN = TCRI/TTRP=1 3 YN = ALOG{DCRT/OGAK) 
X = (TOCRT-T)/XN $$ OXOT = -120/XN 

Z = (TORT/T-1)/ZN $ ODZDT = -TCRI/SZN/T/T 

Q = CUBERT(X) 5 Q1T = Q*OXDT/3/X 


XE = X**EG $ XE1T = EG*XE*DXDT/X 

Y = AV(1)*Z + AV(2)*XE $$ YiT = AV(1)*OZDT + AV(2)*XE1T 
DO 14 K=3,NFG $ QK = Q¥*K $ Y = Y + AV(K) *QK 

YiT = YiT + AV(K) *K*QK*Q1iT/Q 

DENGASF = DCRT*EXP(-YN*Y) $ ODSDT = -YN*Y1T*DENGASF 
RETURN $ END 
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FUNCTION DENLIQF(T) 
DEN = DCRT + YN*(X # (XE=X)*Y). 
VS PALE ACK S(T) 4 TAs eX, YN = DTRP=-DCRT. 

COMMON/5/ DOSDT 

OIMENSION AW(3) 

DATA (TTRP=90 6348) 5 (TCRT=305033) y (OCRT=60 80) y (OTRP=21.68) 
DATA(AW = 0.719684501,5 0.2818660182, -0.289937306) 


PePORMATCMAO 9X *OENETQR = 0, T EXCEEOS TGRI. *-7 ) 
2 EL = 1-6/3.) $ IF(TCRT=T) 354,55 
3 PRINT 12 $ STOP 
4 DENLIQF = OSRT $ OOSOT = -1.0£+10 $3 RETURN 
2 XN=TCRT-TTRP $ YN=DTRP=DCRT $ X=(TCRTI-T)/XN §$ OXOT=-1.0/XN 
eee STE) OS XEN = EL*XEZX $ V= X/XE .$ Vax = (1-EL) /XE 
9 Y = AWM1) + AW(2)*V + AW(3)*X $ Y4AX = AW(2)¥*VLX + AW(3) 
10 W= XE =- X 8 W1X = XEL - 4 $ ODENLIQF = DCRT + YN¥(X + W*¥Y) 
15 DOSDOT = YN¥*DOXDT*#(1 + W*Y4X + W1X*Y) $ RETURN $ END 


FUNCTION ESATF(T) 

FORMULATED TO SAVE COMPUTER TIME. VALIO FOR TTRP £ T § TA. 
SSayeline = si= TAX 4+ BeX2 + C#XS 4°. . 
COMMON/10/ NFE, ETRP, EK(9) 

DIMENSION A(7) 

DATA (TTRP=90 6348) 5 (XN=214.982) 5 (ETR=53170541) 

DATA(A = 26749373135 0620467040, -0.64844727, 1.98167119, 
1 -1.77279598, Ue448017265 0431994366) 

WSR PIAXN de 4 9S OO) 2 J=4)5 7, 
E + A(J)*X¥**J § ESATF = E*ETR $ RETURN $ END 


4 X 
2 & 


SUBROUTINE =SATFIT 
GET ESAT DATA BELOW TA VIA ESATSUM(T) 5 FITe X = (T=TTI/{TC-TT). 
BSATAZETRP = d= AMX Be Xi2 4 CXS + 3  « 
COMMON/1/GK»AL »BB,BEyGAgDEgEP, DCRT,TCRI yPCRT »DGAT,O0TRP,TTRP»PTRP 
COMMON/7/ TAsTBePBsRBe EAsHAsSAy EByHBSB 
COMMON/10/ NFE»s ETRP», EK(9) 
COMMON/999/7 NFUNY oF (30) 
DIMENSION TT (35) 5XX%35) »ESX(35) 
1 FORMAT(1H4 16X *FIT ETHANE ESAT(T), ETRP =*F1003//3(10X3E172.9/)/ 
1 17X SHT>K 7X3HESX SXSHCALCD 6X4HPCNT ) 
2 FORMAT{10X Fiu.35 2Fide1, Fi0.3) 
3 FORMAT(L1HG 12X GHNP =I3, 10H, RMSPCT =F7e4) 
7 ETRP = ESATSUM(TTRP) $ MN = TCRT-TTRP $ ODO 8 K=159 
8 EK(K) = u $8 NP = 34 $$ NFUN = NFE = NF = 7 
10 DO 20 J=1,NP $8 IF(JeEQeNP) 13,14 
Psa le=— TAM’ S) GO M0 15 
144 T = 96 + 5*J 
45 T1(J) = T BS X = (T-TTRPISXN $ XX) = X 
46 FC(1)=X $ FI(2Y=EX*K BF FC SD=EX*X*X $F FCGI=F(S)*X SF FISD) =F (4) *X 
17 Flo) = F(5)*X 8B FI(7) = F(6)*X 
18 ESX(J) = ESATSUM(T) $ Y = ESX(J)/ETRP = 1 
20 CALL FIT $ CALL COEFF $ DO 24 K=1,NF 
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21 EK(K) = FC(K) $ PRINT 1» ETRPs (EK(K) »K=1,9) 
24 SS = 0 $$ OO 30 J=1»yNP 

25 X = XX€J) $8 EC = 1 $ ODO 26 K=1,NF 

26 EC = EC + EK(K)*X**K $$ EC = ETRP*EC 

27 PCT = 100.0*(1-ECVESX(J)) $$ SS = SS + PCT*PCT 
30 PRINT 25 TT(J)_5 ESX(J)9 EC, PCT 

$1 SS = SQRT(SS/NP) $ PRINT 3,5 NP, SS 

99 RETURN $ END 


FUNCTION ESATSUM(T) 
START WITH E = EA AT TA = 256.562 Ky THEN - 
USE CSATXF(T), PSATF(T), DENLIQF(T), FOR SUMMATION - 
DEL(E) = SUM(CSAT*#DT) + SUM(P*DD/DS**2). 

COMMON/5/ ODSDT 

COMMON/7/ TAsTBsP3sRBs EAsHAySA, EBsHB9SB 

E= EA $ TR = T-TA $ N = ABS(TRI72 + 2 $ ODOT = TRIN 
00 6 J=isN $$ TJ = TA + (J2025) *DT 

cs CSATXF(TJ) $ PS = PSATF (TJ) 

OS DENLTIQF (TJ) $ DD = DOSOT*OT 

—E = —€ + CS*0T + 100*PS*00/0S/0S 

CONTINUE $ ESATSUM = E $ RETURN $ END 


aang 


OWE W 


FUNCTION ESUMF(N»T,0A,0B) 
C GET DELTA E OVER DENSITY RANGE FROM DA TO D3. 
Cc DELTA E = INTEGRAL(P=T*{0P/0DT)) *OX/0N**2. 
COMMON/3/DPUT ,D2PGT2,U0PSOT,OPMDT,0P0D,0PDR,DTSOR,DOTHOR 
2 DX = (DB=-DA)/N $ ESUMF = 0.0 
3.00 5 J=iysN $8 ON = DA + (J20.5)*DX $ P = PVTF(T,DN,0) 
5 ESUMF = ESUMF + (P=T*OPOT) *O0X/0N/0N 
9 ESUMF = 100*ESUMF $ RETURN $ END 


FUNCTION FINDENF{T,P) 

C ON ISOTHERM T, FIND DEN, MOL/L, TO MINIMIZE (P=PC) VIA EQNSTATE. 
COMMON/1/GK,AL5BB,BEsGA,DE,EP, DCRTsTCRT »PCRT»DGAT ,DTRP»TTRPgPTRP 
COMMON/3/0PD0T,D2P0T2,0PSDT,DPMDT ,0PDD,0PDR,DTSOR,DTHOR 
DATA (DM = 23.0) 

44 FORMAT{1HO 9X *FINDENF 
42 FORMAT{1HO 9X *FINDENF 
43 FORMAT(1HO 3X *FINDENF 

IF(T-TCRT) 25558 

OG=DENGASF(T) $ OL=DENLIQF(T) $ PS=PSATF(T) $$ IF(P=PS) 353254 

D = DG/2 $ GO TO i1 

D0 = (DL+DTRPI/2 $ GO TO 11 

DG=DL=OCRT $ PS=PCRT $ IF(P=PS) 653397 

D0 = DCRT/2 $3 GO TO iif 

D= 2*O0CRT $ GO TO 11 

IF (T LT.400.0) 9519 


Q0, FAILS TO CONVERGE. * / ) 
OCRT, DP/DR ZERO OR NEGe * / ) 
0, DOUBLE-VALUED AT P = PSATe * / ) 


WEAWNP 
wou u 


on OC} 
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9 PC = PYTF(T,DCRT,0) $ IFC(P-PC) 0533,7 

10 D = DCRT 

11 DO 39 J=1,554 $ JP=P-PVTF{TsD051) $ IFCABS (DP/P)-1-0E-6) 31531512 
12 IFC(OPOD) 34534513 

13 DD = OP/OPOG $ IF(ABS (00/D)-1-0E-6) 31531514 
144 D= D0 +# DU $ IF(D-GT.000) 16515 

15 0 = OGAT 5 60 TO 30 

16 IF(DeGTe-DM) 17518 

17 0 = DTRP $ GO TO 30 

18 IF(T-TCORT) 19,24,30 

19 IF(PeLTePS) 20522 y 
20 IF(DeGT.DG) 21430 

21 0= DG $ GO TD 3D 

22 IF(OLTeOL) 23,30 

23 D= OL $ SO TO 30 

24 IF(PelLTePCRT) 25927 

25 IF{D.LT.DCRT) 30526 

26 D = OCRT - 0-62 $ GO TO 30 

27 IF(DeGTe-DCRT)I 30528 

28 D = DCRT + 0.02 

33 CONTINUE $ FINDENF = 0 $ PRINT 41 $ RETURN 


31 FINDENF = D $ RETURN 

$2 FINDENF = 9 $ PRINT 43 $ RETURN 

33 FINDENF = OCRT $ RETURN 

34 FINDENF = DCRT S$ PRINT 42 $ RETURN $ END 


FUNCTION FINDTMF (P) 
C GIVEN MELTING PRESSURE» Py GET TEMPey To 
COMMON/1/GKy AL 9BByBE»sGAyDEyEP,y DCRTyTCRT»PCRT»DGAT,OTRP,TTRPgPTRP 
DATA (A=2559.€5) »(E=2.179) 
1 X = (P-FTRP)/A + 4 = $ FINDTMF = TTRP¥X**(4.0/E) 
2 RETURN $ END 


FUNCTION FINUTSF(P) 

C GIVEN VAPOR PRESSURE Ps» ITERATE T TO MINIMIZE (P-PC). 
COMMON/1/GKyAL»BBsBEsGA,DE,EP, DCRTsTCRT »PCRT»DGAT,OTRP,TTRP»PTRP 
COMMON/ 3/0P6T, 02P0T2»,0PSOT,DPMOT ,0P0D,0POR,DTSOR,OTHOR 


1 FORMAT(41HO 9X *FINOTSF = 05 FAILS TO CONVERGE. * / ) 

2 FORMAT(1HO 3X *FINOTSF = 0, P EXCEEDS PCRT. * / ) 

3 IF(P-PCRT) 4,11 512 

4 T= 200 $8 BOY J=1,50 $$ DP = P = PSATF(T) $ AOP = ABS (DP) 
5 IFCADP/P-1.-J3E-6) 1051056 

6 IF(ADP/DPSOT/T=-1-0E-6) 1051057 

7 T= T + DP/JOPSDT $ IF(T-TCRT) 959,8 

8 T = TCRT 

9 CONTINUE $ FINDTSF = 0 $ PRINT 1 $ RETURN 

46 FINDTSF T $ RETURN 


11 FINDTSF 
12 FINOTSF 


TORT $ RETURN 
G 3 PRINT 2 $ RETURN $ END 
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SUBROUTINE FURTADO 

C COMPARE FURTADO CP(P,T) ON HIS 14 ISOBARS. 

C READ PRESSURES FIRST. SEPARATE ISOBARS(CARDS) BY BLANKS. 
COMMON/1/GK5AL5BBsBE9GAaDEsEP,y DCRTsTCRT sPCRT »DGAT sDTRP»TTRPsPTRP 
COMMON/7/ TAsTBsPBsRBo EAsHAeSAy EBs,HB»SB 
COMMON/8/ INogIKe PoTsDENsy EgHoSy CVsCP,CSATs WyWK 
COMMON/9S9/ TI5EZZ,y EZsSZ9CVZsy HZ,CPZ 
DIMENSION PF (14) 

11 FORMAT(1H1i) 

12 FORMAT(1HO 16X *FURTADO ETHANE CP(P,T) AT P =* F8e35 4H BAR & 
1 17X3HT,K SXSHMOL/L 3X7HI/SMOL/K SXSHCALCC oX4&HPCNT) 

13 FORMAT{I5, 3F10.0) 

14% FORMAT(10X 2F10.3,5 3F1I0.2) 

15 FORMAT(8F1i0.9) 

21 READ 15, (PFC(I),I=1,14) $B PRINT 11 See e= 90 

22 DO 50 T=1514 $ P = PF(I) $ PRINT 125 P $ L = Lt3 

23 00 4&0 J=1,95 $B READ 135 IDDsTyPX,CPX $ IFC(IDD) 24,50 

24 CPX = 1.874335*CPX $$ L =Lti $G$ IF(LGT.55) 25526 

Coa $ PRINT 12 

20 IF(IeEQe1) 27,31 

2 TiS 1 $ CALL IDEAL $ CP=CPZ $ DEN=0 $ GO TO 39 

S31 IF(TLT.TB) 32,38 

32 IF(TeLTeTA) 33534 

33 CALL LIQTHR’ $ GO TO 39 

34 IF(TeLTeTCRT) 35,537 

3p PX = PSATF(T) B IFCPeGTePX) 36,38 

36 CALL THRMOX $ GOTO 39 

37 PX = PVTF(T,DCRT,6) $ IF (PeGTePX) 36,538 

38 CALL THRMOGS 

39 PCT = 106.u*(1-CP/CPX) 

4&3 PRINT 145 TsDEN, CPX ,CP, PCT 

50 CONTINUE 5 RETURN $ END 


SUBROUTINE HOMOTRM ae eae 
VALIO ONLY AT T = INTEGER MUL ° 
GIVEN Pet FOR HOMOGENEOUS DOMAIN, GET DEN ANDO FUNCTIONS. 
SAVE DK(T) 5 ZKCT)» SK(T) 5 CKIT) ALONG EACH ISOBARe 
COMMON/3/0POT»D2POT2 »DPSOT,DPMDT »DPDDsDPDR»DTSORsDTHOR 
COMMON/8/ INgIKo PoT DEN» EsHsSo CVyCPyCSAT>s WsWK 
COMMON/99/ TIyEZZs EZsSZCVZs HZ,CPZ 
DIMENSION DK (7G) 9EKC7 G0) 9SK(70) »CK(70) 
DATA (Q=1.01325) 5 (G=0-0831434) 
1 IF(IKeEQeIN) 2910 & 
2 TI =SU LS. MCALE LOEALR Oe Sue Ib— T/id $ ODA = 0 
3 DK(J) DEN = DB = FINDENF(T,P) $3 N = pa*20 + 10 
4 
5 


a00 


= ZZ + EZ + ESUMF(NyTs9A,0B) 
en) : = SZ + SSUMF(1s9N,T,DA,08) - 400*G*ALOG(G*T*0B/Q) 
6 CK(J) CV = CVZ + CSUMF (N»T»DAs03) $ GO TO 20 
40 J = 1/10 $ DA = OK(J) 
41 OK(J) = DEN = 038 = FINOENF(T,P) 
12 DX = ABS(D0B-DA) °) N = (OX41)*5 


13 EK(J) = E = EK(J) + ESUMF (N9T »0A,DB) 

14 SK(J) = S = SK(J) + SSUMF (O09N9T»0As 08) 

15 CK(J) = CV = CK(J) + CSUMF (N»T»sDA,08) 

ara 12 GoseeRywss) i leis = PVTF(T,DBs1) 

Zt (OP = OW wv 150*T/0P00* (OPDT/0B) *¥*2 

22 W = SGRTCWK*CP*DPOD/CV) $ RETURN $ END 
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SUBROUTINE IDcCAL 


C ETHANE ITOEAL GAS (1 ATM) THERMOFUNCTIONS (CHAO, 1973). 
C EQN LY) = G22 20) 7Ri = So taeralX) gumXe =a A100) Qe Xee 17 Ss 
C F(X) = A1L¥*Q4 + AZ*Q5 + AZ*¥Q6 +t © « « + AN*Q**¥(Nt3) © 
C LZR se VG NING SW GVIZI RY =) VeZ/7R) 700 5eGeZ7Re= 1 * CVZ7Rs 
C SZ/R = AZ + INTEGRAL(CVZ/R/X + 1/X)*DXe 
COMMON/99/T1I4 EZZy EZySZyCVZy HZ,C0PZ 
DIMENSION A(9) 
DATA (R=8. 5143) 5 (AZ=21.705718) 
DATA(A = 652498641, -362.0115914, 853.3408616, -1123.601€22, 
1 9060218447 $459023502545, 143.0300226, -25.07495605 516897540044) 
De nd or Ole CUBERM  GXo) 6d.) OQDX. = (O43 4X VSiike= Fd) = S = 6 
20900 4 J=1,9 $8 K = J#3 $ L = Jt6 
SV = ACIDVAMvMK § FR SF oY 83  iesh (esl EAN @) 
hy S Sosy ca (Sa7ax $° S = AZ + S + 4FALOG(X) 
5) \f iy eS $ CV = Y + X*F1*DQDX 
C CONVERT TO DIMENSIONED RESULTS, JOULES, MOLES, KELVINS. 
Cece Ren LY  S OHZ = RET EAC) CPZ, =. R* (1+CV) 
9 SZ = R*¥S $ CVZ = R*CV $ RETURN $ END 


SUBROUTINE ISOTHRM 
C PRINTOUT THE CRITICAL ISOTHERM. 
COMMON/1/GK»AL58B,BE2GA,DE,EP»y DCRTsTCRTsPCRT,0GAT,0TRP,TTRP»PTRP 
COMMON 3/0P9T , D2PDT2» DPSDT,0PMDT »9PDD0,DPOR,D0TSOR,DTHDOR 
COMMONS4S X3159XB2_y XC1,XC2,y XE1,XE2_5 OX30R,DXCOR,OXEDR 
COMMON/6/ TSAT»s THETA, PSAT 
1 FORMATC1H1 14X *THE CRITICAL ISOTHERM* // 
1 6X4&HTC =Fle2, BH» BC =FO el, 6H, PC =File/7y 
2 *. AT THE CoPey SPS/DT =* Fed, GHy DP/DT =F8e5// 
3 6X4HDJDC 9XSHTSJTC OXSHPS/PC A0X4HP/PC 9XS5HDP/OR &X6EHDTS/DR 
G4 YXEHOTHJOR 4SXOHDPS/JDOR &X6HDXB/DR 4SXOHDXC/OR ) 
FORMAT@I2X FBe2,y SFL4010, Fl4e9, 5F1065) 
PC = PVTF(TCRT,OCRT,0) $ DXCDR = G 
PRINT 1,5 TCRI,DCRIT,PCRIT, OPSOT,DPOT $ DO 8 J=1,51 
= 0074 + 0.01*J & ON = OR*DCRT 
PR = PYTF{(TCRT»0N51)/PCRT $$ OPSOR = DPSOT*OTSOR 
TSN = TSAT/TCRT $ PSN = PSAT/PCRT 
PRINT 25 DORs TSNoPSN» PRsDPDR, DTSDRsDTHOR,DPSOR, 
1 DXBDR» DXEDR 
9 RETURN $ END 


ON OVE Wh 
Oo 
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SUBROUTINE STLOCUS 
DERIVE THE J-T INVERSTON CURVE. USE ROUTINE DELTAF(T,OI). 
DIMENSION TT (99) yPP(99) ,0N(99) 
VATA (BUCRT=6.80) ,(TCRT=305.33) 
1 FORMAT(1H1 1eX *THE JOULE-THOMSON INVERSION LOCUS FOR ETHANE*// 
Ll 17X3HT,K SXSHP,»BAR SXSHMOL/L 7 XSHTyK SXSHP,»BAR SXSHMOL/L) , 
2 FORMAT(10X Tits Fidelsy FLVe2y, Tidy FiGei,y F102) 
6 TA = 240 $$ NP = 72 
( ERIUNG Dds DO, eer haa INR id OX he 
Ss ye VN eo ber 8 we SB Y/ICRU 
J DI = DORT*# (e645 - VUeDS*K + JVe2h/X) 


du IF (T=$TCRT) 11,12,12 

i1 DL = DEWLIQF(T) $ IF (DI-DL) 25512,12 

12 SS = DOELTAF(T,DI) $ DO 20 IT=1,15 

44 D=OI-9X $ SL=E=DELTAF(T,D) $ O=OI+DX $ SP=DELTAF(T,D) 


59 IF(SS*SL) 13,916,140 


to IF(SP-SL) 13517,17 

Lf SS 2 Sk 8 OL oe DH > DM 5 CO TO 2e 

13 IF(SS-SF) 23520,19 

1 SS = SP 5 Dl & Ju > DK 

2s OX = OX/2 

23 Ti(I) = T 5 ON(I) = OL $ PP(I) = PVTF(T,DI,0) 
2> CONTINUE $8 N = NP/2 $ ODO 29 J=15N 

eco IT = T(J) B&B ITT = TT(JUtN) 

23 PRINT 2, ITsPP(J),DN(U)»9 ITT ,PPC(JU+N) ,ON(JUtN) 

36 RETURN $ cNO 


AN Om WHR 


SUBROUTINE LIGTHRM 


FOR DEN > RB=14e5 MOL/L, AND T § TA=2560562 Ke 
GIVEN Pet, GET DEN AND ALL FUNCTIONS. 
START ON SATIeLIGe, THEN INTEGRATE ALONG ISOTHERM Te 


COMMON/3/0PDT ,02PLT2,0PSDT,DPMDT ,0P00,0PCR,0TSOR,OTHOR 
COMMONZ5/ DOSDT 
COMMON/8/ INeoIKso PeTsDENsy EoHoSe CVsCP»CSATs WeWK 


—E = ESATF(T) $, S = SSATFIT): $ DA = JENC TORI) 

DBDLOT = ODSST $ GEN = OB = FINDENF(T,P) 

PX = PVTF(T,0A,0) §$ CV = CSATXF(T) + 100*T*OPOT*DDLOT/DA/DA 
OX = ABS(DB-DA) 8 N = (9X+1)¥*5 $$ -E = E + ESUMF(N,T:DA,0B) 
H = — + 100*P/0B $$ S = S + SSUMF(0yN,T,DA,DB) 

Cv = CV + CSUMF(N,T,0A,0B) $$ PX = PVTF(T,0By91) 

CF = Cv + 100*T/DPDD* (DPDT/0B)**2 

W = SQRT(WK*CPFUDPOD/CTV) $ RETURN $ END 


SUBROUTINE PEEK 

EXAMINE BEHAVIOR OF THE COEFFICIENTS. 

B(S) = B1 + B2*S + B3*S2, E@S) = SN¥SR*1EL + E2*S + ES*¥S2). 
COMMON Bis82.B35 CisC 29035 EL,sE2,E3_5 ER 
COMMON/1I/GK,AL 9BBy BE sGAQDE,sEPy DCRTyTCRTyPCRT,DGAT,D0TRP»TTRP,PTRP 
CCMMON/6/ TSATs THETA, PSAT 


G FORMAT{(1H1 14X *EQUATION OF STATE, COEFFICIENTS* // 


1 15K GHOTRP =F8e%, BH, TITRP =F8e35 BH. PTRP =E1609/ 
2 15X OHODCRT =F 8e4_ 8H, TCRT =F8e3y 8H, PCRT =F13-9// 
3 15X 4HAL =F5.25 6H, BB =F5e2s5 6H BE =FSe2/ 
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& 15X GHGA =F5e2s 5He DE =F5e2y OHe EP =F5e2// 
5 2(13X 3Fi6.11/) ) 
> FORMAT(15X SHMOL/L 6X4HTSAT SXSHTHETA 6X4HPSAT 9X1HB 9X1HC ) 
Oo FORMAT(10X FiGelsy 5710.3) 
72 PRINT %s OTRPsTTRPsPTRP»s DCRT»sTCRT»yPCRTs ALyBB,BEy GAyDE EP, 
1 B1,3829B35 ElsE25c3 
GERUNDS ed) SOOTK, I-41 5445 055 ONG= Deb) 9S) R = Se= DNZOGRI 
We IN BORIS SS SN eS Sol SR eS Sosy ss fe Sf) 
(Seba Bit Belk + BSR $ ERS SNMORE CEAl + TEER + TE Se RZ) 
76 TS=TSAT=TSATFCON) $ TH=THETAF (ON) $ PS=PSATF(TS) 
77 PRINT 6% ONs TSeTHsPSy ByE $ RETURN $ END 


FUNCTION PMELTF(T) 
ETHANE MELTING LINE» BARe STRATY/TSUMURA, 19756 
SIMON EQNeo P = PTRP + A*(X**E21), X = T/TTRPo 
COMMON/1/GKsAL»BBsBEsGA,0E,EP, DCRT»,TCRT »PCRT »DGAT yDTRPsTTRP»PTRP 
COMMON/3/DPUT,02PDT2,DPSOT,DPMOT »ODPDD,DPOR,OTSORsOTHDOR 
DATA (A=2555265) 5 (E=20179) 

ho B UAUIRE Se rts $ PMELTF = PTRP + A*(XE=1) 

3 DPMDT = A*E*XE/X/TTRP $ RETURN $ ENO 


oOG 


FUNCTION PSATF(T) 
NEW ETHANE VAPOR PRESSURE EQUATION, FALL OF 1975.6 P IN BARe 
EQNes LN(P/PT) = PL*¥X + PZ2FU + P3*U2 + P L*US + POFU*(1-U)**ES 
ARGSe, U = (T-TTI/(TC-TT), X = (L°TT/TI/¢01-TT/TC). 
COMMON/3/0P5T,02P0T2,0PSDT,0P MOT »,OPD0;0PDR,0TSOR,DTHOR 
DATA (TTRP=390.348) » (TCRT=305. 33) 
DATA (EPP=1.3) » (PTRP=1.130842618E=-5) 
DATA (P1=18.84522826756) » (P2=-706354151345) 
1 (P3=50 4284431006) 5 (P4=—-1. 36235270362) » (P5=0. 7692492586) 
1 FORMAT@LHS 9X *PSATF = Oy T EXCEEOS TCRTs * / ) 
2 UN = TCRT-TTRP $ XN = 4 = TTRP/TCRT 
3 U = (TeTTRF)/UN $ X = (12TTRP/T)/XN 
4 U1T = 1.6/7UN $8 X4iT = TTRP/XN/T/T $ U2 = U*U $F$ UZ = U*V2 
6) 
U 
8 


aan 


V = 1-U $ IF(V) 75899 
PRINT 12 $ STOP 
Z2= Zi = G6 x) GO TO 106 
9W V¥FEPP $3 Wi = -EPP*W/V $ Z = U*W $§$ ZL = U*¥W1 + W 
SL) fe P1*X + P2*U + P3S*US + PHFUS + PO*Z 
41 FiT = Pi*xX1iT + (P2 # 2¥P3¥*¥U + S¥*P4Y*US + P5*Z1)*U1LT 
12 PSATF = PTRPY¥EXP(F) $ DPSDT = FiT*PSATF $ RETURN $ END 
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Qa9NddAM 


(Sp) (2) (39) Ge 


FUNCTION FVTF(Ts0.™) 

PVTF = P,yBAR. M=0 YIELDS OP/DT,02P/0T2. M=1i YIELDS ALSO OP/OD. 

NOTE, R = S = DEN/DOCRT, ANDO GK = 020831434*DCRT »BAR/Ke 

EGNSTATE, P-PSAT = S¥GK*(T=TSAT) + S*¥S*GKFTCRI*F(R,T)5 WHERE = 

F(Rs7T) = BCS) *#XBF(R,T) + E(S)*XEFCR 5T)5 AND =- 

BiS) = Bi + 32¥*R + BFR, E(S) = SN*¥SR*¥(E1 + E2*R + ES*R2)0 
COMMON 81582583, CisC2sC35 E1,E25E3,5 ER 
COMMON/1/G6K,AL9BBsBE»GA,DEsEP, DCRT,TCRT »PCRT sDGAT,0TRPsTTRPsPTRP 
CGMMON/3/DFUT ,O02P0T2Z,0PSDT,0PMDT ,0PDD,0PDR,0TSOR,DOTHOR 
COMMON/S4GS X3B19XB25 XC1iyXC2_o XEL1,XE2_ OXBER,DXCOR,DXEDR 
COMMON/6/ TSAT» THETA, PSAT 
R=S=D/BCRT $ R2=R¥*¥R $F$ RS=ER*FR2 $ SN = Sei $ SR = SER 
TC = TCRT $& RG = R¥GK $ GKT = GK*TC $ ODSDR = CTRP/OCRT 
TS=TSAT=TSATF (D) $ FS=PSAT=PSATF(TS) $ THETA=THETAF(D) 

XB = XBF(T,)0) 5 XE = XEF(T,D) 
B = BLFR2 + B2*RS + BS*FRFERZ $ SM = SN*¥R2*SR 
eS) 6 [asl <P EBUiR © [Ese $ —E = ESM*SM 

10 F = B¥XB + ZEFXE $ Fi = B*X31i + E*XEL 

sol ipa & BKB © EVES $ PVTF = PS + RG*(T2TS) + GKT*F 

44 DOFDT = RG + GK*F1i $ O2PDT2 = GK*F2/TC $ IF{M) 15330 

15 BD = (2*Bi*¥R + 3*B2*R2 + 4*B3*R3Z)*OSOR 

ale) (2Sab S Be & BZviESVR $ SM1 = SN*R2 + (3*R2 |= 2*RI*FSR 

147 —€D = (ESM*¥SMi + ES1*SM)*OSOR 

23 Fi = B¥DOXBOR ¢ BD¥XB + E*OXEDR + ED*XE 

26 DFDOR = (DPSOT=RG)*OTSOR + (T=-TS) *GK*DSDOR + GKT*F1 

ef OPDG = DPUR/JDTRP 

30 RETURN $ END 


Ome wnror 


FUNCTION QVAPXF(T) 
VIA PROGRAM QVAP-3, 2/9/76, USING DATA OF = 
DOUSLIN» RIEDEL» FURTADOe FOR 49 POINTS, RMS = 0-56 PCT. 
DEFINE X = (TC-T)/(TC=TT) ,y U = X*¥*4/3, WHEN EQNe IS = 
QV = A1¥U + A2*U2 + « « « + AG*U6. 
DIMENSION Alo) 
DATA (TT = 96.348),(TC = 305.33) 
DATACA = 12,.857784%, 524095302, 3346835851, 
1 -972520486, 100684976, -37-933878) 
1 U = CUBERT ((TC-T)/(TCHTT)) $ Q= 0 $ DO 2 K=196 
2 Q= Q & A(K) *U*FFK $ QVAPXF = 1000*Q 3 RETURN $ END 
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SUBROUTINE SATGSTR 

C GIVEN PsT AT SATURATION, GET DEN ETCe - 
COMMON/1/GKsAL»BB,BE,GA,DE,EPy DCRT,TCRT »PCRT »DGAT,O0TRPsTTRP»PTRP 
COMMON/Z3/DPOT, D2PDT2,5 DPSOT,D0PMDT » 9POD,DPOR,DTSOR,OTHOR 
COMMON/8/ INglKe PoTsDENy EgHoSs CVyCP5CSATs WoWK 
COMMON/99/ TI,EZZ_, E2sSZsCVZ5 HZ5CPZ 
DATA (DA=0 64) 5 (Q=1201325) 5 (G=0008514354) 


20 FORMAT(1HS 9X *T cXCEEDS TCRT IN SATGSTRe* / ) 

1 TI = 7 $$ CALL IDEAL $ ODEN = OB = DENGASF(T) 

2 N = O8*¥20 + 10 

See Cte tee SUMENC Nels DAGUE)) eS muhe= G+) LUO PDB 
& S = SZ + SSUMF(19N9T9DA,0B) = 100*G*ALOG(G*T*DB/Q) 

5 IF(T-TCRT) 996912 

6 CP = CV = W= 0 $ RETURN 

9 CV = CVZ + SSUMFIN,yT,0A,08) $$ PX = PVTF(T,0B,1) 

iu CP = CV + 1J0*T/DPDD*(DPDT/03)**2 


141i W = SQRTC(WK*CP*OPDD/CV) $ RETURN 
12 TPO ai) 5 SUG 85 = TaN) 


SUBROUTINE SATLAQTR 
C GIVEN PeT FOR SATe»LIQe, GET DEN AND ALL FUNCTIONS AT T € TA. 
COMMON/3/0P3T ,D2P0T2,0PSDT ,DPMDT ,OP00,0PER,0TSOR,OTHOR 
COMMON/5/ DOSDT 
COMMON/6/ INoIKs PoTsCENs EsHoSy CVeCPyCSATs WsWK 
2 E = ESATF(T) $ S = SSATF(T) $& CSAT = CSATXF(T) 
2 DEN = DENLTQF(T) $ ODLDT = DOSOT $ H = E + 100*P/ DEN 
3 PX = PVTF(T,DENs1) $8 CV = CSAT + 100*T*OPDT*DDLOT/DEN/DEN 
& CP = CV + 100*T/DPDD* (DPDT/DEN) **¥2 
& W = SQRT{(WKFCP*DPOD/CV) $ RETURN $ END 


FUNCTION SSATF(T) 
C START WITH S = SA AT TA = 256.562 Ky THEN - 


C USE INTEGRAL FORM OF RODER CSATXF(T). GET CHANGE FROM TA DOWN TO T. 
COMMON/7/ TAsTBsePBsRB,y EAsHAsSAy EByHB9SB 


CCMMON/9/ AU, BUsCUsDUy ASyBSyCSyDSyESy ETBsST3,STA 
1 F = (305633-T)**0.4 


3 STT = AS*ALOG(T) + BS¥*T = CS*F/00% = DS/T = ES/2/T/T 
4 SSATF = SA + STT = STA $ RETURN 3 END 
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C 


C 


ao 


OWFWNE 


i] 


4% 


5 
11 


12 


140 
141 
142 
143 
144% 
145 
140 
147 
143 
149 
150 
pus al 


FUNCTION SSUMF (L—oNoTo0A,0B) 
DENSITY=DrPENCENT CHANGE OF S FROM ODA TO OBe 
DELTA S = INTEGRAL(GK=-(DP/DT)/0N)*OX/DNe 

COMMON/3/0F0T ,02PDT2,0PSOT,DPMDT ,vPODD,D0PCR»DTSOR,DTHOR 
DATA (GC = J2v831434) 

SSUMF = @ $ OX = (DB-DA)/N $ IF(LeEQe0) 4o2 

DO 3 J=lyN $8 ODN = DA + (J-0-5) *DX $ P = PVTIF(T,0N,0) 
SSUMF = SSUMF + (GC-ODPOT/ON)*OX/ON $ GO TO 9 


DG 5 J=15N 8 ON = DA + (J2005)*0X $ P = PVTF(T,0N,0) 
SSUMF = SSUMF = DPOT*OX/ON/ON | 
SSUMF = 100*SSUMF $ RETURN $ END 


SUBROUTINE TABLIQ 

TABULATE THE ETHANE SATURATED LIQUIO FUNCTIONS. 
COMMON/1/GK,AL5BByBE 5GAgDEgEPs DCRT,TCRT yPCRT ,DGAT »DTRPoTTRP»PTRP 
COMMON/3/DPDI ,02P0T2 5 D0PSDT,DPMDT »,DPDD,DPDR»,OTSOR,DTHDR 
COMMON/5/ OUDSDT 

COMMON/7/ TAsTBsPB»RB»y EAgHAySAy EB,HB,SB 

COMMON/8/ INeIKy PoTsDENs EsHoSo CVyCPyCSATs WyWK 

DIMENSTON TSA(46), PSA(46) 

FORMAT¢1iH1 13X *PROPERTIES OF SATURATED LIQUIO ETHANE* // 

1 24K2HT 10X1HP SX3HDEN 4&X5SHV,LIQ OXSHV5GAS SXOHDPS/DT SX6HDDL/DT 
2 &XSHDP/OT cXSHDP/DD 2X5HQ,VAP 2X5HQ,XPT / 

3 10X5HODEG K BX3HBAR SXSHMOL/L 4XSHL/MOL OXSHL/MOL 6XSHBAR/K 
4 2@X7THMOL/L/K &XSHBAR/K 2X9HBAR|=L/MOL 2XSHJ/MOL 2XSHI/MOL ) 
FORMATOSXF1ise 3, ELleds FBe3S_y F9e5sy Ell eds EL103y 2F904,E11-3,2I7) 
FORMAT(1iH1 13xX *PROPERTIES OF SATURATED LIQUID ETHANE* // 

1 14X1HT 11XiHP 9X1HE OXLHH 9X1HS 
2 ©X2HCV BX2HCS OX2HOCP OXLIHW 4XOHCS,XPT / 

3 18X5SHDEG K 9X3HBAR SXSHJ/MOL SXSHJI/SMOL 3X7HJ/MOL/K 
G AX7HIJSMOL/K ATX7THISMOL/K AX7HI/SMOL/K 2XDHM/SEC SX7HJ/MOL/K ) 
FORMAT@5X Fi0deSs E12 04s 2FLID0e1, Fi0035 S3F8e2s I75 Fide?) 
FOR PAGE ONE OF TABLIQ. 

PRINT & 3 NP = 4&5 

DO 151 J=1i,NP $ IF (JeEQe1) 1425143 

T = TITRP $ GO TO 147 

IF(JeEQeNP) 144,146 

T= TCRT $ OG = BDL = OCRT $ ODOLOT = 0 

VG = VL = 1/0G $ GO TO 149 

Y a BH oc Su) 
= CENLIQF(T) $ DOLOT = DOSDOT 
0G = DENGASF(T) F) VG = 1/0G $ VL = i/7DL 

TSA(J) = T $ 2SA(J) = PS = PSATF(T) 

TQC = 100*T*DOPSOT* (VG-VL) $ IQX = QVAPXF({T) $ PX=PVTF(T,DL 51) 
PRINT 55 ToPSsOL»y VL»VGsy OPSOT,»,O00LOT, DPOT,DPDD, IQC, IQX 
FOR PAGE TWO OF TABLIQ. 
USE SATLQTR AT T < TA, USE THRMXL AT T > TAe 

PRINT 11 $ ODO 165 J=1,NP 

P=PSA(J) $ T=TSACJ) $ CSX = CSATXF(T) $B IF(TeLToTA) 162,163 
CALL SATLQTR $ IW =W § GO TO 165 

CALL THRMXL $ IW= W 

PRINT 125 ToPy EosHoSy CVsCSAT,CP, IW, CSX 

RETURN $ END 
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SUBROUTINE TERRES 

C COMPARE SOUND VELOCITIES IN HOMOGENEOUS OCGMAIN. 

C ID = (3) TERRES ET ALey BRENNSTOFF CHEMIE 38, 129 (1957). 
COMMON/1/GK,AL>8BsBEsGA,DE,EP>s DCRT,TCRT sPCRT »,DGAT »DTRP,TTRP,»PTRP 
COMMON/7/ TAsTBePByRBy EAyHAySAy EB yHB9SB 
COMMON/5&/ INeoIKs PoTeDENs EsHoSy CV»oCPeoCSAT, W WK 

1 FORMAT(I5, 3F10-0) 
2 FORMAT{1Hi 7X FTERRES ETHANE SPEED OF SOUND, METERS/SECe * // 
1 8X 2HID 7XSHT 9K 5X5HPyBAR SXSHMOL/L SXSHW,XPT SXSHCALCD 4&X4HPCNT) 
3 FORMAT(SX I55 3F1003, 2F1De1», F8e2) 
6 PRINT 2 $ N =O $ ODO 20 J=1,300 
7 READ 15 IDs TCsPKsy WX §$ IFCID) 38,530 
8 N= Nei $ IF (NeGT046) 9,10 
9 N 6 $ PRINT 2 

10 T 273.15 + TC $ P = 0-69806605*PK 

11 IF(TLT.T3) 125918 

12 IF(T.Lt.TA) 13514 

43 CALL LIQTHRM $ GO TO iQ 

14 IFC(TeLTe-TCRI) 15,517 

15 PX = PSATF(T) $ IF (PeGT.PX) 16,18 

146 CALL THRMOX $ GO TO 19 

17 PX = PVTF(T,OCRT 0) $ IF(PeGTePX) 16518 

18 CALL THRMOGS 

19 PCT = 100.0*(1-W/WX) 

29 PRINT 35 IDy ToP»oDENe WXo9Wy PCT 

33 RETURN 5 END 


FUNCTION THETAF(DEN) 

THETA = TSATFEXP(U(S)). 

LET Q = (S-4)/(ST-1)5 WHERE ST = OTRP/OCRT, THEN - 

IF S < Ly U = AL¥FQ¥**3, IF S$ > 45 U = -AL¥Q**3, 

YIELDS ALSO THE FIRST DERIVATIVE RSPe TO RHO = DEN/DTRP. 
COMMON/1/GK,ALsBBsBEsGA,DE,EPs DCRT»,TCRT yPCRT »DGAT,CTRP»TTRP»PTRP 
COMMON/3/DPDT,02P0T2,0PS0T,DPMDT ,0P00,0PD0R,0TSDR,DTHDR 
COMMON/E/ TSAT»> THETA, PSAT 


aoogn 


1 S = DEN/DCRT $ ODSDR = DTRP/DCRT $ C = DSDR=1 
Adee CS=1)7Ce Seraer="QeQy $e 1UP= AL *a* Oe 

3 U1 = S¥FAL*Q2*0SDR/C $ IF (Q) 55994 

4 U= -U $ Ui = -U1 

5 XF = EXP (VU) $ THETAF = TSAT*XP 

6 DTHDR = (TSAT*U1L + DTSOR) *XP $ RETURN 

9 THETAF = TCRT $ ODTHDR = 0 $ RETURN $ ENDO 
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SUBROUTINE THRMOGS 

C GIVEN PsT FOR HOMOGENEOUS DOMAIN, GET DEN AND FUNCTIONS. 
COMMON/3/0PDT ,D2P0T2 ,O0PSOT,O0PMOT ,0P0D,0POR,0TSOR,DTHOR 
COMMON/8/ INgIKs PoToDEN»y EgHoSy CVeCPyCSAT>, WyWK 
COMMON/99/ T1yeZZ_y EZsSZ5CVZsy HZ,CPZ 
DATA (DA=00u) » (Q=1.01325) y(G=0.0831434) 


1 T1 = T §$ CALL IDEAL $ ODEN = DOB = FINDENF{IT,P) 

2 N = OBF2U + lu 

3 — = EZZ + EZ + ESUMF(NsTs0A,0B8)>.35 H = E + 100*P/DE 
4 S = SZ + SSUMF (1oNoT90A,03) - 100*G*ALOG(G*T*0B/Q) 

7 CV = CVZ + CSUMFIN,T,0A,0B) $$ PX = PVTFIT,DBs1) 

8 CP = CV + iJ0*T/DPDO* (DPOT/0B) **2 

9 W = SQRT(WK*FCP*DPOD/CV) $ RETURN $ END 


SUBROUTINE THRMOX 

REGION X AT JEN > DORT» TA < T £¢ TBe 

GIVEN Pets GET ODEN ANU ALL FUNCTIONS. 

START WITH VALUES FROM CVSUMIT, THEN INTEGRATE ALONG ISOTHERM T. 
COMMON/3/0PUT,D2P0T2,DPSOT,DOPMDT,0P90D,0P0R,0TSOR,OTHOR 
COMMON/8/ INoIKy PoTsDENy EgHsSs CVeCPyCSAT, WyWK 
CALL CVSUMIT $ DA = DEN $ ODEN = DB = FINDENF{T,F) 
OX = A3St(V38-DA) $$ N = (OX+1)*5 

E —£ + ESUMF(N»sT,0A,0B) $ H = E + 100*P/0B 

S + SSUMF (0.0N»sT,»DA,0B) 

CV = CV + CSUMFI(NsT,DA,DB) $ PX = PVTF{T,0B51) 

CP = CV + 130*T/DPO0* (DPDT/0B)**2 

W = SQRTCWK*CP*DPOD/CV) $ RETURN $ ENO 


ang 


NOW 4 WN 
nr 
iow 


SUBROUTINE THRMXL 
C GIVEN PsT FOR SATeLIQe IN REGION X,y GET JEN AND FUNCTIONS. 
COMMON/3/0P0T ,02P0T2,0PSOT ,D0PMOT »0PDD,G0PDR,DTSUR»,DTHOR 
COMMON/5/ OOSDT 
COMMON/5/ IN IKo PoToDENy EsHoSy CVoCPsCSATs WyWK 
OATA (TCRT = 305.33) 
CALL CVSUMIT $ BA = 
DX = A3S(0d-DA) $ N 
- = — + ESUMF(N,T,OA, 0B) $ H= E + 100*P/DB 
S = S # SSUMF (05N5T,0A,08) 5 WECYokLYoVGRUY 758 
T= CV = CP = W = G $ RETURN 


DEN $ ODEN = OB = DENLIQF(T) 
= (O0X+1)*5 


CV + CSUMF(N»sT20A5D38) $$ PX = PVTF(T,0B 91) 
CV + 106*T/DPO0D* (DPDT/DB)**¥2 

CSAT = CV = 100*T*OPOT*ODSOT/DEN/UDEN 

W = SGRT (WK*CP*DPDD/CV) $ RETURN $ END 


OCOANDFE WHF 
QO 
Y 
5S 
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C 


C 


WOON OUFEF WHE 


1 
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FUNCTION TSATF (DEN) 

ITERATE T TO MINIMIZE (OEN-DCALC) VIA DENGASF (T), DENLIQF(T). 
COMMON/1/GK,AL»,8B,BE»GA,DE,EP, DCRTsTCRTsPCRT ,DGAT »DTRP»TTRP»PTRP 
COMMON/3/0PDT ,D02PDT2,0PSOT,DPMDT »,DPDD,DPDR,DTSOR,DTHDR 
COMMON/5/ DOSDT 
DATA (Q=1.5) »(FN=3.48168907) 

NOTE, FN = EXP(Q) = 160.6 
FORMAT(1HO 9X *TSATF FAILS TO CONVERGE. * / ) 

O=DEN 8 S=D/DCRT $ YN=TCRTI/TTRP-1 $ IF 1DEN-DCRT) 3530 54% 

ST=DGAT/DCRT $ F=ALOG(S) JALOG(ST)*((1-S)/(1-ST))**2 $ GO TO 5 
ST=DTRP/OCRT $ VE((S-1)/(ST91))**3 $ F=CEXP(Q*U) -1) /FN 
T = TORT/(1 + YN*¥F) 

DO 15 J=1555 $ IF(DEN-DCRT) 7530,8 
DD = D - DENGASF(T) $ GO TO 9 
DD = D - DENLIQF(T) 

IF (ABS (D0D/0) LTo1-0E-7) 16,10 
DT = DO/DOSDT $ IFC(ABS(OT/T) LT.1-0E-7) 16511 
T= T # OT $ IF(T) 12912513 
T = TTRP $ GO TO 15 

IF(ToLT.TCRT) 15514% 

T = TCRT - 6.10 

CONTINUE $ PRINT 1 $ STOP 

TSATF V $ OTSDR = DTRP/OOSDT $ RETURN 

TSATF TCRT $ OTSOR = 0 $ RETURN $ END 


SUBROUTINE TSUMURA 


COMPARE SOUND VELOCITIES IN SINGLE PHASES. 
IDes (9) TSUMURA (NBS» FEBes 1976).- 


COMMON/1/GK, AL 9885BEyGAsDEsEPs DCRT,TCRT»9PCRT »OGAT,0TRPsTTRPsPTRP 
COMMON/7/ TAsTBsPBeRBs EAsHAsSAsy E£B9HBsSB 

COMMON/8/ INoIKy PoToDENs EgHoSs CVsCPsCSATs WoWK 

FORMAT(I5, 3F10-0) 

FORMAT(1Hi 7X *TSUMURA ETHANE SPEED OF SOUND, METERS/SECe * “4/ 


4 8X 2HID 7X3HT5K SXSHP»yBAR SXSHMOL/L SXSHW,XPT SXSHCALCO 4&XGHPCNT) 


FORMAT(5X I59 3F10035 2F1i0el, F8e2) 
PRINT 2 $8 N=O_ $ ODO 20 J=1,3500 
READ 19 IDs T»Ps WX $ IFCID) 8530 
N= NtL 8 IF (NeS5T254) 9511 

N= 0 $ PRINT 2 

IF(ToLT.TB) 12518 

IF(TeLT0TA) 135914 

CALL LIQTHRM $ GO TO 19 
IF(TeLTeTCRI) 15517 

PX = PSATF(T) $ IFCPeGTPX) 16,18 
CALL THRMOX $ GO TO 19 

PX = PVTF(T,OCRT,0) $$ IF (PeGToPX) 16518 
CALL THRMOGS 

PCT = 100.0*(1-W/WX) 

PRINT 35 IDs TsPoDENs WXoWs PCT 
RETURN $ END 
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SUBROUTINE WSATLIQ 
C IDeg (1)POOLEVAZIZ, (2) VAN GEEL, (3) TSUMURA. 

COMMON/7/ TAsTBePBeRBy EAsHAgSAs EB,HBsSB 

COMMON/B/ INoIKo PoTsDEN»y EoHoSe CV>oCPsyCSATs WyWK 
1 FORMATC(I5, 2F10.0) 
6 FORMAT(1H1 7X *SPEED OF SOUND, SATeLIQe ETHANE, METERS/SEC. * / 

1 8X *(1)POOLE/AZIZ, (2) VAN GEEL, (3) TSUMURA. * // 

4 8X 2HID 7X3HTsK SXSHP»y»BAR SXSHMOL/L SXSHW,XPT SXSHCALCD 4&X4&HPCNT) 
7 FORMATI5X 15,5 SF1IG~e3, 2F1IDelsy F822) 
16 PRINT 6 sy ib Se $ DO 20 J=1,300 
414i READ 1,5 IDs Ts WX § IF{ID) 12,30 


42 LC =Lti BS IF (LoGTe40) 13514 
13 = 0) SRR TNT NS 
14 P = PSATF(T) $ IFC(T.GT.TA) 15516 


15 CALL THRMXL $ GO TO 17 

46 CALL SATLQTR 

17 PCT = 100.0*(1-W/WX) 

23 PRINT 7» IDs TsPyDEN> WXoWsPCT 
30 RETURN $ END 


FUNCTION XBF(T,D) 
C XBF = SQRT(T/TCIFLNIT/TS) = Q(T) *Z(RsT) 9 
C Z{RoT) = LN(U), UCR»T) = T/TSIR) © 
COMMON/1/GKs AL sBBy BE sGAyDEsEPs DCRT,TCRT yPCRT yDGAT»OTRPsTTRP »PTRP 
COMMON/3/0PDT yD2P0T2,DPSOT,OPMDT,DP0D,0P [R»DTSORs,DTHOR 
COMMON/4G/ XB15XB29 XC19XC2»y XELXE2, OXBODR»DXCOR,DXEDR 
COMMON/6/ TSATs THETA, PSAT 


1°16 = IGRI 3 pus = USAT See xe 7c 

Ssrys 8 Ube Ss wes7rs GS UlR = CUS SOW US 

3 Z = ALOG(U) $ ZiR=ULR/U $ Z1X=U1K/U_ §$ Z2X=-Z1X*Z1Xx 
4 Q = SQRT (X) $ Q1 = 065/Q $ Q2 = -Ql/2/x 

5 XBF = Q*Z $8 OXBOR = Q*ZiR $ XBi = Q*Z1X + Qi*Z 

6 XB2) = QEZ2KNe we * Oe EX se aU 2Z $ RETURN $ END 
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C 


OONnN OF WHE 


et ae we 
NH © 


bh be 
Ue WwW 


XEF 
PSI 


FUNCTION XEF(T,D) 

PSI = PSISAT, DEFINE W = (1-TH/T), WHEN = 

DE/X + (1-DE)*#(1 - W +t WFLNOW)), O tS OF 1S 1 
COMMON/1/GKy AL »BByBE yGAyDEsEP5 DCRT»TCRI»yPCRT »DGAT,DTRP»TITRP»PTRP 
COMMON/3/0P0T»,D2P0T2.0PSOT»DPMOT ,DP00,DPDR,DTSOR,DTHDOR 

COMMON/4G/ XB15X325 XC19XC2_y XE1,XE2, OXBOR»DXCDR,DXEDR 

COMMON/6/ TSATs THETA, PSAT 


IF(GeEQ-OCRT) 1592 

A=0DE $ B=i1-A $ TC = TCRT $ TH = THETA $ TS = TSAT 
X= T/TC $8 X2 = X*X $F$ XS = TS/TC $$ XSi = OTSOR/TC 
W= 1.0 - TH/T $ ITFCW) 1591555 

WiR = -OTHOR/ZT $$ WAX = TH/T/X $$ W2X = -2*W1iX/X 

G = ALOG(W) $ U= W*G => Wt ft $$ UiR = GFWIR 

U1X = GFWiX $$ U2X = G¥W2X + WIX*W1X/W 

P = A/X + BYU $ PIR = BYUIR 

XEL = -A/X2 + BYULX $$ XE2 = 2*A/X/X2 + BUCK 

WS = 12-0 - TH/TS $ IFCWS) 15515511 

GS = ALOG(WS) $ US = WS*GS - WS + 1 

US1 = (TH*DTSOR/TS |= DTHDR)*GS/TS 


PS = A/XS + BEUS $ PSL = -A*XSA/KXS/XS + BUSI 
F=P- PS $ ODXEDR = PiR = PSi $ RETURN 
XEF = XE1 = XE2 = DXEDR = 0 $ RETURN $ END 
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SUBROUTINE FITTEX 
COMMON/3993/ NCOF 5Y 5639) 
DIMENSION A(3u,31),8(30,31) 
TYPE GOUBLE SYsSYYs,RES A588 
EQUIVALENCE (A,B) 

DATA (NTR==-1) 5 (NGIM=3~.) 

37 FORMAT(*1 THE COEFFICIENTS ANG THEIR ESTIMATED ERRORS ARES*/7) 

38 FORPMAT(*S F/*5 */*S ESTIMATCO R=ESTOUAL SUM OF SQUARES =*£17.9/ 

41 * ESTIMATED REGRESSION SUM OF SQUARES =*E147.9/ 
2 * =STIMATED TOTAL SUM OF SQUAFES =*¥617.9/ 
3* VARIANCE OF FIT =*£417.9/7* DETERMINANT OF THE MATRIX =*E17.9/ 
&* CORRELATION COEFFICIENT =*e1i7.9/* NUMBER OF PCINTS =*15) 

45 FORMAT(*1 THE ARRAYS IN THe FITTING PROGRAM ARE TOO SMALL TO HOLG T 
tHE NUMBER OF CONSTRAINTS AND FUNCTIONS ASKED FOR IN THE CALLING PR 
20GRAM*) 

SHAS TE ORM ATG SANS ele Gree tO eraser) 
ENT SR HESS SOM Bitie THe DATA 
ENTRY FLT 
TFCNTF) 159353 

41 NP=e 
NF=NC OF 
IF(NF.eGTeNDIM) GO TO 44% 
NCON= }} 
SY=u56 
SYY=i. 
NY=NF +14 
09 2 IT=1,NY 
CO 2 J=1,NF 
2 ANGI =r. 

IFC(NTFLEQe9) GOTO Lt 

NTR=8 

S/S ar Si 

SYY=SYY#Y*Y 

DO & J=1,NF 

ACJSosNYI=ALTI NY) FY #605) 

90 4&4 J=i,NF 

G4 ALT ,J)=A(L,J) +G(1)*G(J) 

NP=NP41 

-ETURN 

ENTER HERE TO CONSTRAIN THE EQUATION 
ENTRY CONSTR 

IPN Ur) Bb it} aeab setae 


Gl 


{iz NTR=2 
CO YO 3 
44 N=NY-1i 


IFC{INY#NCON+2).GT».NOIM) GO TO 44 
DO i2 I=1,N 
A(I,NY+i)=A(I,NY) 
A(NGIM=NCON 51 )=5(T) 
A(NY,1)=G(T) 

42 A(I,NY)=G(T) 
NCON=NCON+t+4 
N90 13 I=NF,N 
A(NY,I+i?=2.) 

13 Al(I+1,NY)I=2.3 
NY=NY +1 
RAINY=1,NY)=Y 


44 


Appendix D. Continued. 


R=TURN 
EN Pee REPS FO UNVeRY MANS AND) Gey WW Rie sGishys 
ENV RAG GOEIAlF 
N=NY=4 
BO) 2e L=19NF 
230 G(I)=A(I,NY) 
UD 22 I=25N 
DO Zt Deaton 
21 ACT=1,I)=A0T=-1,))/A(T-1,1-1) 
DO 22 J=I,N 
DO) 22. Weal ony 
22 AMI gKIV=HA0CS, KI HAC I, I-41) FA(TH1,K) 
AIN,NYI=ACINNYISA (N,N) 
DOTS P= 2 5)N 
L=N-I+2 
CO 24% J=LyN 
24 Al(L-1isNYI=A( Lat sNYI-A(L=1,/5)*ATI,NY) 
RES S51 Vf 
v9 25 T=1,NF 
RP=S=RKRZS“ACI NYID*GCTY 
25 Gt(TI=ACI,NY) 
DF=NP-NF+NCON 
Y=NCON 
NTR=-1 
R= TURN 
ENee RENE [ROM Swipes ywGs OlF CORR RGIS 
ENTRY STAT 
6G ef T=1,NCON 
[0 27 J=1,NF 
27 R=S=PESSA(NOIM=-IT 44155) *A(J,NY) *ACNFHTI, NY) 
TO —SYV= SV SY ANP 
REG=TOT-FES 
SYY=8t S/0F 
ST=1. %6+2.72/0F+8 .04/0F**3 
DeET=1. 
LO 32 I=i,NF 
Dz Sieh (0g ag) 
Su A(T ,T)=1-5/Ai 1) : 
SO 32 L=2,NF 
C0 S2 Jeeou 
SY=C. 
DO) Gel Kew} 5 ae 


Ree BS) ES ead Wa OT ts Wi) ia 0 Ged GE I 
fd A(T, J-1)=SY*A(I,T) 

IPRIUINL Be 

DO 36 I=1,NF 

L=NF-T 

DO $3 J=1,yL 

K=NFeJ 

NO £3 M=1,J 

N=NF-F44 


ZZ AIK,I)V=A(KyI) -ACKyN) *L(N,T) 
00 $4 J=2,l1 

$4 A(S-1,1I)=A(1T, J-1) *SYVY 
YO SS Y=ls. 

35 A(T ,JI=A(I,J) *SYY 
68=ce(1,1) 
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Gy 


BS US ile WARTANCe OF VAS CCSFFICIENYS 
LE Ce Bie: «ul B= — 3 

PRES) SOR) CB 3) 

32R=R(1T,NY) 

|= Bat fh Dn 7a Wt Je LG GS la 

LE CSV VokUotot) si VecSvyiv 
COL{LSESE/7 TOT 

BDOTNT S83 =2S, 525, l0},SYY, 087 W602, NF 
Sex = RES bh Vf = SMe Y (SE R/ GF) 

fC eTURN 

PEINT 45 

STOP 

ENO 
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anNgan 


a 
2 
5 
10 


12 
14 
15 
16 


17 
18 
19 


PROGRAM ETHRM8A (INPUT,OUTPUT) 
ETHANE THERMOFUNCTIONS, BRITISH UNITS, REVISED NOVe 20,5 1976.6 
ENTER COMPR. LIQ. STATES VIA NBS(RODER) CY ON ISOCHORE RB, AND CSAT. 
NOTE THAT TSATF (DEN) ITERATES ON DENGASF(T), DENLIQF(T). 
NOTE TTRF = -297.0436 Fy, TCRT = 89.9240 F. 
NOTE DCRT = 12.7650 LB/CUFT, PCRT = 706.536 PSI. 
COMMON B1,82583,5 CisC 2,03, EL1,E25E35 ER 
COMMON/1/GK,AL5BB;BE »GA,DE,EP, DCRTsTCRT,PCRT ,DGAT,OTRP»TTRP»PTRP 
COMMON/3/DPOT ,D2P0T2,0PSDT»O0PMDT ,DPDD,DPDR,DTSOR,DTHDOR 
COMMON/4/ XB15XB25 XC19XC25 XEL,XE25 DXBDR,DXCDR,DXEDR 
COMMON/5/ OOSDT 
COMMON/6/ TSAT»s THETA, PSAT 
COMMON/7/ TAsTBsPByRBs EAsgHAsSAy EB HB 9SB 
COMMON/8/ INygIK»o PoToDENo EsHoSo CVoCPyCSAT> WyWK 
COMMON/9/ AU,BUsCU,DU, AS,BS,CS,DS,ESy ETBySTBySTA 
COMMON/11/ JT, TR, TF 
COMMON/99/ TIsEZZ, EZsSZ,CVZ5 HZ,CPZ 
DIMENSION PSI (99) 
FORMAT(I5, 2Fid.0) 
FORMAT(I5, 3F10.0) 
FORMAT (1X) 
FORMAT(1H1 9X *ETHANE FUNCTIONS ON ISOCHORE RB = 14¢5 MOL/L* // 
1 10X 7HAT TB =F8e35 6Hy PB =F8e35 7Hy DEN =F723// 
2 10X 4HEB =Fi0e2s 6H, HB =Fi0e2s GH, SB =F1iDe4// 
3 10X 4GHCV =F8.35 10H» CVyXPT =F8e3/// 
4 10X 7HAT TA =F9e4%s GHs PS =F8e3_9 THe DEN =F7.3// 
5 1GX4HEA =Fi0.2, GH, HA =F1I0e2s GHy SA =F1I004//10X 4SHCV =FBe3S/// 
6 10X SHETB =E172-95 7H», STB =E17-¢9y 7H, STA =E17.9) 
FORMAT (8F10.5) 
FORMATC(1HL 18X *ETHANE ISOBAR AT P =* F9.5y 4&H PSI / ) 
FORMAT¢1H1 18X *ETHANE ISOBAR AT P =* I6y 4H PSI / ) 
FORMAT(14X 1HT 8X1HT EXSHOEN OXSHVOL 4XSHDP/OT 6X5HOP/00 
2 7X1iHE 7X1LHH BXLHS 7X2HCV 7X2HCP 5SX1HW / 
3 10XSHDEG F &XSHDEG R 2X7HLB/CUFT 2@X7HCUFT/LB 4&XSHPSI/R 
& AX10HPSICUFT/LB 2X6HBTU/LB 2X6HBTU/LB 1X8HBTU/LB/R 
5 1X8HBTU/LB/R 1XBHBTU/LB/R 2X4HFT/S ) 
FORMAT (OX 2F90.3_9 F9e35 F905 FIeSs File2sy 2F8e29 F905, 2F903y 16) 
FORMAT (6X 2F9.39 F9e4s FIe4s FI04s File2s 2FBe2y F905y 2F923516) 
FORMAT (6X 2F9e35 F905s F903 FI05s File2sy 2FBe2y F905,y 2F903,I16) 


CONSTANTS FOR EQNSTATE MARK VII, 4/12/76. 
WM = 30-07 $ WK = 100000/WM $ GKK = 0.0831434 
Q = 1.01325 $ QP = 14.69595/Q 


TTRP = 90-6348 $$ OTRP = 21-68 $ PTRP = 1.130842618E-5 

TCRT = 305.33 8 OCRT = 6680 $ PCRT = PSATF(TCRT) 

DGAT = 1.5153661E-6 $ GK = DCRT*GKK $ EZZ = 4.18638*4827.2 

AL =1 $ DE = 2.9/3.0 $ ER= 2 $ BB=BE=GA=EP=6 

B1 = 0248752227313 $ B2 0-33198750982 $ B83 0.06854249828 


El =-02043113918548 3$ E2 0.0 $ ES 0.03779460468 
INTEGRATION CONSTANTS FOR CVSUMIT, SSATF(T). 

TB=330 $ RB=14.5 $ TA=TSATFC(RB) $ X=TB/TCRT $ G=ALOG(X) 
AU=64.0253 $ BU=-163. 36173 $ CU=205.55275 §$ DU=-54.77598 

ETB = AU*G + BU*K + CUFX*FX/2 + DUFX*X*¥X/3 

STB = -AU/X + BU*¥G + CU*X + DU*X*X/2 

AS=22.64822683 $8 BS=0.08796160711 $ CS=1.090640627 
0S=5371.201054 $ ES=-212338.994% $ T=TA $ TDO=C(TCRTI-T)**0.4 
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36 STA = AS*ALOG(T) + BS*#T = CS*#TD/0-4 = DS/T = ES/2/T/T 


C ESTABLISH TBsTA»y EBsSBe EAsSAy ON ISOCHORE RB = 14-5 MOL/L. 
39 T= TB $$ P = PB = PVTF(T,RBs0) $$ CALL THRMOGS 
40 RC=DEN $ EB=E $ HB=H $ SB=S $ CVB=CV $ CVX=CVXF(T) 
41 T= TA $ CALL CVSUMIT 8 EA =E $ SA=Z=S 
42 PSTA = PS = PSATF(TA) $ HA = EA + 100*PS/RB 
43 PRINT 105 TBsPBsRCy EBysHBySBy CVByCVX, TAsPSsDENy EA sHAgSAsCVy 
1 ETBy STBy STA 
NOTE RESULTS FGR COMPUTED DATA = 
AT TB = 330 Ks RB = 14.5 MOL/L, PB = 422.095 BAR, 
EB = 21321.0% J/MOLy, SB = 16504420 J/MOL/K,y 
CV»CALC = 566108, CVsyXPT = 54.053 J/MOL/Ke 
FOR CVSUMIT, cTB = =-74.58239 J/MOL/Ky STB = 1186423635 J/MOL/Ky 
AT TA = 256.5624 Ky RB = 1465 MOL/L» PSTA = 152616 BAR, 
AND, FOR ESATF(T)s SSATF(T)s EA = 17628286 J/MOL, 
SA = 152.8227 J/MOL/K,y STA = 115.97498 J/MOL/K. 
40) CALL JTLOCUS $ CALL TABLIQ 


AINNDAMNND 


COMPUTE THERMOFUNCTIONS ON ISOBARS. 
EACH ISOBAR STARTS ON THE MELTING LINE. 
ITSOBARS AT P UNDER PCRT TRAVERSE THE DOME. 
110 IN=1 $ NI = 69 $ READ 125 (PSI(I), J=1,NI) 
111 00 300 I=IN,NI $ IK = IT\sSe  EPSTV=FPShCI) SP p= ePShOn Aer 
112 IF(I.EQ.4) 113,114 
113 PRINT 14, PSI(I) $ GO TO 115 
114 PRINT 15, IPSI 
115 PRINT 16 
120 T = FINDTMF(P) $ CALL LIQTHRM $ CALL CONVERT $ V=1/DEN $ IW=W 
122 PRINT 17,5 TFsTReDENsVs ODPOT,O0PDD, ExyHsyS»y CV,CP,IW 
123 IT = TF/20 - 1 $ IF (PeLTePCRT) 125,200 


qgqngnggn 


TWO MAIN CASES FOR P UNDER PCRT. 
NOTE 20 DEG.F STEPS FROM =280 THRU 800 Fy (55 POINTS). 
125 TS = FINOTSF(P) $ TX = TS + 201-8 $ K=L= 0 
130 DO 199 J=1,99 $ JT = 20*(IT + J) 
131 T = TT =-273.15 + (JT=32)/1.-8 
132 IFCP=PSTA) 133,133,150 
C MAIN CASE A FOR P NO GREATER THAN PSTA. 
133 IF(T LT.TS) 134,136 
1434 CALL LIQTHRM $ CALL CONVERT $ V=41/0EN $ IW =W 
135 PRINT 175 TFoTReDENoVs DPDT,»,DPDD», EsHsSs CV>CP,IW $ GO TO 199 
136 IF(T LT-TX) 137,180 
137 T= TS $ CALL SATLQTR $ CALL CONVERT $ vV=1/DEN $ IW=W 
138 PRINT 175 TFsTReDENsV>» OPDT,DPOD, ExsHoSy CVsCP,IW $$ PRINT 5 
139 CALL SATGSTR $ CALL CONVERT $ V=1/0EN $ IW=W 
1440 IF(P.GT.10.9) 141,142 
141 PRINT 189-TFoTReDENoV>o OPDT,DPDD, EsyHsSy CV>oCP,IW $ GO TO 143 
1442 PRINT 195 TFoTR»eDENsV5 DPOT,»,OPDDs EsHsSy CV,CP,IW 
HS 1 SOUL $ GO TO 180 
C MAIN CASE 8 FOR P GREATER THAN PSTAe 
1450 IF(T-TA) 151,151,153 
151 CALL LIQTHRM $ CALL CONVERT $ V=1/DEN $ IW=EW 
152 PRINT 175 TFeTR»DENoVy DPOT,0PD0D, ExHsS»y CVsCP,IW $ GO TO 199 
153 IF(TeLTeTS) 154,156 
154 CALL THRMOX $ CALL CONVERT $ v=1/DEN $ IW=W 
155 PRINT 175 TreTReoDENoVs DPDT,DPO0Ds EsHsS»y CVsCP,IW $ GO TO 199 


aan 
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1560 IF(T.LT.TX) 15751560 

157 T= TS $ CALL THRMXL $ CALL CONVERT $ V=1/DEN $ IW=W 

158 PRINT 175 TFaTRoDENoV> DPDT»,DPDD»y EsH»Sy CV>CP,IW $$ PRINT 5 

159 CALL SATGSTR $ CALL CONVERT $ V=1/DEN $ IW=W 

160 PRINT 18,5 TFseTR»sDENeoV> DPDT,DPDD, ExHsS»y CV,CP,IW $ T = TT 
C CASES FOR HOMOGENEOUS DOMAIN», DEN S$ DCRT, OR T 2 TBe 

180 IF(JT.GT.600) 181,183 

181 K = KtL $ JT = JT + 20*K $$ ITF(UT¢GT 800) 3005182 

182 T = 273.15 + (JT-32)/1.8 

183 CALL HOMOTRM $ CALL CONVERT $ V=1/DEN $ IW=W 

184 IF(P.GT.1i0) 185,186 

185 PRINT 185 TFeTReDENosV>y DPDT,»,DPDDs EsHsS»y CVsCP,IW $$ GO TO 199 

186 PRINT 195 TFeTReDENsV> DPOTsOPDD» EsHsSy CV,CP,IW 

199 CONTINUE 


C CASES FOR P > PCRT. 
200 K =L=0 $ ODO 230 J=1,99 $ JT = 20*1(IT+IJ) 
201 T = 273415 + (JT-32)/1-8 
202 IF(T.LT.1TB) 203,220 
263 IF(T.LT.TA) 204,206 
204 CALL LIQTHRM $ CALL CONVERT $ V = 1/0EN $ IW =W 
205 PRINT 179 TFeTReDENgV,DPDOT,DPDD,E,HsSy,CV5oCP,IW $ GO TO 230 
206 IF(TeLTeTCRIT) 211,207 
207 PX = PVTF(T,OCRT,0) $ IF(P.GT-PX) 211,225 
Cc CASE FOR REGION X AT TA < T < TBy P ABOVE CRITICAL ISOCHORE. 
211 CALL THRMOX $ CALL CONVERT $ V = 1/DEN $ IW =W 
212 PRINT 17% TFsTReDENeoVsDPDT »sDPDDsEgHsSsCVe,CP,IW $ GO TO 230 
Cc CASES FOR HOMOGENEOUS DOMAIN, DEN ¢ DCRT, OR T 2 TBe 
220 IF(JT.GT.600) 221,225 
221 K = K+i $$ JT = JT + 20*K $$ IF(UTGTe800) 300,222 
222 T = 273-15 + (JT-32)/1.8 
225 CALL HOMOTRM $ CALL CONVERT $ V=1/DEN $ IW =W 
226 PRINT 17,5 TFeTReDENoV>sDPOT ,DPDD oE »HoSoCV5oCP,IW 
230 CONTINUE 
300 CONTINUE 
999 STOP $ END 


SUBROUTINE CONSAT 
C CONVERT SI TO BRITISH UNITS» ETHANE MW = 30-07 G/MOL. 
COMMON/3/DPDT ,02P0T2,0PSDT,DPMDT ,OPDD,0PDR,DTSOR,DTHOR 
COMMON/8/ INsIKe PoTsDENy EsHoSo CVoCPsCSATs WoWK 
COMMON/11/ JT, TR, TF 
COMMON/412/ TT5PSsDL> ODLOTs, VGyVL5Q, QX,CSX 
DATA (EK=69. 8960223) 5 (SK=125. 612840) » (DK=0253270577) 
DATA (PTK=0.12410562) , (POK=0.129428985) 
TR = 1e8*TF § TF = 32 + 168¥*(TT-273.15) 
PS = PS*¥14&.503775 $ DPSOT = DPSDT/PTK 
DPOT = DPOT/PTK $ OPODD = DPDD/POK $ OOLDT = DOLOT/DK/1.8 
DL = DL/DK $ VG = VG*DK $ VL = VL*DOK 
E=E/EK $ H=H/EK $ Q=Q/EK $ QX=QX/EK $ CSX=CSX/SK 
Ss S/SK $ CV = CV/SK $ CP = CP/SK §$ CSAT = CSAT/SK 
W 3.280 84*W $ RETURN $ END 


ONNEWNHFE 
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C 


OS wnere 


NaN00D 


6 CK(J) 


ait 
2 
3 DK{(J) 
4 
5 


SUBROUTINE CONVERT 
CONVERT SI T9 5RITISH UNITS, ETHANE MW = 30207 G/MOL. 
COMMON/3/0PO0T ,02PGT2,0PSDT,DPMOT ,OP0D,D0PER,0TSOR,OTHOR 
COMMON/8/ INsITKy PoTsDENsg EsHsSo CVsCPoCSATs WoWK 
COMMON/411/7 JT, TRo TE 

DATA (EK=69.8960223) » (SK=125. 812840) , (DK=0.53270577) 
DATA (PTK=0.12410562) » (POK=0.129428985) 

VWs Sosy 8 We Ss se @ tobe Uscasoss)) 

DEN = DEN/OK $ £& = E/EK $ H = H/EK 

S = S/SK $ CV = CV/SK $ CP = CP/SK $ W = 3.28084*W 
DPOT = DPDTVPTK $ OPDD = OPDD/PDK 

RETURN $ eND 


SUBROUTINE HOMOTRM 

VALID ONLY AT T»F = INTEGER MULTIPLES OF 20 OEGeF. 

GIVEN Pot FOR HOMOGENEOUS DOMAIN, GET OEN AND FUNCTIONS. 

SAVE OKT), EK(T)s SK(T),5 CK(T) ALONG EACH ISOBAR. 

NOUS Sealy Ss oh SS CON elas (eve) - he) S&S 5S POuNySo 
COMMON/3/0PO0T , 02P0T2,D0PSDT,0PMDT ,0PDD,DPOR,DTSOR,OTHOR 
COMMON/8/ INogIKy PoToDENs EsHsSy CVsCPyCSATs WyWK 
COMMON/11/ JT, TR» TF 
COMMON/99/ TIyEZZ,y EZsSZsCVZ_ HZ,CPZ 
DIMENSION OK(70) sEK(70) »SK(70) »CK (70) 

DATA (Q=1-91325)5(G=0-0831434) 

J = JST/20 +15 $ IFC(IKeEQ-IN) 2510 

TI = T 3S GALL IDEAL $ ODA = 0 
DEN = DB = FINDENF(TsP) $ N = OB*¥20 + 10 
= —EZZ + EZ + ESUMF(NsT,DA,DB) 

='SZ + SSUMF(1isNoT,»0A,08) - 100*G*ALOG(G*T*0B/Q) 
CV = CVZ + CSUMF (NoTsDA,DB) $ GO TO 20 


EK (J) 
SK (J) 


How ww 


10 DA = DK(J) $8 OK(J) = DEN = DB = FINDENF(T,P) 
42 DX = ABS(CB+GA) $$ N = (DX+1)*5 


C 


ot 


WFAN 


SE 
= § 

= CV = CK(J) + CSUMF(N,T,DA,0B)~ 

+ 100*P/0B $ PX = PVTF(T,DB, 1) 

V + 130*T/DPOD* (DPDT/DB) **2 

RT (WK*CP*DPDD/CV) $ RETURN $ END 


EK(J) + ESUMF (N,T,DA,0B) 
SK(J) + SSUMF (05N5T»0A, 0B) 


SUBROUTINE JTLOCUS 
DERIVE THE J-T INVERSION CURVE. USE ROUTINE DELTAF(T,DI). 
DIMENSION 17 (99) ,TF(99) »PP(99) ,ON(99) 

DATA (DCRT=6- 80) 5(TCRT=305-33) 

FORMAT(1H1 16X *THE JOULE-THOMSON INVERSION LOCUS FOR ETHANE*// 
1 12X SHTsF SX3HTyR 2X5HP,yPSI 2X7HLB/CUFT 
2 9X SHToF SX3HTsR 2XSHPyPSI 2X7HLB/CUFT ) 

FORMATC3X 2(F12e2, FBe2s I7y F902)? 

PRINT 1 $ NP = 64 

DO 25 I=1,NP $ TF(I) = -20 + 10*I 

TT(I) = 459.67 + TFC) $ T = TT(1I) 71-8 
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Appendix D. Continued. 


6 X = T/TCRT $ OX = 166 
9 DI = DERT*(2e40 - 0658*X + De 2usx) 
10 IF(T-TCRT) 11,12,12 
142 OL = DENLIQF(T) $$ IF (DI=DL) 25,12,12 
12 SS = DELTAF(T,DI) $ DO 20 IT=1515 
44 D=DI-OX $8 SL=DELTAF(T,D) $ D=DI+DX $ SP=DELTAF(T,D) 
1S IF(SS=SL) 18316,16 
16 IF(SP-SL) 13517,17 
iio Sl Ss: OL = DE = DX 1S GOMTO 20 
18 IF(SS-SP) 20,20,19 
19 SS = SP 3 ODI = DI + DX 
20 DX = OX72 $ ONCI) = DI/0-53270577 
21 PP(I) = 14.503774*PVTF(T,01I,0) 
25 CONTINUE $ N = NP/2 $ ODO 29 J=1,N 
26 IP = PP(J) $ IPP = PP(J+N) 
29 PRINT 25 TF(J) oTT( J) os IPyDN(J) 5 TFECUtN) o TT (UN) gp IPP, ON (UN) 
30 RETURN $ END 


SUBROUTINE TABLIQ 
C TABULATE THE ETHANE SATURATED LIQUID FUNCTIONS. 
COMMON/1/GK,AL»BB,BEsGAsDE,EPy DCRT,TCRI sPCRT »DGAT » DT RP» TTRP»PTRP 
COMMON/3/0P01T ,D2PD0T2,D0PSOT ,OPMOT ,DPDD,DPCR»,OTSDR,DTHOR 
COMMON/5/ DOSET 
COMMON/7/ TAsTBsPBsRBe EAsHAsSAy EB,HB>SB 
COMMON/8/ INeoIKs PoTsDENs EsHoSs CVeCPsCSAT>s WoWK 
COMMON/11/ JT, TR, TF 
COMMON/12/ TT,PS_yDL, DOLOTs VGsVLsQs QX_CSX 
DIMENSION TSA(46), PSA(46) 
1 FORMAT(1H1 13X *#PROPERTIES OF SATURATED LIQUID ETHANE* / ) 
2 FORMAT @(14X4HT 10X1HP SXSHDEN 4&XSHVyLIQ 6XSHV,GAS SXOHDPS/DT 
1 SX6HDDL/JOT SXSHOP/OT 6XSHOP/DD 3X5HQ,VAP 2X5HQ,XPT / 
2 10X5HDEG F 8X3HPSI 1X7HLB/CUFT 2X7HCUFT/LB 4X7HCUFT/LB 6XSHPSI/R 
3 IX9HLB/CUFT/R SXSHPSI/R 1AXLOHPSICUFT/LB 2X6HBTU/LB AX6HBTU/LB ) 
3 FORMAT(5X FiGeSs EL11¢4%, F8e3y FIe6y 2EALL]e4%s FID06, FID3, 
1 Elicets F8e2y Fred) 
4 FORMAT(14X1HT 8X1HT 10X1HP 8X1HE 8X1HH 8X1HS 7X2HCV 7X2HCS 7X2HCP 
1 6X1HW &XOHCS,XPT / 
2 10X5HDEG F &X5HDEG R 8X3SHPSI S3XOHBTU/LB 3X6HBTU/LB 
3 iX8HBTU/LB/R 1X8HBTU/LB/R 1X8HBTU/LB/R 1X8HBTU/LB/R S3X4&HFT/S 
& 2X8HBTU/LB/R ) 
5 FORMAT(OX 2F9.35 EAle4sy 2F9eS5 F9055 3FI03, I7y F1I0.3) 
C FOR PAGE ONE OF TABLIQ. 
C NOTE, TTRP = -297.0436 F, TCRT = 89.9240 Fe 
10 PRINT 1 $8 PRINT 2 $ NP = 39 
11 DO 25 J=1,NP $ IF JcEQe1) 12513 
12 TT = TTRP & GO TO 17 
43 IF(JeEQeNP) 14516 


144 77 = TCRT $$ OG = DL = DCRT $ ODOLOT = 0 

15 VG = VL = 1/0G $ GO TO 19 

16 FT = -300 + 10*J $ TT = (FT-32)/1-8 + 273.15 
17 DL = OENLIQF(TT) $ ODOLODT = DDSOT 

18 DG = DENGASF(TT) 3$ VG = 1/0G $ VL = 4/0L 


19 TSA(J) = TT 8 PSACJ) = PS = PSATF(TT) 

20 Q = 100*TT*OPSDT*(VG-VL) $ QX = QVAPXF (TT) 

21 PX = PVTF(TT,DLo1i) $ CALL CONSAT 

25 PRINT 35 TFsPS» DLoVL>VG5 DPSDT,DODLOT, DOPODT,DPOD, Q,QX 
C FOR PAGE TWO OF TABLIQ. 
C USE SATLQTR AT T < TA, USE THRMXL AT T > TA 

30 PRINT 14 $8 PRINT & 

$1 DO 40 J=1,NP $ TT = T = TSA(J) $ PS = P = PSA(J) 

32 CSX = CSATXF(T) $$ IFCTLTeTA) 34,35 

34 CALL SATLQTR $ GO TO 36 

35 CALL THRMXL 

36 CALL CONSAT $ IW = W 

40 PRINT 55 TFsoTRoPSs EsHsSy CVsCSAT,CP, IW, CSX 

99 RETURN $ END 
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Figure 1. 


TEMPERATURE 


Generalized density-temperature phase diagram, 
paths used for thermal computations. 
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Figure 2. Deviation plots for densities. 
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Deviation plots for densities, continued. 
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Deviation plots for specific heats. 
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Figure 3. Deviation plots for specific heats, continued. 
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Table 1. Experimental and calculated vapor pressures. 


Data sources for this table are-- 


ID Author Reference 
A Pal [61] 
7 Ziegler [93] 
9 Pope [61] 
10 Douslin [18] 
90 Straty [79] 
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Table 1. Experimental and calculated vapor pressures. 
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Table 2. Experimental and calculated saturated liquid densities. 


Data sources for this table are-- 


ID Author Reference 
5 x 

10 Douslin [18] 
11 Sliwinski [75] 
72 Chui [14] 
13 Klosek [45] 
14 Rodosevich [67] 
16 Tomlinson [85] 
20 Haynes [37] 
21 McClune [52] 
22 Orritt [58] 
23 Gugnoni [33] 


*Via our vapor-pressure equation, and isochores of Pal [61]. 
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Table 2. Experimental and calculated saturated liquid densities. 


SATO. LIQUID DENSITIES, F = . 353 


TTRP = GC2348, TORT = 365233), DCRT = O68c0, OTRP = 21-6084 


0 7196 845E1 ec 81866182 2289937306 
I0 WT T»K MOL/L CALCL PONT JOS/O0T 
14 1 325.2150 72827 7.316 215 -1.9553 
16 1 3052150 7.830 7.316 Outs -1.9.559 
11 1 354.65% 8.397 8-396 eul -.8t5u 
11 al Su4e15f 8732 8.736 eu2d 2567) 
41 at 303.166 9.196 92136 eve -2 3875 
11 1 302.16 953% G$e.531 oc -2e 5071 
ake 4 3022150 9-544 G.534 elt 2 9J65 
abot 1 356 e665 9.937 9-939 =o 2 -ec4#31 
1€ 0 320.600 9-989 92939 eke -e c4S1 
11 1 298.175 19.472 1c 2468 euh -21887 
di 1 298 .15e 10.499 e472 eco -21i884 
16 ) 298.17¢ 10.502 1L.408 032 -.1887 
11 at 295.67) 103.399 12.3939 -eud -.21585 
16 g 2952675 13.925 15.8939 023 20 SIs) 
abel 1 293-2180 11.c7a 14.266 Ou -213389 
16 G C2995 a oe 11.2950 11.258 02d -2ei589 
if at 295 e15u 11.297 11.272 ede 21387 
23 0) 283.150 12.438 12-.4%c8 -.16 -2t 968 
141 1 28522235 122423 122423 ~eif -.C990 
16 sy 283-29 120434 120423 ouJ -243990 
i¢ 1 233.150 12.458 1224cb 02h -.°988 
1¢ 1 272.1 51 13.342 13.315 ecu -2 6804 
g 3 Cee CLG 13.549 13-502 © e539 -~oS 764 
20 3 26c 0958 144.261 14.275 onan 216638 
2s ry 255 e372 14.473 14.576 -e71 -. 5634 
5 c 2552963 14.554 14.538 edi -2638 
1c a 2552150 14.753 14.715 220 -.C62u 
20 3 254 0 BUG 14.899 14.3938 -eid -evoGl 
5 s) 2472962 15.C50 15.929 eid =. 1596 
23 U 2412483 152.344 15.400 -.37 -25558 
5 J C45 eal 152455 15.446 oul 20555 
2c 3 240 ecu 15.464 152483 Coie =.055 
22 i Z282okur 15.904 15.906 -e Id -eJ5c1 
- 26 3. -250.05¢ 15.997 1o.u14 SO ebah -2t514 
5 g 229.917 16.037 16.019 ei -2u513 
Ze 1 2252058 16.220 16.232 =ef@ -. 0508 
~Bon ne --  -2 22 2 48 16-423 16-385 ec3 ~e 3491 
2c 8 220.950 10.499 16.512 -.t8 2 0484 
5 0 2142942 16.754 16.754 0 lS ~et471 
22 1 212-6e1 16-836 16.862 -.16 -20465 
2c 3 210 .de0 16.9714 162984 etl ei 468 
5) 3 2572941 17.1¢5 17.078 oc& -2l455 
22--. 4 203.191 417.275 172232 ~eit -e 0446 
5 C 1972888 17.529 17.5c6 232 -2 0436 
22 1... .192.-723 17.819 172834 -2f8 -~204%25 
20 3 19C.d0CC 17.861 17.305 -.C2 =o 0424 
5 t 138.451 17.941 17.93 06 -e 2422 


63 


Table 2. 


Continued. 
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Table 3. 


Experimental and calculated saturated vapor densities. 
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Table 4. Experimental and calculated second virial coefficients. 


ID: (3) Michels; (6) McGlashan; (10) Douslin 


GC ON Gs (Or Or Se 


6.552071 -41.106244 -0.5°2947 -0.041944 0.000000 
) Folk WAYS 3+ GALE DIFF ECNT 
LS 6 a 6.4912 -5.183 -5.185 0.002 6.03 
160.930 6.524) “4.469 -4.437 -0.001 -0.03 
170.090 0.5567 -3.2°35 - 3.933 -0.002 -0.04 
136.998 0.53854 - 32485 ~- 32483 -0.001 -0.04 
130.906 0.6222 -3.112 S_goliZ -0.001 -0.02 
206% UIE 6.6549 -2.°00 -2.800 0.000 0.01 
210.099 6.6&/77 =-22525 -2.536 0-001 0.04 
BS) 226.610 067202 -2.3507 -2.309 0.002 0.07 
b 230.5)6 057522 -2.110 -2e111 0.002 0.09 
) 246.000 0.7653 -1.°937 -1.939 0.002 0.10 
6 250.0)% Wa taal ir/ “1.765 -1.766 0.002 0.10 
“9  @26Go0I2 6.65144 -1.650 -1.651 0.0014  @owel 
270.999 0. 6842 -1.529 -1.530 0.004 0.05 
© PSV 5 13)5 0.9163 -1.421 10421 9.000 0.00 
b 290.606 0.9497 -1.523 -1.323 -J.001 -0.06 
5 390.090 0.9824 -1.235 -1.233 -0.002 -0.14 
3 273-150 0.6945 -1.403 -1.494 0.002 012 
1D 275219% 0.6945 -1.498 -1.494 -0.003 -0.21 
3 236.135 0.9703 -1.251 -1.249 -0.002 -0.14 
1 a 239¢.2153 G.c7e4 -1.252 -1.249 -9.003 -0.26 
1u 303.150 0.9927 -1.299 -1.207 -0.002 -0.21 
3 3220045 1.0563 -1.058 -1.058 -0.000 -0.01 
16 323.150 ‘ko G52 -1.056 -1.055 -0.001 -0.14 
3 347.952 1.1385 -0 2.898 -0.900 0.002 018 
10 BA 6 51559) 12.1401 -9.896 -0.597 0.001 0.07 
3 SV Acd22 1.219) -0.769 -0.770 0.001 A oil 3 
10 Siaciolioll Lo 222) -0.766 -0.767 0.004 014 
3 337.544 » 3023 -§.659 -0.659 -9.900 -0.02 
19 336.159 kG BSS) -GFFE -0.657 0.001 G23 
3 KER oT DO Lo BRA2 -£.566 -0.566 0.000 0.02 
10 423.156 163657 -0.563 -0.564 0.001 0.12 
10 446.1350 1.4075 -004E3 -0.484 0.001 0.13 
16 473.156 1.549+ 00415 - -00414 -0.000 -0.09 
10 49é.156 G6 BGLS -0.353 -0.353 -0.000 -0.07 
10 623.150 Lo PLS2 -0.300 =0.299 -0.001 0.37 
10 Bt 158 1.795) -0.251 -0.251 -0.001 -0.30 
10 673.159 1.0763 -3.cC8 -0.208 -0.001 -0.36 
Lu £396021535 1.9503 -0.168 -0.169 0.000 0.10 
1) 623.153 2.0405 -0.132 -0.134 0.001 1.06 
NO S°S3o MEANPCT = fo 139s 


66 


ID: (1) Eucken; (2) Lambert; (4) Hoover; (5) Pope; (8) Gunn; (9) 


Table 4, Experimental and calculated second virial coefficients (Continued) 
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Table 5. 
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Experimental and calculated third virial coefficients. 
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Table 7. 
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Coefficients of the equation of state. 
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Table 9. Interpolated ideal gas functions. 
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Table 10. Experimental and calculated heats of vaporization. 
ID: (10) Douslin; (11) Wiebe; (13) Riedel via Furtado; 
(20) Furtado. 


012857784E+02 054095302E+01 ©33835851E+02 

-297520486E+02 ©10084976E403 -237933878E+02 
Id T9K KJ/MOL CALCD PCNT 
10 248.15 aloe 11.235 057 
10 253.15 10.925 10.863 057 
10 263.15 10.072 10.039 232 
10 273.15 9.073 9.072 2001 
10 283.15 70844 72877 - 043 
10 293.15 6.176 60.252 -1.32 
10 298.15 52034 52090 -1.10 
10 302.15 3.€77 3.710 -.90 
10 303.15 3.207 3.210 -.07 
10 304.15 20546 20544 016 
13 100.00 17.154 17.311 “91 
13 111.11 16.928 17.068 =.82 
13 133.33 16.455 16.493 -.23 
13 155.56 15.853 152797 035 
13 166.67 15.4E0 152403 250 
13 160.00 14.957 14.886 048 
12 164.10 14.70€ 14.716 -.07 
13 184.11 14.802 14.716 059 
13° 190.00 14.52E 14.464 043 
13 200.00 14.049 14.012 027 
13 210.00 13.530 13.524 2°04 
13 220.00 12.999 12.998 204 
13 230.00 12.358 120425 -.30 
13 240.00 11.742 L107 97 - 46 
13 250.00 11.044 11.100 -250 
13 260.00 10.287 10.312 - 02h 
13 270.00 9.416 9.396 222 
13 280.00 8.353 8.286 81 
13 290.00 6.863 60835 042 
13 300.00 4.557 42532 055 
20 100.00 17 544 17.311 1.34 
20 111.411 17.117 17.068 028 
20 133.33 16.497 16.493 203 
20 155.56 152761 15.797 -.23 
20 166.67 15.351 152403 = 034 
20 180.00 14.832 14.886 - 036 
20 184.11 14.677 14.716 -.27 
20 190.00 14.422 14.464 -.30 
20 20.00 13.953 14.012 = 042 
20 210.00 13.4686 13.524 Sot 
20 220.00 12.999 12.998 201 
20 230.00 12.434 12.425 07 
20 240.00 11.614 11.797 015 
20 250.00 11.065 11.100 -o31 
20 260.00 10.287 10.312 -024 
20 270.00 9o41E 92396 ee2 
20 280.00 8.240 6.286 = 046 
20 290.00 6.863 66835 042 
20 300.00 4.612 4.532 1.76 
11 190.00 13.777 14.464 “4.75 
11 195.00 13.635 14.242 -4.26 
11 200.00 13.501 14.012 -3.64 
11 205.00 13.338 13.772 -3.16 
ii 210.00 13.179 13.524 "2.56 
11 215.00 12.999 13.266 -2.02 
ii 220.00 12.80€ 12.998 -1.47 
11 225.00 12.584 12.717 -1.05 
11 230.00 12.3€2 12.425 -.50 
11 235.00 12.111 12.119 -.06 
41 240.00 11.844 11.797 039 
11 245.00 11.551 11.458 081 
141i 250.00 11.233 11.100 1.19 
44 255.00 10.881 10.719 1.51 
11 260.00 10.492 10.312 1.75 
11 265.00 10.03€ 9.873 1.65 
11 270.00 9.584 92396 2.00 
411i 275.00 9.002 8.872 1.47 
11 280.00 8.429 8.286 1.72 
11 265.00 7.700 7-619 1.06 
11 290-00 60834 6.835 -.O01 
41 295.00 5.779 52865 71.045 
11 300.00 Go29h 4.532 “5.26 
NP = 49, RMSPCT =  5E 

NP = 72, RMSPCT = 1.41 
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Table 11. Calculated P(9) critical isotherm. 


The following page gives a high-resolution examination of the critical 
isotherm of ethane as computed by equation of state (5). Column headings 


have the following interpretations-—- 


D/DC = d/d_, density reduced at the critical point. 
IS /IKG = T,(O)/T reduced coexistence temperature. 
PS/PC = P(P)/P > reduced coexistence pressure. 

P/PC = P/P > pressure reduced at the critical point. 


DP/DR = dP/d9 slope of the critical isotherm, bar.* 


The last five columns give the density-dependence of functions used 
in the equation of state, where R = 9 = d/d, is density reduced at the 


liquid triple point-- 


DIS/DR = dT,(0)/dp, K. 
DTH/DR = d6(c)/dp, K. 
DPS/DR = dP, (p)/dp, bar. 
DXB/DR = 96(p,T)/d0. 
DXC/DR = d¥(p,T)/dp. 


# 
Note: 9 = d/d.. density reduced at the liquid triple-point. 
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Table 12. Experimental and calculated P-po-T data 
The following pages compare experimental P-0-T (compressibility) data 
with densities and pressures computed by the equation of state (5). The 


first column identifies sources of the data (as in Table 6): 


ID Authors Reference 


2 Virial equation This report 
8 Reamer et al. [64] 
9 Michels et al. [54] 
10 Douslin and Harrison [18] 
16 Tomlinson [85] 
* A. K. Pal, via Pope [61] 
* Straty et al. [79] 


Data of Reamer et al. and of Tomlinson were not used for adjusting 
the equation of state. Pressure deviations for compressed liquid at low 


temperatures are not significant because (dP/dp) is extremely large. 


* 
Note: Data of A. K. Pal are identified by ID = 1201 through 2414 in uniform 
sequence. Data of Straty and Tsumura begin with ID = 3501. Their 


data at densities below 11.0 mol/% have been omitted. 
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Table 12. Experimental and calculated P-p-T data. 


EQUATION OF STATE VS. PVT DATA 


19 T,k MOL/L CALCO D,PCT Py SAR CALCG 
2 Z 350 5 MEW e405 04 20t 20.42 be 657 52695 
2 240.059 e4UUC 04301 2658 Ogee ary Mar 
2 Ae enone 0405 e4ulc -2t5 76445 (o44e 
2 ZQilo wie w o4sdl e4ucd 25 26 72614 7e8l1u 
2 270 .5Cu o4lus 4a 265) 82179 32175 
2 Elsie (Gig e4COC e4uIc =e 80 O43 32538 
2 AE GWE ee Gs vet e45ucd =o 9) 8.304 3.899 
iz 390. OCU 4905 04 302 -2u5 9.264 92259 
2 Sab GS Gis eG 5GL 4 lr e Soe) 9.622 9.617 
2 RIG o Sitl oh 08 eo42uc -e 4 20 973 90975 
2 SIGE OBES e480 e402 eis 35.5 O85 5) Lo GAs 
@ S42. 503 oH C55 e401 -.t3 AROS 15.636 
2 S56. 760 0 FE e450 4a =pne 8 11-6244 diet 4s 
2 S6f.e SEL ae rere 04324 252 alo BS) Y 11.594 
2 370.09u 04905 eGull SO we Labo DS 11.748 
2 $8L.05C eG2us °40ui OU 2 i2eiu2 12.100 
2 398.506 e4CI0 »40G1 SVepsed 122454 122.452 
@ 38. 0ff eG 8 eH Jui 26 wil 122865 £2. 853 
2 41C.cf8 Gute 0 40CF iy sal 86 ADS 13.154 
2 422,920 e4UCS e4Uue secl 136506 15.505 
2 430.2005 e4tC3 e4JCL eid 13.356 HSH WS 
2 Ghee tle 0420) o4ull ects 14.255 14.255 
2 Oo rear res Oe 042S5C 04 JOC -.209 14.555 14.554 
2 4oledgtt 04 lui aCe -eul 14.904 14.94 
2 470.CC e4CCe o%40UL Jota) L359 2 SZ 55 252 
2 48. e JLb e400 04 jlo veld 20 fs) st 15.601 
2 49L.0CC 24508 45d diene d 15.949 1529443 
2 596.505 oe JEL o4Uss 26 23) 16.297 16.297 
2 DAs Otte 040 7 04 Cue =e 16.993 16.993 
2 53h. 3eE ehudl 0459C -.539 17.345 17.345 
2 DET Oo ee 0 450F e4)UF S38 472058 17.687 
2 550.566 04 CUL 04 aut =e ut 183535 13-6534 
2 560.000 04 U7 Oe ens Soi at 13.582 130381 
2 ST as is 0405 e459 315% ob ASG FES 13.726 
2 580.9% 04 O0QT eGdce Po sd 19.576 196.7% 
2 590.29CC 04 70L e4ste Soe al 19.423 19-421 
2 68le bol eh 30E o4t St Soa 0 19.769 Sefb/ 


NER seesis  ONRM SPUN b= se 3e8, PMCANECR = «218 
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EQUATION OF STAT= VS. 
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Ow DOO POOH ODD OWOWDDOBDO DODD DOW ODD WH ODDO WDWWDDABBDOBDADOAAMWBDAAABMDABOBDWNXO 


T,K 
Salle 5 Sais) 
342.928 
Bale 6 Seis 
31%~292e6 
SLE. 928 
sjat2o Vas 
3153-928 


= abit 5 SK 


31°.2928 
310.928 
310.928 
31%. 928 
310-928 
344.261 
344.261 
344, 261 
344.261 
344.261 
344.261 
S44. 261 
344.261 
244,261 
316.928 
344,261 
377.594 
4104928 
444.261 
477,594 


510-928 


316.928 
544.261 


377-2594 


41%. 928 
44Gb, 261 
477-4594 
512.928 
Sif. 928 


Sh AG 


377-594 
410.928 
444.261 - 
477.594 
510-925 


- 3474 928- 


344.2614 
3772594 


She OUD sa 


GGG. 261 
477.594 
512.928 
ayabEo ClaZis) 
344.261 


MOL/L 
325245 
4eSuel 
525178 
7.2097 
8.5632 
8.8480 
9.1394 
92-4447 
9.6593 
9.6423 
9.99c1 

10.1321 
10.2579 
4.9251 
508407 


607547 | 


7.5579 
8.2287 
8.7434 
921754 
9.5504 
9.8579 
05889 
05165 
04627 
e41@G 
23836 
03544 
03295 
1.35466 
aoe eal Z 
29762 
08711 
07 89c 
0/233 
26084 
204961 
1.8567 
1.5546 
1.36005 
-1.2ife2 
1.1u66 
1.0165 
#22697 
2e8315 
2e213° 


1.3964. 


1.66085 
1.59548 
1.3733 
13.3736 
4.5959 


Pvt OATA 


CALCO 
325624 
4e5I2E 
526 4E2 
72-2149 
8.1555 
Be f29u 
9.106% 
9.3882 
Jo lsjats) al 
9.81" 
9.9028 
10.1054 

10.233% 
4.9022 
5.8642 
6.7262 
72494 
8.1448 
8.6913 
921453 
G.5267 
9.8525 

25928 

o Bab 7/ at 

e402c1 

04190 

25838 

03546 

23297 
1.5532 
al oA ewe 

095766 

28738 

0789 

07236 

26692 
205124 
1.86114 
1.5562 
1.3599 
1.2175 
1.1072 
1.9133 
722149 
228123 
2eci7e 

= 41.891 
1.66932 
1.5352 

1.3703 

10.3488 
4.1267 
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Py oAF 
48.263 
51.2711 
53043% 
5521556 
500037 
572316 
59.2295 
63.2574 
620 Ibs 
656432 
64-4011 
66646197 
67.569 
752 842 
82.757 
892652 
9625c7 

153.421 
110.316 
117.211 
124.16 
131.2892 
13.796 
13.79C 
13.79% 
15.796 
13.79%. 
15-790 
43.792 
27.579 
27.5979 
27-579 
27.579 
2H o BIS 
27.579 
27.579 
41.3569 
41.2569 
41,369 
41.369 
41.369 
41.369 
4i2309 
-554158 
550158 
552158 
552158 
550 1 8 
552158 
656 948 
664948 


CALCD 

+ Oe 284 

512498 

53.2353 

550223 

260950 

58.3)1 

59-6603 

6.2994 

62.372 

63-766 

65.081 

66-466 

67.854 

152553 

822557 

8S. 549 

97-1938 

164.2433 
111.046 
i117. 714 
124.2572 
131.353 
13.752 

13.77% 

12.784 

13.737 

13.787 

13.784 

13.779 

27648 

272545 

270574 

272586 

276577 

272509 

276547 

41.225 

41.296 

4126334 

41.384 

41,365 

41.548 

41.301 

550223 

552026 

55eL 77 

55.1649 

552135 

550117 


pier nsw 


09.260 
608.673 
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Wolk 
377.594 
410.928 
G4, 264 
477.594 
542.928 
310.928 
344.251 
377.594 
41.2928 
&4L, 261 


477.594 


Dale Sars 
SiS. 928 
344.2614 
377.594 
417.928 
44h, CEL 
477.594 
510.928 
31° .928 
344, 264 
377-594 


4ice 928 
444,261 
477-2594 


512.928 
Seo S23 
344.261 
3772594 


41°.9¢8 


444,261 
4772594 
op Wn ae Yate) 
310.928 
344. 261. 
377.554 
410.928 
444,261 
477.594 
519.928 
31° -928 
344.261 
377.594 


41se 928. 


444,261 
477.594 
512.928 
310-928 
344.261 


- 3772594 


410.928 
444, 261 
477.2594 


Continued. 


PvT DATA 

MOL/L CALCD 
Bo Siri 2 2.9782 
204626 204648 
Col 4¥2E 2e1 SHE 
1.9142 1.9167 
LA UBUY 467424 
11.2320 11.2977 
6.3137 6234559 
4.2652 4.U357 
Bo 23/8 $e 2462 
2075945 2el7s2 
2 oh 434i 204482 
(Aas) 202596 
sho BILAL 11.8662 
8.2287 bo1448 
5.2776 5.23591 
4,c686 4.v655 
3.4145 324245 
2.9861 2.9951 
226754 226854 
12.2182 12223845 
9.2076 9.3439 
6.4713 664317 
409226 4 e972 
%.07381 43954 
32.5569 3.5559 
321538 321653 
12.5479 12.6198 
14.1366 1u.1356 
7.5147 74914 
5.7575 5.7811 
4.7294 4.7632 
4eiL SB4 4e1G75 
3.6322 3.6452 
12.8333 12.9619 
10.6669 10.711€ 
8.3549 8.3304 
6.529% 605476 
5. S796 Se 4157. 
4.6299 4.6548 
+eid4i 41196 
43.6738 43.1465 
11.9985 66 LOSE 
9.5258 9.3157 
722289 72346 
5.9919 6.92¢2 
5.215627 521844 
4.5058 4.5839 
13.2369 1323634 
11.4646 11.5346 
9.5652. G15 7/7/51 - 
728461 72846° 
6.5583 62.5845 
526568 56888 
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P,3AR 
682948 
56506 948 
666946 
O5e J46 
632948 
862184 
502184 
862.134 
860154 
862184 
862.184 
600184 

Luce 421 
123.421 
1U3se4e1 
103.4214 
Au30e4e1 
103.421 
1use42i 
125.658 
12S e058 
129.058 
120.658 
12..058 
126.2058 
1222658 
137.895 
137.395 
137.895 
157.395 
LS 70) 31915 
137.395 
1372395 
155.132 
155.132 
155.132 
155.2i32 
1552132 
155.132 
a5 576) 13/2 
172.369 
172.369 
172.369 
172. 369 
172.369 
172.369 
172.369 
189.0u6 
189.606 
189.606 
189-5u6 
189.606 
189.606 


CALC) 
630611 
68.99) 
6828945 
68. 863 
68.775 
B56 794 
Becks 
65.887 
85.993 
386-032 
8625413 
85.2939 

1343-479 
194.493 
162.98% 
103.81 
253.162 
1432142 
103.361 
117.633 
121.189 
12u.495 
123.151 
W25e222 
125.228 
1242245 
133.916 
1372995 
138.338 
137.391 
137.259 
1372396 
137.426 
15uU.695 
153.616 
155-669 
154.741 
154.284 
154.338 
1542558 
1672436 
169.671 
172.653 
172e219 
1712545 
171.2425 
171.637 
183.2508 
1860124 
183.293 
189.639 
188.786 
183.483 
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T,K 
512.928 
315. 926é 
S440 261 
377.594 
410.928 
444,261 
477.594 
512.9238 
$1.928 
S44. 261 
377.594 
&10.928 
444.2614 
477.594 
510.928 
S12.928 
S44. 261 
377.594 
440.928 
444.261 
477.594 
oA. 928 
31.928 
344.261 
3772594 
41.928 
444. 2614 
477.594 
oye Mt ato} 
317.928 
S44. 261 
2772594 
410.928 
GGu.e 261 
477.594 
51°. 328 
310.928 
344.261 
277-594 
410.928 
KGL. 2614 
477.594 
510.928 
3124928 
344-261 
277-594 
412.922 
GG4. 261 
477.594 
Sli SIAb 
310.928 
344.261 
377.554 


Continued. 


PVT DATA 

MOL/L CALCD 
Deh S'S 50325 
13.4768 13/0 59:85 
12.7816 11.85u9 
19.9144 10.9452 
Bree 3175 8.433867 
7.0812 7eilu2s 
60.1325 6.1639 
5 eo4 446 54635 
13.816°¢ 13.8996 
12.3155 1203717 
19.7385 10.789° 
9.2546 9.2346 
8.9748 Be41SS 
6.9989 725238 
662%32 6.2677 
14.1621 14.1951 
12.7400 L257 I3S 
11-3.44 11.3655 
9.94527 9.9934 
8ef521 Bel77t 
7 of 544 Tef(3E 
629554 606973% 
14.3527 14.459 
13.0835 132.1468 
11.708C 14.8354 
10.5117 10.5092 
9.2799 9.4175 
8.4.52 8.4275 
725955 726489 
14.5077 14.6317 
13.3764 1324572 
12216004 i2e232c 
13.9866 11.0507 
959595 9-96¢~ 
3.9742 8.9974 
8.21552 Rel77’ 
14.9359 15.0328 
13.8663 13.9751 
12.8594 12.8796 
bi DICT 11.8242 
10.78c¢%4 16.8365 
9.9u49 9.9376 
9.1371 92.1406 
1502794 15-4248 
14.2981 14.4923 
13.3124 13.3984 
12.3491 12.4327 
11.4576 11-5232 
19.6337 16.6835 
9.38992 9.922e 
15.5277 1567242 
14.6722 14.7072 
13.7271 13.8326 
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Py, 3AR 
189.06 
2U6-. 343 
2062343 
20 be. 843 
22606343 
2066343 
206-343 
256. 843 
241.3515 
241.3516 
c4i.31Ft 
244.2516 
241.316 
241.516 
2412516 
2750798 
2752199 
2750196 
21501796 
2756798 
275079% 
2752798 
313-2264 
S$15.264 
3136264 
31556264 
3132.c64 
$17. 264% 
310-264 
3442738 
344.738 
5442738 
34420758 
344.2738 
344.758 
344-2738 
413.685 
413-2685 
413.645 
4132685 
413.685 
415.685 
413-2685 
4822632 
4822533 
4822635 
43822039 
432.033 
432.633 
4822633 
DLO DBA 
551.2581 
551.581 


CALCD 
1832645 
1399-414 
2)de 81d 
295.2615 
2362867 
236.115 
2052009 
clbeus7 
23223593 
2372154 
233.647 
2400014 
2422921 
249.251 
2452517 
20425708 
271.568 
c71./76l1 
2736501 
274.566 
2742815 
274.893 
2962649 
30352521 
504.923 
306.527 
306.098 
$U92tU5 
309-356 
527-261 
335.259 
3362599 
334-2311 
S41. 312 
S430 24H 
S430 S42 
387.2931 
3972958 
455-2351 
446.399 
4332378 
410.933 
41324195 
4520014 
464.707 
4740251 
4722268 
4752331 
477.553 
480.398 
519.2393 
53200078 
5332797 


3.76 
1.72 
1.11 


Geil 
Le7i 
ie to 


4239 
2017 
1.72 
Wercen 

e7C 

041 
229 
Oth 
2.75 
1.93 
1.43 
099 

243 

241 
644 
Se 8s 
coll 
1.62 
1.16 

05 

oul 
6234 
SG 77 ab 
ZaD it 
2e15d 
1.50 
1.95 

046 
5289 
3043 
3e22 
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EQUATION OF ST4Tc VSe PVT DATA 


0 Ty5K MOL/L CALCL OsPUT P, BAR CALCD P yPCT 
8 41-.9¢8 12.6375 1229346 -.2/6 551.5814 537.261 20 6 
8 444,261 12.0126 12.085¢ -e61 551.581 5412330 1.76 
8 4772594 11.2338 11.2959 SO 55 551-581 544.216 1.37 
8 516.928 1005 29% 12.5795 =e 39 551.584 546.796 087 
8 310.928 1526659 15.9910 -./9 6202523 587.282 5236 
8 S44. 261 14.9988 15.0866 5 29) 6232528 600.729 3219 
8 3772594 14.4°21 14.2067 -e74% 620.2528 6UuU e236 326 
8 41° .928 1332995 13.3021 -./8 0c 70528 642-9638 oe83 
8 444, 2614 12645U1 12.562¢ -.66 6202526 607.893 2eU4 
8 477.594 11.7416 11.8135 -.o1 6202528 61.301 1265 
8 519.928 11.6763 141.2.1205 2 4G o2beddét 614.586 ° 96 
8 Sale o Vas 16-1346 16.2321 =.69 0892476 6605819 4.16 
8 S44. 261 15.2747 15.3712 -.63 6892476 6652c51 351 
8 377.594 14.4584 14.536° mre Oe 6592476 668.593 3010 
8 410.928 13.636¢ 1327348 Ore 6892476 67u.37% Zell 
B 444,261 12.8838 12.9748 So i 6392476 6732486 2032 
8 477.594 1252548 idecoii =26$ 689.476 677.051 1.390 
8 51°. 928 14.5520 41.5966 =. 39 689.476 682.604 1.96 
176, DNRMSPCT = .509, PMEANPCT = 12071 
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TK 
273. 156 
298-142 
327.14" 
S4Re14e 
SU G5 sie 
39%.16C 
NEG SLE 
BiH Bo 158 
c€98.142 
32350140 
346.143 
Bil Bio slo) 
SoBe LEG 
&25e176 
AU So hag 
298.142 
323.145 
348.143 
BS So LS. 
396.16° 
423.217¢ 
273.150 
29&.1h42 
325014 
342,143 
3732156 
398.109 
427,176 
2982142 
322.143 
346.143 
$7 3.159 
398.160 
G22.170 
2922142 
BAS 0 BUS 
S4F.143 
SU Solas 
398.166 
422.175 
298.142 
$23. 146 
345.143 
SARS AGE 
398.16¢ 
We Genie 
ZANE N Se) 5 BUN 
3232146 
343.1453 
Bho sss 
398.166 
&23e170 


MOL/L 
08538 
0&536 


BOIS 


evuyvv 


28528 

ols Bais 

08526 

8538 
1.f672 
1.0672 
1.0672 
1.307<c 
Ye a0" 2 
1.c672 
1.8072 
1.2812 
Loe sale 
122312 
1.2&1¢ 
Le2si2 
Lo CBZ 
1.2812 
1.4376 
1.4575 
1.48772 
16487: 
1.487% 
1.437% 
1.4876 
126354 
1.6354 
41,6354 
1.6354 
1.6354 
1.6354 
1.7032 
1.7)32 
louse 
NSS 2 
1.733¢e 
tole Se 

9 FAs 
Lo PVAeD 
1.9165 


1.9165 


1.9165 
1.9165 


1.9754 


1.975% 
1.9754 
1.9754 
1.975% 
1.9754 


PvT DATA 


CALCD 
08546 
08543 
0 854. 
28538 
C8535 
OSS 
08532 
1.j686 
4.0630 
1.3676 
4.3672 
1.0668 
1.0565 
io Hie Ss 
1.2333 
Lo Zags 
1.2817 
Lo Zeit 2 
Lo@ees 
SLO BENE at 
1.28¢° 
1.4397 
1.24338 
1.4377 
1.487¢ 
41.4862 
1.4355 
41.4355 
1.25376 
1.6365 
1.6357 
1.6346 
4.26539 
1.63538 
41.7558 
1.7342 
1.71334 
15% Ze 
1.7314 
1.7814 
1.9299 
1.9179 
1.9166 
1.9152 
41.9142 
1.9145 
1.97 BE 
1.9769 
1.9757 
1.9745 
1.97338 
1.9731 
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540 581 
59.738 


CALCU 
15.862 
17.969 
ZCH Sz 
2206.04 
C+e C74 
Coeid/ 
23-046 
15.785 
21-527 
24.19% 
20.612 
292296 
312953 
342489 
Gao 825 
242749 
232057 
5122935 
S4e479 
37.627 
40 745 
236415 
274-959 
342496 
352353 
332143 
42.883 
462.5384 
29-372 
33.814 
S8elS7 
420577 
400.2758 
539.6913 
33016. 
3425632 
59.363 
43-2619 
43.2205 
520541 
320449 
372854 
43-091 
Q3ec1/7 
532264 
532249 
3300314 
382045 
44eS7Y 
43.396 
540024 
59273965 
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SQUAT LON OF SuAws WSo 


I TyK. 
298.3242 


S230 146 


e 

= 
s 

ioe 


e e e 
(oar fois Cee} (iy 


eo }¢ 


HNN FE Ot WW & ¢ 


NONFMNMONE MN. 
PRPPRPRrPHRE 
NOUF FN OWE 


aN OD ON OO ON GI CO OTN CO OG Gt 


Nf WC 
eo e 
*) ©) Ca Ce Gy ce 


42Se17¢ 
5250142 
548,143 
37S2158 
$96.16 
425.170 
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Bos Bal 6 


e 
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8. 25c4 
822.24 


NP = i161, ON=MSPCT = 


PVT DATA 

CALC O0,PUT 
cei2bi 39 Eat 
Zolgees) 25 .dh i 
Zoilees -~el2 
ZolZae evo 
ceii34 eil 
Zoreucns 6 bal 
2a S392 292 
223367 =ei2d 
223349 -e 4 
22-3329 eid 
223516 els 
223514 eid 
204135 ecli 
204111 So hat 
204598 -2 4 
Co4O7K ott 
204359 e135 
204356 eit 
200474 ~ec/ 
209455 Tectl 
229426 So Jbl 
Ze 939/60 So 2h 
So V3 U O15) 
2209372 euf 
325955 -.38 
325904 -.23 
$258 62 [0 L2 
3258356 -e.4 
325827 =e 
4.4551 -.59 
Ge 4454 Teron, 
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Table 13. Heats of vaporization and closure computation. 


LOOP CLOSURE CHECK FOR SATURATES LIQUID, 
ENTHALPY HLR VIA NBS CVs HLC VIA CLAPEYRON EQN.» J/MOL 


TK HLR HLC DIF QeVAP DH/Q,PCT 
95 563+ 5687 23 17696 230 
100 5975 6016 40 17553 023 
105 6313 6352 34% 17402 20 
110 6661 069% 33 17246 019 
ARLES 70354 7039 35 17086 20 
120 7349 1387 58 16924 023 
125 7694 (737 43 16761 026 
130 80359 3068 49 16597 30 
135 8386 B441 a4 16432 033 
140 8734 8793 59 16266 03 
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155 5784 9855 70 15761 045 
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ETHANE SPECIFIC HEATS CV(DoeT) , 
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Comparison with specific heats, C_(p,T), from 
Vv 
H. M. Roder. 
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Table 14. 


RODER 


19 
SL atic) 
119 
120 
121 
qA2 
2S 
124% 
125 


Continued. 


ETHANE SPECIFIC HEATS CV(0_T);, 


MOL/L 
14.2016 
14.604 
14.592 
14.58u 
14.578 
14.567 
14.555 
14.2543 
142532 
14.6047 
14.6037 
14.631 
142024 
14.618 
14.012 
14.2606 
15.u359 
152330 
15.035 
15.027 
15.2uié 
15eyt9 
15.001 
14.993 
14-976 
14.968 
14.959 
14.9469 
16.144 
16.129 
16.1ic1 
166113 
16.105 
16.2101 
1oeU8b6 
16.071 
16.356 
16eU41 
16.029 
16.914 
16.898 
16-864 
16.872 
10.864 
16.852 
16.841 
16.8350 
16.818 
16.807 
17.990 


TK 
efiedao 
277.105 
233.039 
2392438 
290.020 
2952963 
$01.2879 
337.773 
313.2655 
Co ebTT 
260-618 
203.689 
266.736 
209.783 
272.825 
2752-855 
2530445 
2020425 
200eC21 
204.2005 
2o7.gfl 
(Cosel 
2752845 
279.768 
233.079 
2872-579 
29106462 
2952349 
300.0625 
231.182 
2354 e160u 
237.2158 
2400143 
243.120 
2442528 
230-439 
2560503 
202.2126 
2672928 
2720240 
214.213 
213.031 
221.935 
225-787 
228 e421 
2322253 
2360004 
2592848 
2432012 
2472355 
139.118 


Py3AR 
108e+c/ 
141.461 
174.s1¢ 
206.152 
211.466 
2432438 
2742397 
Sue s52 
337280 

$3.517 
500939 
67.577 
B4e+16 
101.239 
117.360 
134.575 
T5057 
1056244 
832774 
112.764 
130.369 
1596747 
183.026 
2060147 
223.936 
2512 dbs 
274ec25 
296-522 
$23.33 
352066 
5b elb42 
8Ue7114 
103.186 
1256439 
1352831 
1792-88 
2M 
265-419 
346.735 
3374540 
292552 
63.15% 
Joe/u8 
130.375 
153.276 
130.272 
219-403 
25123uU5 
2822880 
314-275 
26.913 


107 


3/MOL/K 


CV, XPT 
48.17 
48.61 
49.212 
49.76 
49-84 
50.248 
Bel 5 eZ 
51-82 
52046 
47.20 
47230 
47.38 
7674 
47.92 
48.20 
430438 
47.15 
47.46 
47.20 
47252 
47.97 
48.29 
48.66 
49.05 
49041 
49275 
56.27 
50.62 
51.18 
4531 
45.16 
45253 
45.79 
462.06 
4608 
46.68 
47.34 
47.92 
48.53 
48.291 
44.89 
44.88 
44.97 
45216 
45231 
45.70 
45.98 
46.636 
46.75 
47.207 
44.55 


CALCU 
48.10 
482664 
49e.ec1 
49079 
49.89 
50.48 
51.09 
51.71 
P40 Bie 
47.00 
47.23 
47.47 
47.73 
47.399 
48.25 
48.53 
47.19 
47049 
47.31 
47262 
47.94 
48.28 
48.63 
49.00 
49.37 
49074 
50.13 
50-52 
51.00 
44.92 
45213 
45.235 
45.57 
45.79 
45.90 
46.34 
46.80 
47.69 
48.34 
48.70 
44.25 
44,50 
44.76 
45.02 
45.19 
45245 
45.71 
45.98 
46024 
46051 
43.60 


PCNT 
014 
-. U7 
-.18 
-.05 
-.08 
-.01 
o 06 
022 
025 
042 
e 16 
=e ag 
C2 
-213 
59 ab (0) 
-.09 
-.U8 
- 0: 
-223 
-.21 
o OE 
e 02 
004 
e1l 
208 
200 
029 
eci 
050 
e 66 
007 
039 
eo SE 
058 
e40 
72 
1.13 
C6 
© 40 
043 
1.43 
085 
e 48 
032 
025 
055 
0DE 
0 84 
1.10 
1.20 
2.13 


Table 14. 


RODER ETHANE SPECIFIC HEATS CVCUO0oT) » 


ID 
609 
610 
611 
612 
613 
614 
615 
616 
617 
507 
520 
Sen 
522 
523 
524 
525 
526 
527 
409 
414 
410 
411 
415 
416 


Continued. 


MOL/L 


17.974 
17.959 
17.950 
17.938 
17.927 
17.915 
17.903 
17.892 
172881 
18.966 
18.919 
18.901 
18.885 
13.8671 
138.857 
18-843 
16.830 
13.816 
19.657 
19-643 
192643 
19.624 
19.624 
19.6069 
19.594 
19.579 
19-565 
19.553 
192043 
19.624 
19.609 
19.594 
20-289 
20.270 
202.254 
230239 
200224 
20.209 
20-996 
20.981 
20.969 
20.957 
202944 
20.930 
20.915 


TK 
132.300 
195.613 
198.035 
201.281 
204.5195 
207-717 
210.97 
214.081 
217.237 
168.242 
166.260 
109.533 
172.795 
176.58 
179-34 
182.531 
135.736 
138.925 
145.9293 
148.039 
143.064 
151.130 
IB} al 5 al Dal 
154.270 
157.360 
100.437 
163.2445 
165.634 
148.027 
151.152 
154.c657 
157.364 
131.088 
133.939 
136-8Ui 
1392647 
1422469 
145.2790 
110.494 
112.461 
114.268 
1166167 
1182324 
120.762 
123.173 


P»BAR 
61.619 
960/76 

125-238 
1582435 
193.434 
2272542 
261.151 
294.2023 
3272615 
58 e067 
320769 
150487 
118.262 
161.393 
232/66 
245.2390 
286.077 
3262368 
15.3576 
472809 
480265 
94.575 
94.2354 
143.5%6 
191.105 
2372626 
263.740 
31306239 
472589 
G+e5ie 
143.492 
191.363 
612192 
110-.cvce7 
161.5¢2 
211.945 
261.2213 
3092415 
S0ecI9t 
742051 
114.167 
153.89 
197./398 
248-208 
296.666 
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J/MOL/K 


CV,XPT 
44.21 
44.09 
44.27 
44.50 
44.71 
44,94 
45.21 
45.246 
45.74% 
44.19 
44.15 
44.21 
44.17 
44.27 
44.43 
44.63 
44.83 
45212 
44.78 
44.89 
44.98 
GU41 
44.39 
44.48 
44.66 
44.84 
44.93 
45.10 
44.78 
44.52 
44.55 
G4.65 
45.99 
45.215 
45223 
45.225 
45-40 
45246 
44.79 
46.07 
46.20 
46021 
46228 
46041 
46.43 


CALCD 
43.78 
43296 
44.10 
44.29 
44.48 
44,66 
44.85 
45.04 
45.22 
43.50 
43e41 
43.50 
43.70 
43.86 
44.01 
44,16 
44.31 
44.47 
43.57 
43.63 
43.63 
43.73 
43.73 
43283 
43.294 
44.04 
44,15 
44.24 
43.63 
43.73 
43.83 
43.94% 
44.15 
44.20 
44,25 
44. 31 
44.37 
44.43 
45-11 
45.211 
45.11 
45e11 
45012 
45.13 
45214 


PCNT 
98 
29 
0 39 
0&7 
051 
e614 
079 
093 

1.14 

1.55 

1.67 

1.47 

1.05 
093 
° 96 

1.04 

1.14 

1.45 

2071 

2.79 

2298 

1.53 

1.49 

1647 

1.62 

1.77 

1.74 

i. ci 

2.56 

1.77 

1.62 

1.61 

3299 

2210 

2017 

2.08 

Coc? 

2028 

-e71 

2. 09 

2.37 

2037 

2051 

2076 

2.78 


Table 15. Comparison with specific heats, C whe T), from 
(33) Bier, Kunze and Maurer. 


BIER ETHANE SPECIFIC HEATS CP(P,T)» J/MOL/K 


ID P,BAR T+K MOL/L cP CALCD PCNT 
33 0.060 233.159 Uedul 50249 50.69 -.39 
33 0.000 298.156 Oe GUE 52047 52047 e 00 
33 0.000 305.235u UetOU 53.28 53.35 -.12 
33 0.000 323.2150 Uev00 55 048 55250 -.15 
33 Qeuud 348.155 de OCU 58.16 58.75 -1.02 
33 Oeuud 373.2150 Geuul 61.16 61.99 -1.36 
33 0.000 398.150 0.J0U 64.59 65024 -1.01 
33 0-000 423.150 0.900 67.78 68.47 -1.02 
33 0.060 473.1506 0.uu00 74.09 74.75 -.89 
33 Leuut 2832159 eU4LS 50.88 51.10 =-.43 
33 1.U00 298.150 e041 52077 52083 -.10 
33 1.u00 305.2350 0U40 53255 53.68 -.23 
33 1.000 323.2150 0037 55675 556 8D -.17 
33 1.000 348.158 0335 56.37 53.98 -1.06 
33 1.u0G 3730154 euse 61.37 62.19 -1.32 
33 1.000 398.150 0030 C4e71 65041 -1. 08 
33 1.000 423.150 e029 67.90 68.61 -1.05 
33 1.000 473.150 025 74.18 74.80 2 Set 
33 2eUCuU 285.150 ou 87 51.30 51-52 -o43 
33 2-000 2980156 0082 53.13 53.19 -.10 
33 2-U00U 305-2350 0080 53.95 54.202 -.13 
33 2euU00 3232150 e075 55299 56013 59 25 
33 2edUU 343.150 eu7U 58.58 59.21 -1.09 
33 2-060 3732159 065 61258 62.57 -1.28 
33 2-060 398.2150 e061 64.86 65256 -1.08 
33 2.000 423.2150 e057 68.02 68674 -1.07 
33 2.0060 4732150 e051 74227 74.96 20 9)8) 
33 5.000 283015t e223 52.80 52-96 -.3) 
33 520 0UG 298.2150 ec18 54043 54.39 e 08 
33 5.Uu0 $052350 ecu 55-09 55212 -.0o 
33 5.000 323.1590 e192 56095 57204 -.16 
33 52500 343.154 0177 59.33 59.93 -1.02 
33 5.UU0 3730150 0104 o2ce21 62296 -1.20 
33 5.000 3982015) 0153 65631 660-05 -1.13 
33 50000 423.150 0144 68.44 690615 -1408 | 
33 5.000 4732159 2123 74.63 75220 -. 84 
33 45.000 283.156 0469 56.02 56007 ~.08 
33 10.000 2980150 °438 56295 560 8% 020 
33 10.900 3052550 0+2> 57225 570 34 -.15 
33 10.000 $23.2150 0396 58.79 58.81 -.04 
33 10.0060 343.159 0363 60.68 61.26 - 298 
33 1u.ug0 3732159 0535 63436 64.02 -1.04 


33 10.000 398015 0511 66.09 06691 -1.23 


Table 15. 


BIER ETHANE SPECIFIC HEATS CPI(P ,T), 


Continued. 


P»yBAR 
10.060 
10.uu6 
2d euUUU 
293-006 
2u.u00 
20-0900 
2u.000 
2ueubG 
20.006 
20.000 
20.000 
$53 e-0UU 
30.006 
30.0060 
33.G0G 
30.uCG 
30 .UU0 
30.2000 
30.000 
32.000 
35-000 
362000 
39.0U0 
40.360 
49.u0G 
4U.2UU0 
40.000 
40.000 
40.000 
4j.U00 
46.000 
438.2800 
48.800 
48.80U 
48.80 
48.800 
48.800 
48.800 
48.806 
48e8UG 
52.006 
60.006 
60.u00 
63.000 


T,K 
423.2150 
473.156 
283.155 
298.150 
$05.2350 
3232159 
343.215) 
3732153 
398.150 
423.150 
473.1598 
29382155 
39523550 
323-150 
343-2150 
573.150 
398.150 
4232150 
4732150 
323.2153 
4732158 
3732153 
296.150 
3052356 
323.150 
3480154 
373.2150 
3982159 
423.150 
473.2150 
3050350 
307.2150 
313.2159 
316.159 
$232150 
3482158 
373.2150 
398.150 
4232155 
47321530 
3052350 
3052350 
39492159 
$13.215%6 


MOL/L 
0291 
wee 

1.0675 
0973 
0533 
e352 
0/65 
0097 
e242 
0097 
0724 

1.0589 

1.584 

1.395 

1.22138 

1.092 
©390 
0918 
0/99 

1.513 
0939 

1.347 

20D55 

2ou8u 

1.738 

1e652t 

1.373 

1.255 

1.J81 

30586 

4-362 

30403 

3.244 

20897 

2ec/il 

1.946 

1.723 

1.567 

1.336 

92288 

10-0047 

JedU0 

8e3UC 


110 


J/MOL/K 


CP 
69222 
75.18 
67045 
64.68 
63.99 
63254 
63-99 
65-91 
68205 
70.78 
76226 
81.64 
76286 
70-84% 
68250 
65292 
70.27 
(2.53 
77231 
72.98 
77-94% 
70-78 

153.57 
114.69 
84262 
74278 
72041 
72.83 
74.45 
78254 
276243 
401.61 
157.99 
134041 
108.07 
82.06 
76259 
75.36 
760141 
739257 
273.85 
162-62 
204.30 
341.87 


CALCD 
69-86 
15.77 
67-37 
64.49 
63.89 
63.57 
64.60 
66249 
08-84 
71242 
76265 
81.28 
76.46 
70.99 
69-04% 
69254 
71.09 
732156 
78.00 
72.99 
78-60 
71.70 
144.14 
112.39 
84.08 
75021 
73031 
73-70 
75.10 
79-23 
229259 
372236 
156.36 
133.4% 
107.80 
828d 
77240 
76.34 
76299 
80.30 
290.58 
169.00 
223-12 
392.90 


PCNT 
-.92 
-./S 
eil 
230 
e015 
-.04 
-.95 
-.88 
Fal) 
-.89 
-. 78 
044 
052 
-.<0 
-.80 
-.90 
-1.16 
-.87 
-.89 
-.02 
-.85 
-1.30 
6014 
2.00 
0 64 
-.57 
-1.25 
-1.19 
-.87 
-.87 
16.94% 
7228 
1.C3 
of5 
rar4s) 
-.97 
-1.06 
Sol nejel 
-1.16 
-1.00 
-6e11 
-3.92 
-9.22 
-14.93 


=THANE SFECIFIC HEATS CF(P5T), 


Table 15. 

BIER 

IO PyBAR 
aS 6ue000 
33 60.400 
33 69.-0u0 
33 bu etuL 
33 60-000 
33 00.006 
33 60-0CU 
33 6u.uU0 
33 6Ueuut 
33 60e-uUG0 
33 63-006 
33 66.000 
33 68euGu 
33 70-0060 
33 70-600 
33 70-000 
33 73200 
33 70.COC 
33 70.006 
33 70.000 
SS 70-UG0 
33 73-uUUG 
33 70.0u0 
33 70.000 
33 7uUe066 
33 73-300 
33 6u20uU0 
33 8u.U00 
33 80.0uU0 
33 8C.0G0 
33 6Ueu0U 
33 85 ul 
33 80.000 
33 80-000 
33 832400 
33 92.006 
33 10uUeuvUU 
33 100.068 
33 LuceJvil 
33 1U9-UUL 
33 106.000 
33 19u-006 
S13 156.000 


Continued. 


TsK 
3142655 
$10015d 
317.2150 
$190158 
3232150 
S$e30e155 
373.2150 
3340155 
4232155 
47320150 
3232150 
3232150 
$23615u 
3052350 
$13215) 
316.2150 
32321538 
3232153 
325215) 
3292150 
3392150 
34301545 
373.2150 
$98.159 
4232150 
473.159 
3U5035U 
323.150 
3292150 
3390150 
3480159 
$7 32150 
398.150 
4230158 
G473015u 
3480153 
$050359 
$23015¢6 
3730155 
3980150 
4232150 
4732150 


MOL/L 
Veo-rl 
6e07Y 
bece/ 
DoD25 
40083 
SeJIGs 
2eor4l 
e215 
1.935 
Loe li 
50456 
6.319 
6.38/71 
11.253 
1uetcd 
9.3531 
becb4 
7123578 
beSS1 
523939 
4.038 
4eo17 
Soi4d 
20080 
20575 
1.476 
11.564 
Jeu9r 
7844 
oellid 
50137 
3810 
3179 
cel B34 
2ec 89 
e241 
12e2ec 
162545 
7.363 
50 SLU 
&eltc 
32051 
2-932 


J/MOL/K 


cP 
27254 
458.33 
429.55 
326035 
206.76 
96-40 
32.99 
79.29 
78072 
81.07 
258.96 
296258 
2868.07 
133.00 
108.60 
198.01 
243.24 
272-70 
266206 
225240 
141.18 
115.38 
90.15 
82.96 
$1.13 
82.56 
120-61 
180.39 
198.88 
176.96 
138.89 
98.18 
87.17 
83.56 
83074 
152.33 
106.12 
126.35 
149.48 
114.27 
96043 
89.07 
36.69 


CALC!) 
44335 
416.84 
3384.67 
305645 
203-58 

97.48 

83.92 

80.22 

79264 

81. 87 
250.03 
273.15 
268.96 
136.14 
178.14 
209.31 
258.78 
262238 
2472-63 
215.12 
144.71 
117.uu 

91.17 

84.20 

62225 

63-29 
121.80 
187.15 
203.42 
170.66 
138.75 

99.57 

88.59 

65.0 

84.75 
148.65 
167.77 
130-02 
147.28 
114671 

97.07 

90.78 

87. 74 


PCNT 
-3./70 
G05 
10.59 
6.40 
1.54 
Sal oatel 
Slo abeZ 
oo ary BA 
-1.16 
-.99 
3045 
7.90 
Oe €3 
250 
-5.66 
-5.71 
6. 39 
3.78 
7. 16 
4.56 
-2.50 
-1.40 
Al tS) 
-1.49 
Slo 
-1.12 
-.98 
-3. 75 
=-2.28 
2.59 
e10 
-1.42 
-1.63 
S35 
-1.<0 
2042 
-1.55 
-2.90 
1.47 
-.39 
-1.28 
-1.92 
-1.20 


Table 16. Comparison with specific heats, C (9,T), from 
A. Furtado. P 


FURTADO ETHANE CPIP,T) AT P = Ge000 BAR 
T9K MOL/L J/MOL/K CALCD PCNT 
99-817 0.Gdd 352685 3572 037 
110.928 02039 36049 36034 041 
113.372 Qeuda 3E. 95 36.76 259 
122.039 Jetuy 37224 36297 0/3 
1332150 0-C09 37.99 37.61 1.00 
144.261 G6 .000 38.74 38.28 1.19 
155.372 0.600 39025 38.93 069 
106.463 Jevus 4det2 39.72 074 
1772594 OeGdu +uel7/ 40.51 264 
180-872 0eCdd 41.39 e1e21 043 
188.706 0.0006 441.65 41.235 ef2 
199.317 Jeuts 42.53 42.24 067 
210.928 Uetuy 43252 43.20 0/4 
222-039 GeCid 44.53 44.21 0/3 
2332150 0eCuG 45.55 45.27 60 
241.761 02-000 +6.43 %6e14 062 
244.261 Jetud 46.07 46.39 059 
2556372 0.000 47293 4756 0/6 
2660483 Ge000 +9.20 48.78 085 
2772594 0.090 50.46 50.04 e82 
232.706 JeuJd Sue 396 30.63 265 
288-706 Detdu Sie Ti 51.34% ef2 
299-817 d006 520 BY 52-667 e Si 
305-2t1 0eCUG 53221 53-34 -.23 
316.926 0.ubG 53.35 54.03 = 34 
322.039 Jeuud 55.24 55.42 -.33 
$242617 JelcS 55-61 55 eT =0 25) 
333.2150 Deddu 5ue02d 56-83 =.30 
3442261 0.000 58.12 58.25 -.22 
3556372 Je030 00.15 59.68 of7 
3662463 Uetud 61.65 61-12 085 
3662817 3-005 01.78 61.17 099 
3772594 0-090 o3.1/ 02057 0 94 
388e7U6 0.000 64.53 64.02 e 6b 
399.2617 0.u90 65.é1 65246 053 
410.926 Jetdu 66.56 660-90 ~251 
+222039 0-000 27.81 £8233 -./0 

FURTADU ETHANE CP(PsT) AT P = 172237 BAR 
TK MOL/L J/HMOLAK CALCD PCNT 
99.817 210557 08e&> 68-33 017 
110.926 20-953 o8.c9 68.63 210 
118.372 20.682 08.33 c8.-80 003 
122.59 200547 b6034 68.90 206 
133.153 230137 O69ecu 69.24 -.05 
144.261 19.722 69282 09.70 eid 
155.372 19.3014 70.70 70-33 052 
1066485 132371 71.47 71.14 240 
177.594 136430 72246 72.16 041 
180-872 18.051 73234 73-20 ecu 
183-706 17.97% C3e4l 73243 207 
199.817 17.552 74.99 74298 ef2 
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Table 16. Continued. 


210.928 17.C0G6 70.75 76.89 
222.039 16-481 738.89 79.29 
2332150 15.917 320.2 382-37 
241.761 156444 35e3u 85.45 
244061 15-299 85.94 86250 
2522-594 14.785 91.34 90.76 
2554372 214.651 D4oi1 92-51 
258.1 5t 14.6409 98e7o 95.15 
260.928 1.053 042395 71.58 
262-594 1.U37 72072 70229 
2652372 1.011 7e0.ecl 68.54 
2666483 1.6u2 09044 07295 
2772-594 e921 65.17 64.18 
2826706 890 64.05 63.23 
283-706 0857 63.28 62249 
299.6817 085 O2e4U 61.81 
3052261 0/83 o2e15D 61.671 
310.926 ef 61 62.02 61.72 
3222039 0/23 622c8 62201 
324.317 071% 62.40 62.12 
333-150 2689 62077 62256 
344.261 0659 63.05 63029 
3552372 0632 64.53 64.16 
3662483 0603 056 3U 65.1% 
366.817 e607 056.41 65017 
3772594 0585 660.18 06019 

FURTALO ETHANE CP(P,T) AT P = 28 

ToK MOL/L J/MOL/K CALCD 
2772-594 12.988 115.37 114.59 
278.706 12.873 116.13 117.07 
279.3017 120/754 118.05 119.88 

FURTADO ETHANE CP(P,T) AT P = 34 

T9K MOL/L J/MOL/K CALCD 
233.2156 16-032 81.78 81.22 
241.761 15.581 84.42 83.86 
26020463 14.2046 96.88 97.235 
2772-594 13.137 109.70 110.07 
282.706 12.620 119.26 120671 
285.928 12.238 128.82 131.31 
2872594 12.2013 136.88 139.10 
288.706 11.859 142.81 145.77 
2902372 204238 131.69 133.24 
2912433 20375 122.23 124.61 
2922594 26323 115.63 117.93 
2942261 20257 197.19 116.26 
2972059 20164 980.52 101.27 
299.817 22085 98.38 95-04% 
305.261 1.959 387.U6 67.05 
310.928 1.853 B2ecs 861.90 
3222039 1.0693 76.25 76.20 
324.817 1.059 752-61 75230 


ICIS) 


© 269 


0474 


Table 16. 


Continued. 
33321506 1.57% 73-72 
3442261 1.472 (2.22 
3552372 1.389 71.47 
3662483 1.319 71.47 
3662317 1-317 71.47 
FURTADO ETHANE CP(P,T) 
T 9K MOL/L J/NOLIK 
282.706 12.313 1ii.6u 
288-706 12.1548 129.59 
294-261 11.326 1Lol.evod 
2952928 10.994 179-256 
2972039 10727 195.01 
293-706 Se247 Ceue Bu 
299-817 30125 189.98 
299.817 $ei25 189.35 
300.650 3.049 168.60 
301.2206 3.653 157.76 
30226594 22905 137.89 
305-261 20749 1220.41 
3102926 Cope 104.68 
322.039 20213 85.2) 
3240817 20156 86.93 
FURTADO ETHANE CP(P,1) 
T9K MOL/L J/MOL/K 
296 2428 10.917 174.38 
299-317 10.682 198.79 
300.6050 10.379 221-45 
301.206 10.209 24U031 
302.035 3.9ue 2735.31 
3022594 QGeo4d $2u.83 
303.2150 9.275 401.35 
303.983 4ecl7 425225 
304.261 41738 352.26 
3052261 30922 2990 ee 
3SUS5e3/2 3.900 295012 
305.2928 32798 261.69 
306.761 3.670 223.95 
308.2156 3e5U2 184.94 
309.539 3eSt7 Lod. 3t 
310.928 32255 148.47 
310.928 3.255 146.84 
3$13-7U6 Set 74 130.85 
3162483 22930 12.28 
$222039 20769 119.02 
3240517 22626 97.38 
FURTADO ETHANE CP{P,T) 
T9K MOL/L J/MOL/K 
282.706 13.002 106.07 
233.156 12.962 1J0.69 
2886/76 126416 119.15 
2942261 41-751 135.50 
297.039 11.342 149.08 
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73232 
71.85 
71.19 
71-04 
71204 


= 441 
CALCD 
114.18 
129.92 
163.70 
186.22 
211.22 
212.03 
181.20 
181.20 
165.94 
157.93 
142.72 
124.23 
104.03 
87-64 
85-38 


= 46 
CALCD 
179.57 
264-66 
228.19 
250-70 
305.73 
$7635 
438.30 
546226 
464.43 
399259 
284-31 
276269 
245-94 
214.56 
181.93 
161.39 
147.11 
147.11 
128.37 
116.49 
162.19 
97.55 


= 49 
CALCD 
108.90 
109.54 
119.63 
137.12 
152.45 


0369 


°678 


e160 


Table 16. 


298.264 
299.817 
300-928 
$020035 
3020594 
3030150 
303.706 
304.261 
304.4817 
305.261 
305.261 
3060463 
3072039 
307.594 
3086156 
3090261 
310-928 
313.706 
3160483 
322.439 
3246817 
327.594 
333.150 


FURTADO 

T 9K 
3220039 
324.2817 
3330150 
344-2261 
3550372 
366.483 
30660817 
3772594 


FURTADO 

T sk 
310.372 
310.2058 
310.928 
J1ii1.cto 
311.463 
311.761 
312.039 
312.156 
$12.261 
3120428 
3120594 
312-706 
312.817 
3122928 
313-150 
3132428 
313.706 
313.963 
3140201 


FURTADO 
TK 


Continued. 


11.2156 156-64 162.16 
10.839 1/f7.39 1792.34 
10.595 193.13 197.29 
10.313 217.4 224.33 
10.146 C5124 243.96 
30.969 240048 270293 
Qe 4+ 204040 310.68 
3e4B1 566.10 576223 
32.133 091.97 9141.01 
Beli Ludo. du 817.38 
8.710 TASS 2 817.38 
4o9G5 553.58 713.00 
4.2604 453215 475.27 
4eeie 375256 372044 
4e248 332-39 $13-42 
4eO1o 201.83 247207 
32/768 200.05 196-72 
32433 158.14 156.03 
Secl@7 13525u 134.66 
22934 166.82 Ul Zo ake! 
20872 101.91 185249 
Coll 162.54 1u0.40 
Zora 95.37 393.219 
ETHANE CF(P,T) AT P = 51 
MOL/L JSMOLAK CALCD 
32363 HES hal 125.97 
32162 116.88 115.98 
228358 99.26 98.88 
205438 86.82 88.29 
22338 83.03 82296 
2eifod 79.27 69.03 
2oili 79213 79.97 
2.043 iisasgs 7841 
cTHANE CP(P»T) AT P = 56 
MOL/L JS/MAOLIK CALCO 
8.565 4u6.89 456084 
8-323 usGeic 4960.21 
8124 483.13 536074 
7.9ii 53U.3e 574.01 
72659 577248 662.16 
7e4E0 020055 310-85 
7232 666.51 617.46 
70143 6360205 614.49 
72055 689.46 ob1uel8 
60925 693.c2 601.93 
aed U0 O9Ge +h 449.55 
60718 693.22 591.17 
02635 631.97 584244 
00565 679. 3g 577229 
00408 044.15 361041 
6.230 595.10 538.36 
6.063 552.31 512.46 
52910 544.58 485213 
56168 485.64 457.55 
ETHANE CP(P,T) AT > = 68 
MOL/L JSMOLIK CALCD 


115 


-3.53 
-.61 
-2216 
-3.19 
-5259 
=9. 92 
~1'3'. 22 
33054 
26215 
18.79 
27281 
-28. 80 
2004 
0 83 
Sell 
be b+ 
1.66 
1.34% 
e061 
-3.03 
-3.52 
2e0ug 
2028 


efii BAR 
PCNT 
-.13 
077 
0 38 
e614 
°f9 
-.97 
-1.06 
B20 Albis) 


°4+63 BAR 

PCNT 
-11.73 
11.73 
-11.10 
-8.1i2 
“4.27 
1.55 
7240 
Jefe 
11.50 
13.17 
35027 
14.72 
15.54% 
15.03 
12.85 
9254 
7.22 
5ef2 
De 78 


e948 BAR 
PCNT 


Table 16. 


68.55 
66269 
05.63 
69.u7 
69.71 
7046 
71-21 
(ie o& 
72046 
73.23 
73.98 
(CDG SN/ 
77.030 
79227 
81.14% 
S2e41 
&€5.66 
89.46 
95.7% 
G8.89 
194.93 
111.34 
137.15 
pve mers) 
176.33 
2Iju.1E 
222.69 
Cote cl 
276.62 
277.¢u 
284.97 
285259 
2385222 
231.53 
25ee2u 
277.10 
2273.30 
240.93 
20u.55 
170.51 
138.91 
143.095 
123.29 
Lu4e4e 
94.49 
Joe 24 


ETHANE CP(P 471) 


Continued. 
110.928 21.02t 
118.37¢c 20.759 
ZAHM) 20.623 
133.2150 2522238 
144.261 19.82% 
155.372 192415 
160.483 19.Gu) 
177.594 18.577 
1862872 18.215 
188.706 13.143 
199.817 17.696 
210.926 Ib 0B 38) 
2222539 16.753 
2332150 L60e«c4i 
241-761 15.825 
2442261 15.706 
2552372 BUG hak ey 
2600483 14477 
277.594 13.757 
282.706 132389 
283276 12.9io 
299-517 11.866 
305-2261 11.213 
310.928 10.549 
314-261 92657 
316.463 92146 
31720594 8-835 
319ecel 626334 
320.372 7.980 
320.928 7e3ul 
3222039 70451 
322-7606 72248 
323./U6 02559 
3242817 620d 
324.617 02664 
SAD 00525 
3202483 Beeld 
327.594 66614 
330.372 22503 
333-156 5e1d2 
33529c8 4784 
338-7 4er2s 
S44e2cl1 G4ei27 
3552372 30594 
3660483 30 24G 
506.817 35 2 sial 

FURTADO 

TK MOL/L 
282670E 13-603 
288e/ut 13024%+ 
299.817 aa SC/ st 
305.261 SUSU 6 teh 7/ es) 
$102928 11-298 
3102483 Ag .€41 
322.039 9.866 


J/MOL/K 
94.86 
93.649 

1065.21 
114.62 
123.42 
1$+.61 
149.97 
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Table 16. Continued. 


3242817 
3272594 
33.372 
333.1556 
334.817 
337.2039 
338.706 
341.463 
344.261 
349.817 
355-372 
3062483 
366.817 
3772594 


FURTALO 

T9K 
241.761 
2442261 
2550372 
20602483 
277.594 
282.76 
288.706 
299.817 
305-261 
313.926 
3162483 
322.039 
324.817 
$27.594 
333.150 
338/706 
344.261 
3442528 
345.926 
3482150 
349.817 
3522594 
SE 5 BZ 
360.928 
3062483 
360.817 
37724594 


FURTADO 


344ec61 
347.2039 
349-317 
351.761 
3526317 
3550372 


9426 loGl.i€ 
8e9b2 108.20 
8 e456 176.13 
UoYS9 lou. 54 
7.674 161.17 
72305 13Uec9 
7.950 177.65 
Be6538 17U.23 
6eS06 161.79 
Do (fall 145.47 
ryan) 154.06 
42556 114.2% 
4eoS4L abl Sig Shs} 
4.86 1t1.66 
ETHANE CP(P,T) 
MOL/L JSMOL/K 
166b4e 79204 
152924 80.15 
152390 82.79 
14.818 35.54 
14.157 59.04 
13.892 91. 34 
13515 94.73 
12.749 1u1.66 
12.333 1U5.03 
11.856 Live 34 
11.367 115-03 
1U.818 LZCo2 
10.523 126044 
10.214 129.97 
9.552 137.89 
8.8438 144.61 
8.145 150.65 
32063 151.73 
7.342 15164%E 
7268C 151.11 
7%91 1-93.22 
7.194 144.81 
00917 14u.9i 
be4ey 132.35 
5993 126. 8d 
52969 126.44 
520509 A Glena> 
ETHANe CP(P,T) 
MOL /L J/MOL/K 
11.183 AGA STen tas 
10 .%34 LAG Bie 
9.899 docu. 68 
96344 129.76 
Geos 131.60 
30753 132074 
60590 133.12 
0535 133.23 
Be259 LSiZ5 Se 
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165.39 
174.72 
182.94 
183.61 
161.46 
170.53 
172.04 
164.36 
156.94 
142.53 
129.88 
111.9% 
112.513 
101.33 


= 103 
CALCD 
73-57 
80.10 
82.73 
Boe dL 
530645 
32-60 
95047 
102-08 
1U6.15 
111.20 
abo 
124.28 
1265.29 
132.53 
142077 
148.47 
148.77 
143.40 
147.70 
145.62 
145-72 
140.18 
136.50 
129258 
123.19 
122.82 
RZ 6 e/ 


= 120 
CALCG 
113.63 
122.53 
126.84 
13uU.20 
131.25 
131-73 
131.69 
131.62 
130.82 


-3.26 
-3.83 
-3.8/7 
Salo (8) 
36 ol{S) 
2208 
3e1o 
3245 
3-66 
400 
3-07 
2001 
2014 

032 


0421 BAR 

PCNT 
0fJ 
elo 
003 
-.67 
-.68 
-.82 
-./3 
-o41 
~04D 
-./38 
alo Sat 
Sil6 7S 
-1.46 
P3Lo 97 
-3.53 
S40 58) 
1.38 
2.19 
2050 
3063 
3268 
3e20 
3.13 
2209 
2280 
2.60 
3.93 


e058 BAR 

PCNT 
° 09 
046 (eh) 
“1.70 
-.338 
eco 
0/9 
1.U7 
1.21 
1.15 


Table 16. 


358 e15u 
300.98 
5660483 
$662817 
377.594 


FURTADO 


T 9K 
110.928 
118.372 
122.9039 
133.156 
144.261 
155.37¢ 
1666483 
177.394 
186.872 
133-706 
139.517 
210.928 
222.0389 
2332156 
241-761 
244-261 
2552372 
2662-483 
2772594 
282-706 
288.706 
299.817 
305.2261 
310.926 
322.339 
S240 817 
333.2159 
344.261 
555-372 
306-483 
3662817 
377-594 


Continued. 


7.3979 130.34 
7/33 128.3% 
tt OE I 124.04 
72240 123.84 
GekSe pa areal) 
cTHANE CRP(P,T) 
MOL/L J/MOL/K 
21.119 08.19 
23.2855 06045 
2007 31 68.56 
2023543 08.94 
19.952 O94 
139.5593 7b.u6 
19.161 70.59 
18./538 71.u9 
18.416 71.55 
132347 7i.7i 
17.926 WaCrelae 
17.499 732965 
17.057 PS OI 
16.606 77.13 
10.233 7E,03 
162124 79.0C 
15-6¢5 81.14 
15.102 83641 
14.547 So. 31 
14.2749 ot 
13.953 B9e11 
13.2314 92.97 
12.981 34.30 
12.624 37.88 
11.8€2 103.86 
11.to02 105.43 
11.033 110.45 
10.136 115.62 
Ge2uu 424u223 
Be296 120.05 
32269 118.368 
7491 112.61 
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129.54 
127.86 
123.84 
123.58 
115.87 


= 137 
CALCD 
66.23 
68435 
68-40 
68462 
68.94% 
69-38 
6996 
70.08 
71-39 
71.54% 
72057 
73277 
75215 
760074 
76014 
78258 
80.72 
83.80 
36.09 
88.19 
90.09 
94.09 
396.28 
98-76 
104.19 
102.07 
112.09 
117.406 
120.93 
120.05 
119.95 
115.59 


°62 
037 
017 
17 
1.15 


°895 BAR 

PCNT 
7. 67 
015 
ec3 
045 
0f2 
039 
289 
058 
027 
023 
edi 
Ore) 
029 
050 
062 
053 
052 
-047 
-043 
-.50 
-.43 
-1.20 
-1.5u 
-.90 
-.38 
-e22 
1.43 
7 e UZ 
-e 54 


a 
e- 


-1.33 
-2. 65 


Table 17. Comparison with speeds of sound, Wo(T), from 
(1) Poole and Aziz; (2) Vangeel; (9) Tsumura. 


SPEED OF SOUND, SAT.-LIQs. ETHANZ» METERS/SEC. 
(1)POOLE/AZIZ, (2) VAN GEEL» (3) TSUMIJRA.] 


19 ToK P,BAR MOL/L Wy XPT CALCD PCNT 
at 902.183 -050 21ecbe 1994.38 203928 2025 
9 90.350 eCudu 212-086 200629 203866 -1.58 
9 91.006 ocuu 21.056 2002.4 2033.9 2 Ab Br? 
1 92.955 eUuu 21.585 197723 2019-66 -2014 
g) 95.000 0000 21.511 1975.3 2004.3 -1.47 
al 98.723 e0Cd 21.375 1933.6 1976.1 -2.c0 
2 99-7uG 0008 212339 1934.0 1968.6 -1.79 
) 100.600 e690 212328 1939-6 1966.3 -1.38 
1 100.153 e000 21-323 LI Zit o§) 1965.1 -2025 
1 103.151 e000 21.2213 1904.5 1941.9 -1.96 
2 134.576 euue 220151 189€.5 1928.5 -1.56 
9 105-0CC ebU5 21.146 190367 1927.4 -1.25 
1 108.298 e001 21eu25 1867.2 1901.5 -1.84% 
2 109.800 0001 202.970 1862.0 1889.6 -1. 48 
.) 110.000 eCUl 202963 18674 1388.0 -1.10 
1 411.524 etulL 250 5t7 1644.24 1875.9 -1.71 
1 114.943 e002 260/81 1818.0 1848.7 -1.69 
S) 115.000 eG02 CUe779 1831.2 1648.2 =-.93 
2 115.010 else 2uel7y 1827.0 1846.1 -1.16 
i] 120-006 eb04 20.595 1794.9 1808.3 -e74 
1 121.u08 e004 202558 177426 1800.2 ~1. 44 
g 125.000 007 2be4rlt 1759.0 1763.2 -.53 
2 12526490 oul Cuess2 175228 1764.3 -.66 
1 1276247 0409 200326 173004 1750.3 -1.15 
2 130.110 el13 20.220 1721.0 1727 64% -.37 
4 132.325 2017 SOL Be 1693.8 1709.7 =. 94 
2 134.89 of22 206041 1685.90 1689.2 025 
9 135.006 e023 202U37 168529 1688.3 26 28 
.) 1402006 0633 19-849 16494 1648.0 2 C5 
9 145.000 eUo2 19.559 1612.7 1609.0 023 
2 145.086 e062 19.056 1608.0 1608.3 -.02 
1 147.123 0t75 19.579 1587.7 1592.2 -.28 
1 147.61 el Bu 19.55¢ 156324 158607 -eci 
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Table 17. 


SPEED OF SOUND, 
(1) POOLE/AZIZ;5 


I 


ONNWOUONONRFNOUONPRPONEFRPRFPUONF ONFUONFP RF ONNNNWORrFNWOrFN WO O 


T 9K 
150.000 
151.266 
154.047 
155.000 
155-9406 
159-772 
160.0060 
161-250 
161.2260 
164.860 
164.896 
165-0CC 
165.519 
168.945 
109.840 
175.006 
1742266 
174.590 
175.000 
177.140 
1792856 
180.000 
180-458 
433.076 
184.571 
184.670 
135.000 
188.832 
169.710 
194.006 
195.000 
195.090 
193.761 
199-640 
200.u00 
204042 
205.2003 
210.000 
21.351 
214.560 
215.000 


Continued. 


SATeLIQe ETHANEs METERS/SEC. 
(2) VAN GEEL, (3) TSUMURA. 
P»BAR MOL/L Ws XPT 
e097 19.469 157565 
0103 192419 156567 
e135 19-313 153607 
e146 19.270 1538.3 
e157 192240 1531.2 
e211 19.91 1494.7 
e215 “19-482 150161 
0235 19.033 1491.2 
0235 19.032 1491.4 
eo SU 13.891 1464.0 
0 3U5 16.55 146301 
03L7 1324685 146366 
0313 18.365 145302 
049090 18.729 14305 
0425 18.093 1426.5 
0429 182.087 1426.0 
e561 18.515 1391.7 
SBE 18.50¢ 1391.0 
0567 132485 138601 
2667 16.3359 13686 
0781 15.288 135009 
of 3/ 18.282 1350.3 
2808 18.265 1343.22 
0935 16.155 1325.4 
1.2015 182.93 131269 
123 18.89 131622 
1.038 18.075 ej 2q at 
1.269 17.31% 1230.3 
1.323 17.877 127661 
12347 17.865 1273-6 
1.723 17echl 123520 
1.730 17.647 1234.3 
20456 17.+37 1202.2 
20139 4170449 119961 
20174 170433 1196.0 
2et+d 17.237 1166.2 
20710 Noes 115629 
3234) 10.564 1117.4 
32382 16.369 111501 
+ e035 loe/fe 1081.5 
4 e074 16.751 1077.7 
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CALCD 
1569.5 
1559.5 
1537.8 
1530.3 
1523.0 
1493.0 
1491.3 
1481.5 
1481.4 
1453.5 
1453.2 
145264 
1448 64% 
1421.8 
1414.9 
1413.7 
1380.8 
1378.3 
137561 
1358.6 
1337.8 
133607 
1333.1 
1313.1 
1301.6 
1300.8 
1298.3 
1269.0 
1262.3 
1260.0 
1221.8 
1221.1 
1192.9 
1186.3 
1183.6 
1149.7 
1145.3 
1107.96 
1104.6 
1071.9 
1068.5 


PCNT 
2 38 
049 
e 07 
052 
054 
eil 
2 66 
065 
067 
72 
067 
e777 
e 33 
e061 
e 81 
e 86 
0/8 
©9392 
94 
73 
© 97 
1.01 
075 
0 93 
2 86 
1.17 
1.05 
2 88 
1.68 
1.06 
1.07 
1.07 
077 
1. 7 
1.04 
90 
1.00 
093 
094 
89 
285 


Table 17. 


SPEED OF SOUND, 
(1) POOLE/AZIZ, 


I T 5K 
220.006 
220.100 
2242-960 
225.u00 
229.930 
239.060 
234.910 
235.0006 
239.900 
240.006 
2442966 
244.2750 
245.060 
250.006 
250-220 
255.000 
255.2110 
259.936 
260.000 
264.970 
265-000 
269-920 
273.000 
2742170 
275.060 
285.006 
280-058 
285.006 
2852150 
288.570 
296.000 
295.000 
2952397 
296.2000 
297.060 
295-0 
299-000 
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Continued. 


SAT»LIQ. 
(2) VAN GEEL, 


Py»BAR 

4.922 

4o94G 

52887 

52895 

6.6957 

72004 

80235 

8.259 

92641 

9.679 
11.238 
11.1608 
shal ol 
13.011 
13.093 
14.904 
15.009 
17.090 
17.i2i 
19.481 
19.496 
22.057 
220161 
240460 
240951 
28-062 
28.095 
312452 
312559 
34.2055 
35 01%2 
3921535 
390493 
40.004 
402805 
41.6741 
420632 


ETHANE » 


MOL/L 
1o6.5ic 
16.588 
10.269 
16.207 
16.u16 
Loeul4 
136/58 
15.753 
15.488 
15-483 
15624 
15.216 
15.201 
14.906 
142894 
14.006 
14.593 
14.280 
14.275 
13.332 
13.930 
13.568 
13.562 
13.232 
13016% 
12.728 
12.723 
12ech4t 
12.225 
11.350 
11.080 
11.343 
10.942 
10.840 
1uU.0/79 
10.568 
1u.du6 


2 


METERS/SECe 
(3) TSUMURA. 


1037.7 
1036.4 
996.8 
997.3 
956el 
95665 
91604 
Mail Bo S 
375.3 
873.6 
83164 
83265 
331.23 
78805 
787.7 
745.0 
74406 
700.0 
70007 
65545 
65504 
€09.7 
609.0 
56926 
561.1 
511.5 
509.4 
4594 
459.0 
423.8 
404.2 
343.9 
339.0 
S31ie1 
317.8 
304.0 
28906 


CALCD 
1029.9 
1029.1 

991.4 

Heat oat 

95265 

952.0 

913.3 

912.6 

873.5 

87267 

83207 

83404 

83204 

79125 

78907 

750.0 

74900 

70507 

70501 

662.3 

66201 

618.0 

617.8 

57909 

DU AOE 

52409 

524 04 

475 04 

47309 

438.3 

423.0 

366.3 

3615 

35402 

34163 

32809 

31507 


PCNT 
07D 
0/0 
054 
©62 
058 
047 
034 
0 30 
°c0 
e010 

-.16 

-.23 

-.13 

=e 38 

-.26 

-.67 

-.60 

-.82 

-.63 

cae Wry Bs 
=1. 02 
-1.45 
-1.45 
95 (eal 
-1.98 
O25 (al 
-2.94 
- 3.48 
-32.24 
=3243 
-4,65 
-60e50 
-6.64 
-6.97 
-7e 54 
-8.21 
-9.00 


Table 18. 


TERRES ETHANE 
I TK 
2922650 
292.650 
292-658 
292.050 
2922650 
292.650 
2922656 
292.650 
323.150 
3232150 
32352150 
323-2150 
323.2150 
323.2150 
3232150 
323-150 
3232150 
323.2150 
323.2150 
323.156 
3232150 
$23.2150 
3232150 
348.156 
34382150 
348.2150 
3482150 
343-2150 
343.15u 
348.150 
348.150 
343.150 
3438-2150 
348.150 
348.150 
348.150 
346.156 
348-156 
348.158 
348.2150 
373.150 
3732156 
373.2158 
373.150 
3732156 
373.156 
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(3) Terres et al. 


SPEED OF SOUND, 


P BAR 
0981 
4.9053 
9807 
142/10 
19.613 
242517 
290420 
342323 
09861 
9.3807 
19.613 
290420 
340323 
390227 
44oi13u 
490133 
53.937 
582840 
030743 
680647 
780453 
88.260 
98.066 
0961 
9e8u7 
19.613 
292420 
390227 
490033 
382840 
082047 
73.550 
782453 
832357 
88.200 
93.163 
38.066 
iu2.970 
107.873 
117.680 
0981 
9.867 
19.613 
292420 
39-227 
49.0633 


METERS/SECe 

MOL/L Wo XPT 
0041 30628 
e2lu 30209 
0 44u 29428 
0094 285-2 
2982 27465 
1.519 263.0 
1.733 24920 
2.299 23220 
2037 324.0 
0335 312.3 
0533 29805 
1leSoi 28205 
1.069 27464 
2eu19 266.0 
20426 257.28 
20923 24965 
Sissi 24122 
4425 233-20 
5.007 227.0 
7eb%0 235020 
8-98 296-0 
92854 347.20 
10.450 39609 
0554 334.2 
0555 32528 
of &B 31464 
1.190 30302 
1.695 29202 
2ec87/ 281.2 
2.999 27le2 
32380 261.8 
4.595 255928 
4.954 260.0 
520239 265.0 
60118 272.0 
be007 286.0 
7.175 291.0 
72039 304.0 
beU59 319.90 
8.767 348.0 
euse 34561 
2328 33709 
0083 32965 
1.365 321.0 
12.491 31265 
1.357 SO04el 
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Comparison with speeds of sound, W(p,T), from 


CALCD 
303.6 
301.8 
29267 
282.8 
271.9 
259.7 
24540 
227.9 
32207 
310.3 
296.5 
280.9 
27265 
26367 
254-95 
24449 
23502 
227 02 
227 6% 
24820 
304.2 
35725 
401.3 
333.6 
32402 
$1361 
301.6 
290.0 
27864 
26765 
259 «+ 
25767 
258.8 
26361 
270.9 
231.0 
292.3 
304.7 
3182 
34723 
344.2 
33605 
327.8 
319.0 
310.5 
302.3 


Table 18. 


TERRES ETHANE 


I 
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Continued. 


ToK 
3732156 
3732150 
373.2156 
$732150 
373.156 
3732150 
3732150 
398.2150 
333-150 
3932150 
3982150 
398.2150 
3932150 
3982156 
398.156 
396.150 
398.2150 
3362150 
398.150 
398.150 
4232150 
423.150 
423e15u 
4232150 
4232150 
423-150 
+252150 
423015 
423.150 
4232156 
4232150 
4232156 
4230156 
+432150 
4482150 
4462150 
448.2156 
448.150 
448.150 
443.150 
4480156 
448.150 
448.150 
448.150 
448.150 
4482156 


SPEEDO OF SOUND, 


P»yBAR 
58.849 
2864/7 
782453 
383.269 
98.066 

107.873 
117.660 
0981 

9.807 
19.613 
296-20 
3492227 
490633 
5 be EKL 
68.647 
780453 
88.260 
38.066 

157.873 
117.680 
0981 

9.307 
19.613 
290425 
392227 
49233 
5 Bed 40 
68.647 
780453 
88.269 
398.066 

107.673 
117.686 
0961 

92607 
13.613 
290420 
39.227 
49.033 
53.845 
58 e647 
780453 
3882.26) 
98.066 

107.873 
117.686 


METERS/ScCe 
MOL/L Wy XPT 
Co+l6 29720 
$.U56 281.0 
3.70c 28608 
GUS 2852-5 
52158 28905 
52889 29902 
ce5/8 $13.0 

0J3U 355.0 
505 349.1 
e029 34205 
097K 33€.0 
1-343 32925 
1.738 323.0 
ceolobi 3172 
20015 313.0 
30100 310.7 
3e015 30923 
weld 309.9 
4elib 3128 
BoeD i 3182 
228 36502 
0285 366-8 
0585 35507 
2899 35067 
1.229 345.8 
1.575 340.8 
1.339 33603 
2e3el 33207 
20/20 33022 
30135 32869 
$2505 32901 
Geut4 330.9 
+o448 33403 
eU2c 37564 
e263 371% 
oD47 36608 
0336 36206 
1.136 35Bee 
1.448 35404 
1.770 351.2 
Col 34868 
20448 34722 
22502 3460.8 
32165 347.5 
32534 346.9 
3296 35009 
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CALCO 
29428 
288.5 
284.6 
284.5 
28905 
300.0 
314.5 
354.3 
S48 ol 
S411 
334.3 
$27.8 
$21.8 
316.3 
311.7 
308.3 
30607 
307.5 
31161 
31707 
S64e02 
35901 
35305 
34802 
343.2 
335806 
33405 
S311 
32805 
327.0 
326.9 
$28.3 
331-6 
373.8 
36907 
36502 
$6009 
357.20 
35305 
35005 
348.) 
34601 
34500 
344.8 
345.26 
34706 


°46 


e Si 
052 


0 94 


TSUMURA ETHANE SPEED OF SOUND, 


I 
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Table 19. 


ToK 
130.000 
100.060 
100.uUC6 
10cdeuuu 
100.006 
100.0006 
100.000 
100.060 
109.u00 
120.000 
120.000 
12u.JU0 
120.uU00 
120.000 
140.060 
140.006 
140.006 
140.000 
140.066 
160.000 
160.066 
100.006 
168.000 
160.000 
166.000 
185 euuG 
160.0066 
160.000 
180.000 
180.000 
2032000 
200.000 
206.000 
200-000 
200.006 
2uc.000 
200.000 
206-0900 
209-000 
200.000 
200-000 
220.000 
226.000 
2202006 
220.000 
220.000 
226-000 
22u-u00 
220-000 
240-000 
240.000 
244-0uU0 
240.000 
240.2000 


(9) Tsumura. 


P»BAR 
3190467 
2842948 
251.961 
2172283 
184.270 
149.807 
114-141 

79-623 

41.217 
3034424 
285.840 
189.126 
149-630 

81.052 
3552646 
147.077 

78.032 
284.338 
214.568 
3532276 
2812-447 
218.056 
148.464 

80.133 
3562967 
286.152 
216.001 
147.2105 

78.79% 
115.108 
141.899 
367.545 
321-402 
285.205 
2536809 
21826344 
183.259 
148.313 
115e2e9 

78.806 

492536 
155.130 
359.773 
285-346 
216.788 
147.819 

18047 
4U0.571 
100.292 
3420789 
285.2379 
216.704 
147.2328 

78067) 


METERS/SEC. 
MOL/L Wo XPT 
21./09 2051.3 
21.070 2040.0 
21.033 202E.8 
21.593 2017.3 
216555 2005-7 
21.514 1993.6 
216471 1981.0 
212429 1968e% 
21.381 1954.2 
21.106 1943.2 
c1ieJUod 1914.2 
20.575 1576.1 
262819 1800.9 
202/19 1831.0 
200451 1819-4 
20.1138 1725.2 
132396 1690.7 
202343 1788el 
20.231 175607 
19.805 1698e2 
19.€76 1662.8 
1923550 1629-9 
19.+16 15917 
19.c68 1551.5 
19.i76 1583.4 
19.0652i 1543.8 
18-862 1502.0 
18.593 1457.9 
16.511 1410.3 
462610 1435.6 
172-3534 1321.0 
18.551 147524% 
13.436 1446.9 
16.2346 1423.6 
18.cb62 1402.5 
18.16% 1377.9 
16.061 135204 
17.954 132567 
172848 129961 
17.724 1268.2 
17.017 1241.26 
17.2b7 1200.8 
17.3é1 1361.0 
17-060u 1308.1 
17 6+34 1254.4 
17.179 1194.2 
16.566 1125.4 
16-983 1147.8 
16.982 1147.5 
17.163 1242.2 
10.3957 1195.8 
16.081 113404 
16.357 1063-4 
15.974 980-1 
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CALC) 
2092.9 
2080.2 
2067.8 
205400 
2041.8 
2028.2 
2013.9 
1999.8 
1983.8 
1982.6 
1947.9 
1902.9 
1883.9 
1850.0 
1846-4 
1735.4 
1695.6 
1809.8 
177207 
1714.7 
1673.4 
16355 
1591.9 
1547.0 
1591.6 
1546.7 
1499.9 
1451.3 
1399.9 
1427.7 
1311.5 
1479.5 
1447.9 
1422.3 
1399.3 
137226 
1345.2 
1316.8 
1288.9 
1256.8 
1229.7 
1191.8 
1361.5 
1304.3 
1247 64% 
1184.9 
1115.3 
1138.2 
1137.9 
1240.9 
1191.8 
1127.9 
1055.7 

973.3 


Comparison with speeds of sound, W(0,T), from 


PCNT 
-2.03 
-1.97 
-1.92 
-1.85 
-i. 80 
-1.74 
-1.66 
-1.60 
-1.52 
-2.03 
-1.76 
-1.43 
-1.29 
-1.04 
-1.49 

-.59 

-. 29 
-1.22 

-.91 

-.97 

-.64 

-2 34 

-.02 

°c9 

-.52 

-.19 

014 
045 
074 
055 
72 

-.c8 

-.07 

09 
023 
039 
053 
067 
0/8 
90 
° 96 
075 
-.04 
°c9 
256 
078 
90 
e 84 
0 84 
eil 
33 
057 
0/2 
069 


Table 19. 


TSUMURA ETHANE SPEED OF SOUND, 


I 


OWOOWNDDONOWOWWNOWWWOWWOWWWOWWOWUOWWONOWUOUWOWOWWOWOUKOWWWOWOWWOWNWMOWHOWWUOWUOWWOWWUODNOWOWOW WTO 


TK 
240.000 
260.000 
200euU0U 
260.000 
260.000 
260.000 
26u-000 
260.000 
260.000 
260.000 
260-040 
260-00 
280.006 
280.000 
230.000 
285 -UG0 
283-000 
280.000 
280.006 
280.006 
280 .UbU 
280.000 
280.060 
280.000 
280d 0U 
280-000 
280.060 
300.000 
300.000 
300.UU9 
309.000 
300.000 
300.000 
303.0006 
309.000 
320-000 
320-000 
3202000 
326.006 
320-000 
320.000 
320200 
32u2U00 
329-008 
3202000 
320.000 
320.000 
320.006 
3230150 
323.2150 
3252150 
323.156 
3232150 
3232150 
3232150 


Continued. 


P BAR 
352593 
225.294 
3456413 
3192873 
284 er 8 
25223863 
217.659 
175.814% 
148.346 
1132645 
78.313 
31.344 
3492207 
285-846 
216.544 
147.522 
782243 
308.287 
3202351 
235.347 
25352215 
216.2837 
183.451 
147.380 
113.102 
782036 
492981 
357.883 
2652355 
251-374 
216.716 
95.448 
952463 
7323) 
692997 
103.867 
103.839 
102.914 
344.819 
274.947 
216-491 
147.067 
112.578 
99.688 
852875 
78.019 
73.477 
$10.099 
331.703 
285.364 
2510142 
2162653 
1462943 
100.339 
90.7790 


METERS/SECe 
MOL/L WyXPT 
152084 917.3 
15.537 1026.1 
10249i 1142.5 
10.38% 1120.1 
152229 1087.8 
16.678 105666 
15.836 1019.5 
15.057 97102 
152+82 93t.3 
15.236 887.6 
14.342 830.2 
14.683 786.0 
15.812 1050.9 
15-484 987.7 
15.056 907.7 
142512 806929 
13.728 67660 
15.90e 1068.5 
15.669 1023.2 
156461 967.4% 
15.293 952.0 
15.2uU58 908.2 
14.814 86369 
14.510 809.9 
14.106 75002 
13.725 67507 
13.253 598.0 
152152 97161 
14.704 89226 
14.457 85067 
14.170 803.6 
12.571 566 el 
12.571 56662 
12ei%2 509.8 
11.3383 47702 
12.00 610.9 
12.04%99 610.9 
1226454 60605 
14.2351 B77e1 
13.811 792-8 
13.225 70725 
12.196 573.8 
11.346 475963 
10.3875 43306 
10.145 37205 
90495 S28 el 
50937 297.2 
140145 B44.1 
140141 85061 
13.772 79325 
13.452 74626 
13.u73 69365 
11.983 55605 
10.561 414.6 
1v.t29 37208 
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CALCD 
9135.4 
1019.6 
1138.3 
1114.9 
1081.3 
1049.1 
1011.3 
9626 
927.9 
860.0 
624.6 
777 4 
1046.6 
981.4 
900.1 
802.9 
67365 
1064.9 
1017.8 
980.9 
944 of 
900.5 
85601 
802.6 
744.04 
673.0 
600.9 
967.9 
386.4 
843.9 
T9604 
5629 
562.9 
509.2 
47805 
604.3 
604.2 
602.3 
87263 
78604 
700.5 
567.3 
474.3 
429.26 
370.0 
32702 
298.90 
83867 
845 1 
787.5 
74HU 04 
686.7 
550.2 
410.4 
36905 


PCNT 
043 
62 


047 
e €0 
e771 
° 61 


70 


90 


1.08 
1.10 
1.02 


099 
1.13 
1.04 


-.25 


1.14 
1.02 
°89 


Table 20. Calculated P(T) isochores. 


The following pages give P(T) along isochores, as computed by the 
equation of state. The third column DP/DD is the isotherm slope (dP/9dp) 
in units of the bar and mol/&. The last two columns give the isochore 


slopes and curvatures dP/0dT, 9°P/9T-, in units of the bar and K. 


These tables show that the isochore curvatures are qualitatively 
consistent with a maximum in the specific heat C(9,T) at the critical 


point. 
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Table 20. Calculated P(T) isochores. 


P, JAR DP/Ov BP/OT CeP/OT2 
LV eo 15S 120 442 011i 48 ~eb2u23 
17.994 150457 e132 -.04017 
418.3815 142441 9 db & Sb) ~eJGu14 
19.626 152493 21209 -~.0%042 
22.436 16.347 21003 -.26010 
Gio 2aU 17.27 L993 -.b9289 
Zo Ji ate} aS) S} 029860 -.0908 
Bo SiS 19.399 el 98s ~edcuG7 
25.587 19.996 28975 -.02006 
24. 206 22885 OSS) iret -.00uU56 
252145 c1ie7o06 »C¢%66 -,030¢5 
Ce Fae 2202041 of S63 ~euscliid 
Coe 68 BS8o Babel et 959 ~2JGcC4 
270446 24.374 of£956 “e,bdLd4 
2becitc 254235 5 UES S -.JIC04 
282.971 26.088 0 958 See) S) 
29.735 26.932 ou 948 ~-Gstu3 

£6 487 272765 06945 SOG HBS 
Si. 242 28.628 0f543 -.03083 

12-996 2592468 ef 941 -.GsS02 
Seo 748 Ste 5t4 0939 ~.057062 
33.499 Sie 138 00937 -.00CS¢e 
S4e 248 31.965 ef 950 -evsdG2 
S4&29GE See79E 06934 -.05982 
356.742 SiS 2D 0i 9395 SOU eee 
Si/o GSE SA5 CU 06 SIGs wet iitlt2 
37.976 Soec9l »0929 =e Ue wee 
38.713 6.91: ei G27 ~.35GtL 
3%. 46L 37-726 ef S26 -.U5d0tL 
Go. 200 S6be541 0b 925 jeuslet 
Boe GGL 39.355 0b O24 -eJubbl 
41.67¢ 4eeilbé& 0923 seo focd 
42.517 Gi eS77 et S22 ~.Gu0Ci 
45015 412786 et 921 m~oultel 
435.891 42.595 2€920 -.C3C01 
440627 GSe4tu 20923 -~-0dCT1 
456 302 4ue255 et 919 -eJ3Gu1 
46.8351 45.81¢c C917 =~ .UuteEel 
47.305 4oeolk °C317 ~edtGol 


48.2298 472414 eb S16 e055 


Table 20. Continued. 


Wists USOC Re Wi Zod WGOLAL 


T,K P,SAS 5P/00 BP/OT AIP AG 2 
ZS OS Sie 549 Be 32e uae ji) -.bc574% 
29365 $22235 Ge 544 12534 -eb0054 
Subet 315-6 chr, LEO LAS ©2296 -20cf42 
Siders S&e 714 25 92 eccht -.3s0c9 
Sogek 4&l, 498 14.577 5 areal’) So aE 215) 
S38 0% 42.265 1561606 ec2eus = Deuce 
S4E0eFf 440015 1625¢ 2218+ -.f919 
254.2 &5.76° 172286 022609 =e Ui aer, 
$62.2 G70 49% eG SS i 22156 -.00ui6 
27 Sec 49.229 i9e376 0ci 44 -250014 
S7 Bet 566 920 Cee +7 02135 -eG3t13 
E8600 522 622 210431 0ci2s =e 02042 
59468 54,217 cece 44 22414 -.Gdd4a1 
uo Zed 560.3f5 250 46k 2e1cd =e05Ul10 
440.0 5723526 240467 22092 ~etdid 
4138.0 53.361 256468 e2eUs Spreies 15S) 
4eOeu zo BS 250465 0c. 84 ~200 58 
G34e0r 62.695 274457 2cl78 =euttS8 
442e” 64.355 250446 OBE ME =e UD a ET. 
45 Ges 6602138 C9043 22066 —eJ)087 
453.7 672.661 Steal 2 2 ce ol ~ebuELd 
456.0 69.3¢3 Bo Behe ecl5e ~eJull6 
4746o 7.295% S20 364 Os Dak —e3 505 
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Table 21. Calculated P(9) isotherms. 


The following pages gives P(p) isotherms, as computed by the equation 
of state (5). The third column DP/DD is the isotherm slope (dP/00) in units 
of the bar and mol/%. The last two columns give the isochore slopes and 


curvatures, DP/DT = (dP/dT), D2P/DT2 = (9°P/3aT7) in units of the bar and 


Kelvins. 


These tables show that dP/dp is non-negative, and that it increases 


monotonically with density to pressures about twice those used for adjusting 


the equation of state. F 
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MCU/L 
21235 
21249 
21.45 
21.53 
21.55 
21.602 
21.65 
21.75% 
21.75 
21.59 
241285 
21.95 
21.95 
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MOL/L 
21i.cs 
21205 
21.19 
21.15 
21.24 
21025 
210 38 
21635 
2104) 
21645 
2ie 50 
21.55 
21-6) 
21.65 
24.75 
Zahn (C2) 
21.80 


THE 


Wn 


Table 21. 


ISOTRERM AT 95.06 

Py ,BAR DP/uU0 
32.304 828 .12t 
74-121 8442951 
116.3i¢1 861.843 
165.316 878.777 
204.679 895.731 
249.889 912.6653 


P,BAR 
16.822 
560416 
96.f41 

138.1°7 
18).c17 
223.175 
266-984 
311.649 
Soe datas 
403.548 
456.783 
498.873 
547.615 
5972625 


P,BAR 
26.263 
626422 
99.397 

137.117 
175.652 
215.048 
2552195 
296.220 
338.096 
38%. 811i 
424,391 
468.835 
514.148 
569.3533 
607.295 
655.355 
704.155 


THe TSOTHERM AT 29t.Ud 


DF/O9 
783.592 
Bouel87 
816.904 
823.722 
B54 2.656 
867.061 
884.732 
901.849 
918.992 
936.139 
953-265 
973.341 
987.3238 

10524,220 


WSOV PER (AY bala oe 


OP/SO 
714.469 
730.356 
74650 3c4% 
762.518 
77828382 
795.409 
812.993 
828.926 
845.898 
862.939 
880.219 
897.545 
914.9€&2 
932.456 
950.978 
967.6.0 
965.216 
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Calculated P(9) isotherms. 


UEGe K 


DP/ OT 
292-8533 
BE uy Uf By 
3022997 
SU.5254 
3567528 
3Te9R2u 


DEG. K 


OP/ oT 
2804217 
28. O444 
28.8687 
22.9649 
29035229 
2525526 
2927843 
Steul/9 
3902535 
36.4912 
S$C./7369 
$U297C9 
34122171 
31.4638 
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bP/ET 
2529264 
26.1448 
260565u 
2605871 
2626110 
272369 
2702647 
2704945 
2727263 
2729652 
28.1362 
2824343 
2806747 
28-9173 
2921623 
29.4896 
2520594 


D2P/0T2 
5 dla SIE) 
-e111436 
-.109723 
-.107954 
-e1J60129 
~elG 4244 


DeP/0T2 
-2104420 
-e1032e6 
~0192U74 
~e 15.841 
-2t 99567 
-2 198249 
=e!96835 
~2i 95475 
~e v9 4517 
-.092508 
-2L 95948 
= 7910 3)9)9)3)5 
-.087669 
- 2085946 


C2P/0T2 
—-2u867u5 
2480159 
2335598 
-.085015 
~ C3844 9 
=.,6 83779 
-~2683124 
~2i 82442 
= 081733 
=. 6 85995 
-.t8S226 
-.U/S427 
-.C78596 
-ei77732 
-2076834 
SoU DE)Oe 
-2t74933 
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P,£Ag 
So eset 
35.36) 
68.246 


111.879 
1360627% 
1712427 
267.362 
24400 83 
281.-6:.2 
319.926 
359..67 
399... 52 
439.83 
481.6474 
523.966 
567.318 
611.537 
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262.541 
536135 
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108.132 
136.615 
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495.761 
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MCL/L 
19.1¢ 
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14.133 
323-2616 
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334.517 
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412.374 
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463.364 
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434.728 
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526.819 
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565.222 
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DP/OT 
etu 43 
13.2566 
13635872 
13.5191 
13.6522 
135-7866 
15.9224 
14.0594 
14.1578 
14.3374 
14.4785 
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~2 1 d5499 
-2- 25655 
-~.025820 
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=. 26259 
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P,BAR 
°813 
1.592 
62343 
13.832 
312653 
44.813 
58.315 
722168 
66.376 
150.946 
1152884 
131.196 
146.888 
162.967 
179.439 
196.311 
€132.598 
231.202 
249.394 
26729354 
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206.324 
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383.2552 
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623.2812 
€55.791 
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THI 221) 


AT 260.080 DtGe K 
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259.783 
266-61 
2732535 
285 2.5388 
287.761 
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310.518 
317.589 
3252498 
$53 04%2e 
341 .4F7 
349.636 
$58.9 23 
366.499 
375.116 
$35.376 
392.762 
451.854 
411.036 
420.399 
&29.898 
439.561 
449.383 
459.366 
469.511 
479,322 
4CGe 30 
520.947 
Dile lo? 
5224762 
5332934 
545-265 
5562819 
568.557 
569.442 
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DP/ DT 
eu dS6 
10.8812 
LEo39Us 
11.1005 
11.2117 
1423240 
12064373 
11.5517 
11.6672 
11.7858 
11.9.14 
120euci2 
Teel li 
12.c6ii 
123833 
12.5266 
12.6311 
12.7567 
1205836 
1setilo 
12.1459 
15.2714 
L5e4u S41 
Leo Ds Oak 
13.5703 
13.8059 
13.9427 
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14.2292 
14.3610 
14.5u 31 
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14.7914 
14.9376 
15-6852 
15.2343 
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15.5568 
Bol OW 
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~es03031 
et JU004 
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-e 14529 
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-2 15944 
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7215418 
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-2t18734 
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-.019215 


Table 21. 


Continued. 


THE 


ISUTHEFM AT 


P,EAR 
1.765 
3e4v6 
2UeStL 
4je1ls7 
602447 
61.676 
424.445 
128.211 
153.2138 
179.483 
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452298 
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566.¢ &9 
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658.751 
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17.3c9 
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310.858 
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356.314 
3732167 
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523.776 
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Table 22. 
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The Joule-Thomson inversion locus. 
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Table 23. Thermophysical properties of the saturated liquid. 


This table was computed along paths described in section 3.0. 


Column headings have the following interpretations-—- 


DPS/DT = dP/dt, vapor pressure, 
DDL/DT = do,/dT » saturated liquid, 
DP/DT = (dP/dT), single phase, 
DP/DD 2 (dP/d0), single phase, 
Q, VAP = avant heat of vaporization, 
CV = C(P,T), 
cs = C,(T) 
cP = G. (0,T) 
W = speed of sound 
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Table 24. Thermophysical properties along ‘isobars* 


The following pages give physical and thermodynamic properties along 


selected isobars, as computed by methods of section 3 of the text. 


The first line of each table refers to freezing liquid on the P(T) 


melting line. 


Each table at P < Po contains a blank line for the transition from 


saturated liquid to vapor, as seen by the abrupt decrease of density. 


Table headings for partial derivatives have the following inter- 


pretations-—— 


DP/DT = 9P/OT, 
DP/DD = OP/dp. 


The specific heat interpretions are-- 


WW 


CV 
CP 


C(P.T), 
Cay 


WW 


* 
These tables are extrapolated beyond the range of P-p-T data used for 


adjusting the equation of state (P v 350 bar). 
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Table 25. The Joule-Thomson inversion locus. 


THE JOULE-THOMSON INVERSION LOCUS FOR ETHANE 


T9F TsR Ps»PSI LB/CUFT ToF TyR PsyPSI LB/CUFT 
-10.00 449.67 387 28026 $1060.00 769.67 6803 22025 
0.00 4%59.67 722 28005 320-600 779.67 6909 22.09 
10.00 469.67 1045 27.85 330.00 789.67 7012 21.92 
2U.00 479.67 1356 27.0% 340.00 799.67 7111 21.76 
30-00 %89.67 1654 27 4h 350-00 809.67 7206 21.60 
40.00 499.67 1942 27624 360-00 819.67 7297 21.44 
50.00 509.67 Zens 27-04% 370.00 829.67 7386 21.28 
60.00 519.67 2485 26083 380¢00 839.67 7471 21.12 
70200 529-67 2742 26003 390¢00 849.67 7552 20-96 
80.00 539.67 2989 20 0+4 400.00 859.67 7631 20.80 
90.00 549.67 3228 26024 410.00 869.67 7706 20.65 
100.00 559.67 3457 2604 &%20.00 879.67 7778 20.49 
11U.u0 569.67 3679 25285 430-00 889.67 7847 20234 
120.00 579.67 3893 25065 44C.00 899.607 7914 20219 
130.00 589.67 4099 25465 450.00 909.67 7977 20.04% 
140.00 599.67 4298 250c7 460200 919-67 8037 19.88 
15U.00 609.67 4490 2568 470000 929-67 8094 19.73 
163.00 619.67 4675 24289 480.00 93967 8149 19.58 
170.00 629.67 4854 24011 490.00 949.67 8200 19-43 
180.00 639.67 5026 24.52 500.00 959.67 8249 19.29 
190.00 649.67 5193 24234 510.00 969.67 8295 19214 
2uuU-00 659.67 5354 24016 520¢00 979667 8336 18.99 
210.00 669.67 5510 23-98 530.00 989.67 8379 18.84 
220-00 679.67 5660 23.80 540.00 999.67 8416 18.69 
230.00 689.67 5805 23.62 550.00 1009.67 8451 18.55 
240.00 699.67 5946 23.44 560.00 1019.67 8484 18240 
250-00 709.67 6081 23027 570-00 1029.67 8513 18.25 
200-00 719.67 6212 23.10 580.00 1039.67 8540 18.11 
270-00 729.67 6339 22-92 590200 1049.67 8565 17.96 
280.00 739.67 6461 22.75 600.00 1059.67 8587 17.82 
290.00 749.67 6579 22058 610-00 1069.67 8606 17.67 
300-00 759.67 6693 22042 620-00 1079.67 8023 17.52 
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Table 26. Thermophysical properties of the saturated liquid. 


This table was computed along paths described in section 3.0. 


Column headings have the following interpretations-- 


DPS/DT = dP/dt, vapor pressure, 
DDL/DT = dp,/dT, saturated liquid, 
DP/DT = (dP/0T), single phase, 
DP/DD = (dP/d0), single phase, 
Q, VAP = a ae heat of vaporization, 
CV = C(p,T), 
(ots = C,(T) 
cP = oot) 
W = speed of sound 
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Table 27. Thermophysical properties along isobars* 


The following pages give physical and thermodynamic properties along 


selected isobars, as computed by methods of section 3 of the text. 


The first line of each table refers to freezing liquid on the P(T) 


melting line. 


Each table at P < Ba contains a blank line for the transition from 


saturated liquid to vapor, as seen by the abrupt decrease of density. 


Table headings for partial derivatives have the following 


interpretations—— 


DP/DT = dP/dt, 
DP/DD = dP/dp. 


The specific heat interpretations are--— 


CVe= C (0,7), 
CP CMPD) o 


“These tables are extrapolated beyond the range of P-p-T data used for 


adjusting the equation of state (P v 350 bar). 
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